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“Essentially, this corporation will be 
people—people of the highest quality. 
The United States Air Force recog- 
nizes that men of great scientific and 
technical competence can perform at 
their best only when they can exercise 
their initiative to the full under lead- 
ership which creates the climate for 
creativity. We expect Aerospace Cor- 
poration to provide that kind of 
environment.” 

SECRETARY OF THE AIR FORCE 


Among those providing their leader- 
ship to this new non-profit public 
service corporation are: Dr. Ivan A. 
Getting, president; Allen F. Donovan, 
senior vice president, technical; Jack 
H. Irving, vice president and general 
manager, systems research and plan- 
ning; Edward J. Barlow, vice president 
and general manager, engineering divi- 
sion; and Dr. Chalmers W. Sherwin, 
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are creating a climate conducive 


to significant scientific achievement 


vice president and general manager, 
laboratories division. 


These scientist/administrators are 
now selecting the scientists and engi- 
neers who will achieve the mission of 
Aerospace Corporation: concentrating 
the full resources of modern science 
and technology on rapidly achieving 
those advances in missile/ space systems 
indispensable to the national security. 


The functions of Aerospace Cor- 
poration include responsibility for: 
advanced systems analysis; research 
and experimentation; initial systems 
engineering; and general technical 
supervision of new systems through 
their critical phases, on behalf of the 
United States Air Force. 


Aerospace Corporation is already 
engaged in a wide variety of specific 
systems projects and research pro- 
grams—offering scientists and engi- 


neers the opportunity to exercise their 
full capabilities, on assignments of 
unusual scope, within a stimulating 
environment. 


Immediate opportunities exist for: 
© WEAPONS SYSTEM PROJECT 
DIRECTOR 

¢ SENIOR SCIENTISTS/ SUPERVISORS: 
Propulsion Systems 
Guidance Systems 
Spacecraft Design 
Telecommunications 

¢ SPACE VEHICLE SPECIALISTS: 
Senior Power Systems Engineer 
Sr. Flight Performance Analyst 
Re-entry Aerodynamicist 


Those capable of contributing in these 
and other areas are invited to direct 
their resumes to: 


Mr. James M. Benning, Room 103 
P.O. Box 95081, Los Angeles 45, Calif. 
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The American Rocket Society is pleased to announce the publication of this second volume of 
Progress in Astronautics and Rocketry. Approximately half the papers in the volume were pre- 
sented at the ARS Propellants, Combustion, and Liquid Rockets Conference held at The Ohio 


_ State University on July 18 and 19, 1960. The other papers had been submitted to the Society for 
a publication and were chosen for this volume in order to round out a coherent presentation of 
- several topics of current research interest in liquid propellant rocketry. 
The papers in this volume, listed below, contain timely and useful information for both the 
: research scientist and the development engineer. 
Simulated High-Altitude Testing of Rockets I. J. EBersTEIN and I. GuassMAN, Consideration of 


F. Massrer and E. 
Investigation of Exhaust Diffusers for Rocket Engines E. Stokes Fisuspurne and R. Eps, Detonability of — 


ALLEN SANForD Tick, and Ernest HInck, 
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With just a little bit of armor (against outside nicks and scrapes), we could ship these gyros in a 
dump truck. Nortronics-Norwood GR-H4’s protect themselves and the mission against the 
hazards of human handling. We were the first to offer 20,000 hours of life in a miniature inte- 
grating gyro (GI-K2). We pioneered the development of a fluid-damped rate gyro (GR-H3), a 
standard in the industry. We have delivered the most accurate gyros in production (10-FG). 
Statistics can only suggest the state-of-the-art approach and the quality follow-through at 
Nortronics-Norwood. Now the likelihood of hobbling a mission in handling is countered by the 


subminiature GR-H4, the nearest thing to NORTRONICS 


indestructible in a four-ounce rate gyro. A Divis 
NORTHROP 


PRECISION PRODUCTS DEPARTMENT, NORWOOD, MASSACHUSETTS » Field Offices: Highway #46, Teterboro, New Jersey, Telephone. 
ATlas 8-1750, TWX- Hasbrouck Heights 871-U + 2486 Huntington Drive, San Marino, California, Telephone: ATlantic 7-0461, TWX-Alhambra 9619-L 
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N RECENT years the rate of research and contribution 

to astrodynamic literature has been increasing rapidly. 
In each of the past 3 years we have witnessed the publication 
of about twice the number of astrodynamic papers as in the 
year preceding. The purpose of the present survey is twofold: 
First, it is to present the reader with a comprehensive listing 
of the astrodynamic papers published during 1960; second, 
it is to point out the advances in particular divisions of astro- 
dynamics. There has been a tendency in modern astrody- 
namics to generate manuscripts without a sufficient in- 
vestigation of prior work—a tendency encouraged by the 
fact that the roots of astrodynamics extend through several 
centuries. To aid the reader in the study of more recent 
literature, a listing has been compiled of what is believed to be 
about 80 per cent of the papers published in 1960 prior to 
September, plus a few of older vintage. The following tenta- 


tive subject headings have been adopted for the discussion of 
various advances: Geometry and coordinate systems, astro- 
dynamic constants, orbit determination, n-body problem, 


special perturbations, general perturbations, nongravitational 
and relativistic effects, observation theory, attitude dy- 
namics, orbit selection and transfer, and applications. 
Nineteen sixty has seen the publication of at least two 
monographs in astrodynamics that are general in nature. 
Ehricke’s book “Space Flight, Volume 1”’ (60)* as an authori- 
tative, well written and comprehensive tre: atment of astrody- 
namics; “An Introduction to Astrodynamics” by Baker and 
Makemson (9) is an introductory textbook on the subject. 
Articles of general interest are the series by Van de Kamp 
(255,256) and the first portion of an article by Battin (15). 
Other volumes published this year are compendiums of sym- 
posium papers (e.g., “Space Trajectories”) and will be dis- 
cussed by chapter under the appropriate subject division. 


Received Oct. 27, 1960. 

' This article stems in part frem research sponsored by the 
United States Air Force under Contract no. AF 49(638)-498 
monitored by the Air Force Office of Scientific Research of the 
Air Research and Development Command. 

2 Lecturer, Assistant Professor of Astronomy, and Project Of- 
ficer USAF. Member ARS. 

Numbers in parentheses indicate References at end of paper. 
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University of California 
Los Angeles, Calif. 
Ballistic Missile Div. _ 
U. S. Air Force 


Geometry and Coordinate Systems 


For the purpose of definition, the term ‘coordinate sys- 
tems’”’ is used here only to refer to a system, e.g., of orthogonal 
axes, utilized to specify the location of space vehicles and 
celestial bodies, and is not intended to signify a particular set 
of orbital parameters. This latter subject will be taken up 
aay the division heading of special perturbations. 

A study of an interesting geometrical property of confocal, 
co-planar ellipses was made by Gedeon (76) in which he de- 
termined the intersection points of these conics. Van de 
Kamp (255), in the first of a series of articles, delves into the 
geometrical properties of conics as does Battin (15). Some 
rather clever geometrical constructions involving the cycloid 
are employed by Vertregt (260) to determine the time re- 
quired for an object to move from point to point along an 
elliptical orbit. Other constructions allow for determining 
the direction and magnitude of transfer impulses. 

In a category perhaps more properly classed under the 
subject of observation theory, Murray (185) discusses the 
proper coordinate system (principal direction) for the deter- 
mination of Ephemeris Time. In this same area Schilt (231) 
studied the corrections to planetary precession and to the 
motion of the principal direction (in this case the vernal 
equinox) under the assumption that the space velocity com- 
ponent perpendicular to the plane of the galaxy is, on the 
average, independent of the galactic longitude. In chapter 
5 of Ehricke’s book (60) and chapter 4 of Baker and Makem- 
som’s (9) the subject of coordinate sy stems is" surve ved in 
considerable detail. 


Astrodynamic Constants 


The importance of vt rac constants has been recog- 
nized to a greater extent each year. Nineteen fifty-nine left 
us with markedly improved geoe aiteda constants, and the 
emphasis in 1960 has been on atmospheric constants. For the 
purpose of discussion, we shall subdivide the field of astrody- 
namic constants into heliocentric, geocentric, selenocentric, 
planetocentric and atmospheric constants. A survey of each 
of these areas can be found in chapter 5 of (9). 


The author received the first doctoral degree in astronautics awarded in this country, and is a lec- 


turer, assistant professor of astronomy at UCLA. 


An associate editor (and presently acting editor) of 


the Journal of Astronautical Sciences, Dr. Baker is the senior author of the first textbook to be 
written on astronautics, and has had over 40 papers published in the areas of space mechanics, 


attitude dynamics and rarefied aerodynamics. 


His 11 years of experience in the field began with the 


job of consultant to Douglas Aircraft Co.; at Aeronutronic Systems Diy., Ford Motor Co., he 


wee - managed three space vehicle trajectory contracts for the Air Force, Army and NASA. 


He has been 


selected as the next chairman of the ARS Space Mechanics Committee. Dr. Baker is presently 
serving a tour of active duty with the U. S. Air Force at AFBMD. 
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ments in celestial mechanics to improve all of the astronomical 
constants. In the area of heliocentric constants, he proposes 
an artificial asteroid (similar to that proposed by Wilson of 
Rand Corp. a few years back), but questions its use for a bet- 
ter estimate of the astronomical unit. Clemence feels that the 
new spectroscopic methods proposed by Brouwer and Lilley 
at the Toronto meeting of the American Astronomical Society 
last year would provide a simpler, cheaper and possibly more 
accurate estimate than would an artificial asteroid. In this 
same vein, Herrick, Westrom and Makemson (102a) have 
discussed the determination of the astronomical unit (a.u.) or 
rather the determination of laboratory units, such as the mile, 
in terms of it, and have found that a significant improvement 
in interplanetary navigation would be afforded by an im- 
proved value of the a.u. It was hoped that Doppler data 
taken from Pioneer V would improve the a.u.; unfortunately, 
as pointed out by McGuire, Morrison and Wong (162), there 
was insufficient data to do so. Their most recent analysis 
has, however, favored the Rabe over the Spencer-Jones 
value (a personal communication from Wong gives a value 
of z, = 8.797,38 see of arc). Intermixed with a basic 
argument as to the proper basis for “time,” a new unit of dis- 
tance has been proposed by Golay (82) to replace the funda- 
mental a.u. He feels that the period and wave length of 
atomic radiation can now be measured so accurately as to 
compete with the 10%:1 astronomical accuracy, and therefore 
feels that atomic wave length should be adopted as length 
units by definition, whereas other quantities, such as the 
speed of light and the a.u. should be measured in terms of 
these atomic length units. 


Geocentric Constants 


Brenner [see, for example, (29 and 30)] has felt that both 
the maldistribution of observations (all near the 70th me- 
ridian) and the use of osculating orbital elements gave rise to 
an erroneous value for Earth’s third harmonic (the pear 
shape); in fact, he doubts whether such a term truly exists. 
In rebuttal, O’Keefe (200) refutes these claims and, in par- 
ticular, shows that although the observations of the satellite 
orbit are made on Earth’s 70th meridian, Earth’s rotation 
carries this meridian around so that various portions of the 
satellite orbit come under observation. Furthermore, O’Keefe 
notes that the basic theory of Musen and Herget does not 
make use of osculating elements and that.the misleading 
short-period terms (periods less than the orbital period) have 
been eliminated just as Brenner advocates. O’Keefe’s argu- 
ment seems to have been amply confirmed by the work of 
Cohen and Anderle (42) who verified the third harmonic by a 
reduction of Transit 1B data by a method entirely different 
from that originally used by O’ Keefe, Eckels and Squires (199). 
O’Keefe in another article (198) discusses the hypothesis 
about Earth’s gravitational field made by Heiskanen and 
Vening Meinesz in light of recent satellite findings. Kopal 
(129a) analytically develops the figure of equilibrium for 
Earth and other planets in a very exhaustive manner. Murry 
(186) shows the utility of satellite observations as a new 
geodetic tool—not only to establish the figure and gravita- 
tional field of Earth more precisely, but also to relate ac- 
curately one position on Earth to another. In a field also re- 
lated to general perturbations, Sehnal (236) has discussed the 
influence of the equatorial ellipticity of Earth’s gravitational 
field on the motion of a close satellite and finds that the 
amplitude of short-period terms may be typically 20 sec of 
arc, whereas that of long-period terms may be as much as 200 
sec of arc. Sehnal also finds a third harmonic term. Finally, 
again noting reference (39), Clemence suggests that a nearly 
circular satellite be placed between 1000 and 10,000 miles of 
Earth and that it be utilized to check the time invariance of 
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Earth’s gravitational field (i.e., the difference between “gravi- 
tational’ and “atomic” time). Two satellites exhibiting dif- 
ferent mass to area quotients would be utilized to estimate 
any residual drag. 


Selenocentric Constants 


It is apparent that we will have to await the analysis of 
observational data of a selenocentric satellite before obtaining 
any great improvement in selenocentric constants. Alexan- 
drov (1) has summarized the present status of our knowledg: 
of the triaxiality (or mechanical ellipticity) of the moon and 
of other lunar constants. Few of us, however, fully appreciat 
the urgent need for improved knowledge of this apparent], 
insignificant feature of the moon. The lunar triaxiality is an 
exceedingly important detail in lunar theory and was, in fact. 
classed by Clemence (39) as the “most troublesome constant 
in the lunar theory.” It is hoped that sufficient accurate anc! 
useful data will be taken from the Air Force’s forthcoming 
Able shots (involving a selenocentric satellite) to improve our 
values for the moon’s figure. Since the lunar gravitationa! 
field is a product not only of the moon’s shape but also of her 
internal structure, MacDonald (158) has proposed severa! 
experiments to improve our knowledge of the moon’s con- 
stitution. We hope, therefore, that the 1960’s will bear the 
fruits of improved selenocentric constants based not only 
upon selenocentric satellite data, but also upon the findings 
from actual lunar landings. 


Planetocentric Constants 


In Westrom’s chapter in “Space Trajectories” (268) and in 
(9) detailed surveys are presented of the status of planetocen- 
tric constants through 1959. Ehricke’s book (60) also shows 
in chapter 3 some selected values for planetary masses, diame- 
ters, etc. (not all of which are as up to date as those given by 
Westrom). The picture of planetocentric constants is, how- 
ever, a constantly changing one, and as more and more im- 
proved data become available, a general upgrading of all of 
these constants can be expected. In particular it should be 
emphasized that the values of planetocentric constants given 
in the ‘American Ephemeris and Nautical Almanac’ each year 
do not reflect the latest revisions. On the contrary, they have 
been carried along for decades, if not for centuries, simply for 
the purpose of consistency and should be replaced by newer 
values for any accurate work. As examples of the continuing 
process of re-evaluation of planetocentric constants, which 
arises even through the use of conventional astronomical 
observations, the work of Clemence (38) and Rabe (213) is 
cited. By an analysis of the discrepancy between observed 
and computed values of the longitude of Jupiter, Clemence 
finds an improved mass of Saturn to be 1/3499.7 that of the 
sun (as compared to the older accepted value of 1/3497.64. 
[See p.100 of (9).] Rabe proposes to utilize the large per- 
turbations of 1011 Laodamia by Mars for a more accurate de- 
termination of the mass of that planet. 


Atmospheric Constants 


By far the most active area of astrodynamic constants is 
that of atmospheric constants. Jacchia (110) published a 
more definitive version of his earlier work giving quantitative 
results. A supplement to a more extensive German paper was 
published by Martin and Priester (161a). From the orbital 
elements of satellites 1958 Beta Two, 1958 Alpha and 195s 
Delta, they computed the densities in Earth’s atmosphere at 
heights of 210, 350 and 660 km. They found it possible to 
quantitatively separate the solar effect in the atmosphere from 
the diurnal fluctuations. King-Hele and Walker (122) ex- 
amined the variation of atmospheric density in two bands 
(200 to 300 km and 300 to 700 km) and noted the seasonal, 
latitudinal, monthly and diurnal effects; Cole (45) also studie: 
the theory of the propagation of disturbances from the sun t» 
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Earth, but based his work upon evidence of a proton density 
in interplanetary space lower than that previously assumed. 
Paetzold and Zschérner (204) investigated the variable con- 
ditions of the terrestrial exosphere occasioned by sun spots, 
diurnal and corpuscular influences. Of more general applica- 
tion is the work of Brandt and Chamberlain (25), in which con- 
sideration is given to the density of neutral gas in a planetary 
exosphere in which the “satellite” component of the gas is 
taken into account. Opik and Singer (202) take exception to 
this approach (it was criticism of their earlier work) and 
argue that the discussion of the bound or satellite component 
of the exosphere is premature, and furthermore that Brandt 
and Chamberlain make the invalid assumption of a Maxwel- 
lian velocity distribution in the exosphere. Brandt (26) also 
participates in an exchange with Firsoff (64) concerning the 
lunar atmosphere. Brandt computes the density of the lunar 
atmosphere to be between 10° and 10° for particles and be- 
tween 10° and 10‘ per cubic centimeter for electrons. Firsoff 
argues with Brandt and feels that Brandt oversimplifies and 
that he neglects the component of lunar atmospheric gas held 
by persorption, a portion of which is liberated each ‘“‘day” by 
the sun.- On the subject of planetary atmospheres, 
a paper by Jastrow (113a), who considers sources and 
sinks of the lunar atmosphere, and the influence of Van 
Allen-type belts on Earth’s and other planetary atmos- 
pheres, are noted. Kopal (129a) challenges the astron- 
omer and space scientist to concentrate on the improvement 
of planetary atmospheric data. He also notes that most of us 
utilize the now outdated atmospheric compositions. De 
Vaucouleurs (52b) discusses and summarizes the astronomi- 
cal, astrophysical and “geophysical” data on the planets; 
DeMarcus (50b) considers planetary interiors and the dy- 
namical features of the planets. Perhaps the most complete 
and authoritative discussion of planetocentric atmospheric 
constants is presented in a work by Urey (254b). Table 11 on 
page 413 of Urey’s paper, in particular, gives a succinct sum- 
mary of our present knowledge of planetary atmospheres. 


Orbit Determination 


In 1960 we have seen a considerable revival of interest in 
orbit determination—particularly in the area of the determina- 
tion of orbits from range and/or range rate data. Besag and 
Anderson (21) looked at this problem from the classical point 
of view; Deutsch (52a) applied the method long utilized by 
astronomers in the determination of double-star orbits from 
spectral Doppler shifts when over half of the orbit is visible; 
Guier and Weiffenbach (90) and Guier (89b) discussed a 
Doppler navigation system that promises satellite position 
accuracies to half a mile; Baker (7,9) proposed a series of 
orbit determination patterns for range and/or range rate 
data (see “Electronic” section under “Observation Theory”). 
Experience with optical data taken from the Atlantic Missile 
Range led Duke (57) to show how highly accurate orbits 
could be established there—as an example he utilized satellite 
1958 Delta One. A comparison of differential formulas, which 
yield analytical differences, with the difference between 
Gaussian orbits having slightly different initial conditions, 
was carried out by Weirauch (265). A question was raised by 
Herrick as to the errors that still remained in the Gaussian 
orbits and if, therefore, the differential expression were, in 
fact, more accurate. In the closely related problem of dif- 
ferential correction, Conte (46) proposed a new method for 
inversion of the least square matrix, whereas Herrick (100, 
101) developed a set of universal variables and associated 
analytical differential formulas for the correction of all orbits, 
whether ultimately elliptical or hyperbolic. Application of 
orbit determination to the outer planets was discussed by 
Krotkov and Dicke (138), and machine techniques utilized by 
the Russians in their minor planet orbit determination were 
discussed by Jahontova (111). Olber’s method for the de- 
termination of parabolic comet orbits has also been pro- 
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grammed on a digital computer, as was noted by Herget (98). 
Kranjec (136) proposed a slowly converging method of chords 
that circumvents the problems encountered in the orbit de- 
termination of objects on circular orbits. Batrakov and 
Kulikov in two articles (13,139) suggest a Herrick-Gibbs 
method that includes the possibility of errors in reckoning the 
time of observation—both the generation of preliminary 
orbits and a modified differential correction technique are de- 
veloped. Finally, in the application of orbit determination to 
lunar (or interplanetary) space vehicle orbits, Kizner (123) 
proposed a set of orbital parameters related to osculating ele- 
ments of the terminal hyperbola near the moon. 


The N-Body Problem 


In the past few years there have been at least four Russian 
papers concerned with the n-body problem, perhaps more if 
one considers related papers by the Russians in general and 
special perturbations. Jarov-Jarovoj (112) discusses the 
gravitational potential for two or more bodies; Merman 
(163) carries out qualitative research into the three-body 
problem; Dubosin (55) considers the stability of the regular 
motions of artificial celestial bodies; Safronov (229) examines 
the related problem of gravitational instability in plane rotat- 
ing systems having axial symmetry. A review is presented by 
Leimanis (114) of general dynamic systems, with emphasis 
upon the n-body problem. Three methods of qualitative in- 
vestigations are described by Leimanis including dimensional 
analysis, continuous induction and invariant measures. 
Rauch (216) utilized “Steffensen coefficients” that were ex- 
tended to the n-body problem by Chency, for the iterative 
solution of the n-body problem in real time. A third rather 
sophisticated mathematical formulation of the n-body prob- 
lem is carried out by Grébner (87) who utilizes the LIE series 
and reproduces the known algebraic integrals. 

In the field of equilibrium or libration point analysis, 
Thiring (252) undertakes the numerical integrations of 
small mass about the triangular libration points. For 
moderate variety of initial conditions, he finds nonperiodic 
and unstable orbits. Rabe, of Cincinnati Observatory, feels 
that it is unjustified for him to interpret these few numerical 
results as proof of the nonexistence of periodic orbits of finite 
dimensions. In a closely related paper, Sehnal (237) also 
examines the stability of the libration points in the Earth- 
moon system; he extends the original analysis of Klemperer 
and Benedikt, but like Thiiring finds that solar perturbations 
render the Lagrangian points unstable. In an investigation 
extending beyond these ordinary Lagrangian three-body libra- 
tion points, Klemperer (125,126) considers the interesting poly- 
gonal and rosette constellations of the system of masses re- 
volving in dynamical equilibrium. The equilibrium of the 
system is very sensitive to the correct specification of masses 
and distance ratios between them. (Greater latitude is, how- 
ever, available in the selection of suitable mass values than in 
the distances.) 


Special Perturbations 


Although sometimes associated with the n-body problem, 
special perturbations has the more exact signification as the 
calculation of trajectories and orbits by means of numerical 
integration. In this area astrodynamicists draw upon over a 
century of astronomical experience. Most of this experience 
has pointed up the inadequacy of the sophisticated Hamil- 
tonian approach as applied to special perturbations. This 
position has been formulated on less dogmatic grounds by 
Brillouin (31) who demonstrates the shortcomings of the 
Hamilton-Jacobi method for classical or quantized mechanics. 
Both prior experience and Brillouin’s article seem to contro- 
vert the approach suggested by Castruccio, Bass and Slotnick 
(35); nevertheless their approach may be of some advantage 
in general perturbations and will be taken up later. There 
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exist a number of so-called classical methods of special per- 
turbations; principal among these are the methods of Cowell 
and Encke, and a variety of methods classed as either varia- 
tion of parameters or variation of elements. [Cowell’s 
method seems to have been proposed anew by Knoll (127).] 


Two groups of researchers—Pines, Payne and Wolfe (210), 
and Baker, Westrom, Hilton, Gersten, Arsenault and Browne 
(11)—have attacked the problem of comparing the relative effi- 
ciencies of classical perturbation methods. For the purposes of 
comparison, the former group utilized the nonphysical analyti- 
cal solution of Euler (first proposed in 1958 by Szebehely of 
GE), and the latter group utilized the physically meaningful 
libration points of the three-body problem and an analytically 
soluble low thrust spiral trajectory. Both groups found that 
Encke’s method was superior to Cowell’s, both in point of 
economy of computing time and accuracy, and, provided that 
the parameters were properly chosen, the variation of parame- 
ters technique was shown to offer advantages over Encke’s. 
Makarov (160) formulates the classical Lagrangian equa- 
tions of the variation of elements in vectorial form and 
compares them with the systems of Herrick and Musen. The 
computation of the orbit of Pallas as perturbed by Jupiter 
(from 1955 to 1965) is utilized as an illustrative example. 

Another method of “stepping” the integration of orbital 
motion by extrapolating the parameters from one revolution 
to the next (thus something akin to the variation of parame- 
ters) is proposed by Mace and Thomas (159). An extension 
of the Encke method was accomplished by Hall, Galowicz 
and Wallman (95) who applied the law of the mean to the first 
two Encke terms and thereby generalized the standard fgx 
series to 4 series useful for any arbitrary analytical reference 
orbit. A method similar to Encke’s was proposed by Geyling 
(78) for the analysis of gravitational and radiation influences 
on a satellite. Two papers dealt with the selection of proper 
orbital parameters for special perturbations: Herrick (100,101) 
proposed a system of unified variables that are equally appli- 
cable to all conics and to all values of semimajor axis and ec- 
centricity (Herrick also develops the appropriate differential 
correction formulas); Garofalo (74) discusses a set of new 
orbital parameters that also exhibit no singularities for all 
values of eccentricity and inclination. 


Comparison of Methods 


Atmospheric Entry 


The subject of the integration of atmospheric entry trajec- 
tories attracted the interest of three investigators in 1960. 
Moe (176) considered a rectilinear approximation to an entry 
trajectory that is an extension of the Allen-Eggers approxima- 
tion and has possibilities as a reference orbit. Gedeon (77) 
looked at the determination of the characteristics of rapidly 
decaying orbits by an analytical procedure based on a power 
law variation of the atmospheric density. Baker (9,10) ex- 
tended his earlier two-dimensional work to three dimensions 
and presented a variation of parameters technique which is 
useful for entry vehicles exhibiting area to mass quotients 
that are not too large. : 


General 


In a relatively large number of papers the problems of special 
perturbations in general are considered: Lur’e (157) looks at 
the perturbative change in Keplerian elements; Sperling (245) 
utilizes formulas good for small values of the eccentricity [es- 
sentially a special case of the already cited work of Herrick 
(100,101)]; Brenner, Fulton and Sherman (29) discuss short- 
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period terms, such as in the “osculating eccentricity” that 
might be overlooked by an improper integration technique (see 
section ‘‘Geocentric Constants’ under ‘‘Astrodynamic Con- 
stants’); and Kozai and Whitney (133) utilize numerical in- 
tegration (based on variation of parameters) to predict the 
influence of solar and lunar perturbations on satellite 1959 
Delta Two. In this last paper attention was called to the in- 
fluence of Juni-solar perturbations on the premature demisc 
of satellite 1959 Delta Two. With regard to the application 
of special perturbations to interplanetary trajectory compu- 
tation, Merrilees and Walker (164a) found that the numerica| 
integration of the planets along with the vehicle was just a- 
fast and accurate as was the interpolation of ephemerides 0 
the planets. Conte of STL disputed this conclusion an: 
offered in evidence years of study of this very same problen 
carried out by himself and his associates. 

Separate from the special perturbation method (which deal- 
with the form of the differential equation) is the subject o: 
integration technique. The latter problem is the subject 0: 
two Russian and one American paper. Myachin (190 
presented results of an application of a general theory for th: 
estimation of numerical integration errors to the differentia! 
equations of the undisturbed motion of celestial bodies; 
So¢ilina (243) gave results of the applications of his estima- 
tion method for the accumulation of errors in numerical in- 
tegration in some problems of celestial mechanics; Schlesin- 
ger (233) considered the accumulation of errors in difference 
techniques (including the second-sum process), Runge-Kutts 
and predictor-corrector techniques. Schlesinger concluded 
that the second-sum process of Herrick is best when the in- 
tegration is stable; the Adams-Bashforth method is best if 
the integration is unstable (in this latter case four differences 
are recommended if highly unstable and seven differences if 
moderately unstable). 


Integration Techniques 


General Perturbations 


General perturbations differs from special perturbations in 
that it utilizes an analytical series development to compute the 
orbit (i.e., to generate an ephemeris) of a space vehicle. In 
this area the use of Hamiltonian mechanics appears justified, 
but the arguments of Brillouin (31) show that it is still often 
not the preferable method. 

mA 

A number of very significant theoretical advances in general 
perturbations was achieved in 1960. Perhaps the most sig- 
nificant publication was the doctoral dissertation of Kovalev- 
sky (132) who, under the direction of Brouwer of Yale, ex- 
tended and applied the Delaunay method to the VITIth 
satellite of Jupiter. Of equal importance, and perhaps more 
direct application to satellite orbits, is the exceedingly clever 
approach of Kozai as exemplified by one of his latest papers 
(134) in which he derives osculating elements from mean ele- 
ments by adding short-period perturbations. Vinti (261) 
makes good use of Hamiltonian mechanics in his new method 
for solution of an unretarded satellite orbit. As mentioned 
previously, Castruccio, Bass and Slotnick (35) propose some 
rather sophisticated mathematical techniques for astrody- 
namics that may have application in general perturbations. 
Two papers proposed the use of lunar theory in the analyses 
of satellite motion: Musen in (187) [and in (188a) as re- 
ported by Siry of NASA] discusses the programming of 
Hansen’s method of general perturbations on electronic digital! 
computers, and Lanzano (142) solves Hill’s equations of lunar 
theory to obtain the coordinates of the periodic trajectory 
as a Fourier series of the time with respect to a rotating 
reference system. 
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Aspherical Earth 
Still the most popular subject for general perturbations 
papers is the study of satellite motion about an aspherical 
Earth—1960 has seen at least eight such papers. Hall and 
Galowicz (94) carried out a check of their method which 
neglects terms of order J? and e4J? (where J is the coefficient 
of the second harmonic of Earth’s gravitational field) by a 
special perturbations integration and found a residual of only 
0.03 nautical miles after 60 revolutions. Hori (106) avoids 
the difficulties encountered near the critical inclination by ex- 
panding in terms of /J instead of J. A modification is 
made of King-Hele’s classical work (121) by Message (164b) 
. in which first the secular motion of the node is obtained in 
terms of the more conventional mean osculating descriptive 
orbital elements (similar to Brouwer’s treatment of the con- 
onical equations) and second in the singularities in the 
Je terms are removed. Ewart (62) gives a “direct’’ and 
then a ‘‘variation of elements” approach, and justifies the 
neglect of all terms higher than J in Earth’s gravitational 
potential because of a possible ambiguity owing to at- 
mospheric drag. Anthony and Fosdick (8) find a solution for 
arbitrary eccentricity [somewhat similar to Garfinkel’s (73) | 
and specialize to polar and equatorial orbits. Orbits of arbi- 
trary inclination are obtained through an auxiliary formula. 
Brenner (30), as previously noted, discussed the inherent 
error and variation in orbital elements during a satellite circuit 
and recommends a method that circumvents elliptical ele- 
ments. The Delaunay elements are utilized by Sehnal (236), 
and the second and third harmonic terms are included (see 
the astrodynamic constants section). In three papers, by 
Glitzer and King-Hele (81) and by Grebenikov (85,86), the 
secular perturbations in aspherical Earth satellite motion are 
discussed. Sehnal (238) considers the perturbations that act 
upon the orbit of a 24-hr satellite. 


Luni-Solar Perturbations 


In 1960 particular attention was given to luni-solar pertur- 
bations. As pointed out in a preceding section, Kozai 
and Whitney (133) applied special perturbations to this prob- 
lem, whereas Upton, Bailie and Musen (254a) who qualita- 
tively develop a series, Blitzer (22) who shows that the princi- 
pal effect is an orbital precession about the plane of the eclip- 
tic, Moe (177) and Penzo (206) who solve the problem 
by holding the sun and moon fixed, all attack the problem 
by general perturbations. New short periodic and secular 
solutions of the problem of the moon and satellites were dis- 
cussed by Heinrich (97). 


Satellite Drag ty 


A number of researchers examined the problem of the treat- 
ment of satellite drag. Of interest here are the researches of 
Geyling (79), Baker and Makemson [(9), chap. 9] and papers 
by four Russians: El’yasberg (61), Batrakov and Proskurin 
(14), and Ljah (153). Comments on the first two Russian 
papers may be found in the reviews included with their trans- 
lations in ARS Journat Russian Supplement. 


Nongravitational and Relativistic Effects 


As is the case with certain of the other divisions of astrody- 
namics, it proves convenient to subdivide nongravitational 
and relativistic effects into the following categories: Transi- 
tional drag and sputtering, low thrust forces, electromagnetic 
forces, radiation effects, and relativistic effects. [For a 
definition of terms in the following, see the glossary of terms 
beginning on p. 285 of (9).] 

i 


Transitional Drag and Sputtering 


Although some might consider it an unwarranted trespass 
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into the domain of aerodynamics, the astrodynamicist has 
found it necessary to interest himself in aerodynamic drag at 
high altitudes. The subject of free molecule flow has been 
considered by Schrello (234) through a novel approximation 
that is not restricted to specific values of the accommodation 
coefficient, but does require that the ratio of the vehicle’s 
speed to that of the mean molecular speed be greater than 2. 
Typical satellites having perigee altitudes between 200 and 
400 km were found by Charwat (37b) to spend 5 to 10 min 
within the range of validity of the single collision theory pro- 
posed by Charwat and Baker in 1958. Charwat extended 
this theory and considers the transverse momentum balance 
on asymmetrical surfaces: Cohen (44) restricted himself to 
an accommodation coefficient of unity, but considered free 
molecular flow over nonconvex bodies. (Cohen’s may well be 
the ultimate approach.) This same problem is dealt with by 
Chahine (36) who considered interreflection with various 
degrees of accommodation and utilized illustrative examples 
of an infinite cylinder and a spherical surface. Schaaf, 
Maslach, Maulic and Chahine joined forces in a discussion 
of aerodynamic stability at ultra-high altitude (230); whereas 
Patterson (205b) considered recent trends in the mechanics of 
highly rarefied gases from the Boltzmann equation view- 
point. Actual experiments involving slip flow have been 
accomplished by Nagamatsu, Sheer and Schmid (191) over a 
Knudsen number variation from 0.4 to 86 and a Mach number 
variation from 8 to 25. They found that slip flow increased 
with increasing Mach and Knudsen number in the vicinity of - 
the leading edge of a flat plate. Masson, Morris, Bloxson, 
Rhodes and Charwat (161b) reported on the experimental de- 
termination of sphere drag in the transition between con- 
tinuum and free molecule flow. Sputtering was considered in 
a portion of an excellent work by Benedikt (18) in which a 
number of disintegration barriers to extremely high speed 
flight were discussed. 


Low Thrust Forces =? 


At least four papers considered the effect of low thrust — 
forces on space vehicles. Rider (220) extended the work of 
Lorell and Lass (ARS preprint 872-59) to the modification of — 
the orbital inclination and the longitude of the ascending 
node. The modification of these orientation elements was ac- 
complished by reversing a normally directed low thrust device — 
twice every revolution. Arthur, Karrenberg and Stark (5) 
utilized an approximation based upon ballistic reference orbits 
for the feasibility or design approach to low thrust trajectories. 
The variation of parameters scheme proposed in (11) was 
utilized by Brown and Nelson (32) in the analysis of thrust 
orientation patterns. Their conclusions can be summarized 
in a table (see Table 1), and can be obtained by an inspection 
of the variation of parameters formulas themselves. 

Finally, Szebehely (249) considered the inverse problem: 
Given the gravitational force field and the desired orbit as 
initial conditions, find the requisite thrust force. As an ex-— 
ample he showed that by allowing the speed to change along 
a circular orbit, nonradial forces were admissible. 


Wilson and La Paz on the subject of the magnetic damping of 


the rotation of the Vanguard I satellite. Wilson in his original 
article (270a) [see also (270b)] found an agreement between 
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observations and theory for the retardation of tumbling oc- 
casioned by electromagnetic couples on the conducting and 
magnetic parts of the satellite. La Paz (143) disapproved of 
Wilson’s numerical approach and showed that perfectly good 
analytical integrals existed (which La Paz had published some 
years back); furthermore, La Paz found an error in Wilson’s 
work. In his rebuttal (269) Wilson conceded that La Paz had 
a good point, but that, remarkably enough, errors had also 
entered via the observational data that just canceled his other 
errors. Thus, happily, agreement was again achieved be- 
tween theory and observation. Two papers were published 
by Rand that dealt with the wake of a satellite traversing the 
ionosphere. In the first paper (214), a particle treatment was 
utilized and two conical regions in space were found behind a 
disk; examples of the disk (both conductor and dielectric) 
and a wire were given. In the second paper (215), considera- 
tion is given to the damping of the plasma oscillations that 
constitute the wake of a line charge moving supersonically 
through a low density plasma. Chang and Smith (37a) 
published their work on charged satellite drag that was 
initiated many years ago at Rand Corp. By imposing an 
equilibrium condition between positive ion and negative elec- 
tron currents (neglecting any sheath surrounding the satellite 
and any diminution of electron concentration in the wake), 
they found a very low voltage, on the order of a fraction of a 
volt, for most satellites. Hewitt (104) analyzed the effect of 
electrical charge on the vertical descent of a missile and con- 
cluded that it is entirely negligible. In three 1960 research 
papers charge and magnetic field interactions with satellites 
were investigated. Beard and Johnson (16) in a very clear 
article found an analytical expression for satellite charging 
that, like that of Chang and Smith, yielded a satellite poten- 
tial of a fraction of a volt. Kraus and Yoshihara (137) ex- 

amined some aspects of the coupling between the fluid motion 

and the electrostatic properties of the plasma assuming the 

boundary layer about the satellite was absent. They noted 

that a complete fourth-order system was required for a pre- 

cise analysis. Wyatt (274) considered the satellite as a com- 

pletely permeable spherical shell of charge, thus avoiding 

some difficult boundary conditions that are introduced by 

the usual exact linearized treatment. The effects of per- 

meability were shown to be approximately removed by an 

iterative process, and an approximate result was obtained 

for satellite drag arising only from electrical effects. Wester- 

man (267) extended the proposal made by Fain and Greer last 

year concerning the modification of satellite orbits by arti- 

ficially charging a satellite and allowing it to interact with 
Earth’s magnetic field (as a matter of fact he found a small 
scale error in one of Fain and Greer’s charts). Essentially, 
Westerman utilized an analytical perturbational approach 
based upon weak magnetic field forces and circular orbits. 


In 1960 a number of thought-provoking papers were about 
radiation pressure. Musen (188a), and Musen, Bryant and 
Bailie (189) agreed with Parkinson, Jones and Shapiro (205a) 
on the efficacy of solar radiation pressure on satellites that 
exhibit large area to mass quotients. In fact, even for Van- 
guard I, Musen et al. (189) found that the radiation influence 
was comparable to that of the sun. Parkinson et al. (205a) 
noted that at an altitude of 1000 miles, solar radiation pressure 
could displace the 100-ft Echo balloon satellite some 3.7 miles 
per day. Geyling (78) utilized his Encke-type method (al- 
ready noted) in the analysis of gravitational and radiation 
pressure perturbations. On the related subject of solar sail- 
ing, London (155a) found two classes of closed analytical 
solutions to the problem—the first yielding a logarithmic 
spiral and the second, conic sections. He analyzed the ap- 
propriate trajectories for orbital transfer and found that 
neither the spiral nor the conic are optimum. 
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Relativistic Effects 


One section of Benedict’s paper, which concerned the 


scientific significance of deep space exploration (19), dealt in 
a very clear cut fashion with the interpretation of the rela- 
tivistic effect on space vehicles, and essentially translated the 
somewhat formidable theorems of mechanical equivalence 
developed by Levi-Civita into relatively familiar terms. 
Benedict concluded that definitive orbits were necessary in 
order to detect these second-order relativistic effects. In two 
articles Bogorodskii (23,24) considered the relativistic per- 
turbations in the motion of artificial satellites. Utilizing 
Einstein’s field equations, including the influence of Earth's 


rotation, he investigated the equations of motion of a satellit« . 


on the basis of the method of variation of elements. The 
general relativistic deflection of light in a gravitational fiel: 
was the subject of two publications: Mikhailov (175) and 
Lillestrand (151). Lillestrand developed the magnitude o/ 
the light bending in the neighborhood of the sun and the: 
proposed to utilize a “stellar aberrascope” to measure thi 
deflection accurately (see “Optical” section under ‘‘Observa 
tion Theory’’). Ina propulsion paper, Huth (108) considere: 
the relativistic theory of rocket flight for very high speeds- 
an area which may well be of concern to the astrodynamicist~ 
of the future. 


Observation Theory 


Observation theory deals with the design and management 
of instruments and systems of instruments, and the reduction 
of data obtained therefrom which are utilized for the acquisi- 
tion, orbit determination or guidance of space vehicles and/or 
for the collection and treatment of data obtained from space 
vehicles. In this sense observation theory has much in com- 
mon with the subject of astrometry. Again it proves useful 
to subdivide the field, in this case into subsections concerned 
with time, space range systems (i.e., the management of ob- 
servational instruments), optical observations, and _ elec- 
trical observations, as was done in (9). 


Time 


Time is perhaps the most important coordinate of an ob- 
servation. In 1960 a new standard uniform time, called 
Ephemeris Time (ET), was introduced as argument into the 
tables of lunar and planetary position (ephemerides). Gab- 
bard (72) comprehensively discussed the historical background 
of, requirement for and determination of, Ephemeris Time. 
He described the Markowitz dual-rate camera which is 
utilized to ‘tell time” from the moon, and atomic clocks 
which utilize atomic vibrations for their ‘‘pendulums.”” Mur- 
ray (185) corrected his earlier article and suggested that in the 
observations of the moon for the determination of Ephemeris 
Time, the vernal equinox of Newcomb’s fundamental catalog 
be utilized. 


= 


Space Range 


In a number of papers, the management, implacement or 
station keeping of space range system observatories have been 
considered. Among these were two papers by Koskela, 
Nicola and Walters [the second, (131), was a condensation of 
the first (130) ], in which the station keeping implications of an 
artificial satellite were studied. They attacked the inverse 
problem: Given a very accurately known satellite orbit, from 
observation accomplished at a given station, determine the 
position of the station. Range and angle measurements were 
utilized, and the novel analytical differential correction scheme 
developed by Herrick (101) was applied. Remarkabl; 
enough, a significant advantage of range over angle data was 
demonstrated. As noted previously, Murry (186) also ex- 
amined the problem of fixing station positions from observa- 
tions of a geodetic satellite. Kahn (116), in a specialization 
of the general problem, utilized differential correction formu- 
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las to determine corrections to the latitude and longitude of 
Mark II Minitrack stations. In an effort to estimate the 
relative effectiveness of optical and electronic data collected 
by a space range system, Hertz (103) compared the orbits of 
satellite 1959 Alpha One as obtained from electronic data 
(prime Minitrack) with the one obtained from optical al 


Nunn) observations. 


Optical 


By far the most numerous in 1960 were papers dealing with 
the acquisition and reduction of optical data. In a number of 
these papers the positions of reference stars—the stars that 
are utilized as bench marks in order to fix the position of a 
satellite on a photographic plate—were considered. Danjon 
(47) considered the degree of consistency between the right 
ascension and declination of star positions cataloged in the 
northern and southern hemispheres. Fricke (68) surveyed 
the absolute catalogs for southern hemisphere fundamental 
stars, whereas Nemiro (194) discussed whether or not the list 
of reference stars in the southern hemisphere satisfied the 
rigid demands of catalog accuracy. Oort (201) proposed that 
very accurate positions be sought for selected stars of particu- 
lar interest. Because there are fewer observatories in the 
southern than in the northern hemisphere, particular atten- 
tion must be given to improving the quality and quantity of 
reference stars in this hemisphere (as a number of the listed 
papers have indicated) in order to accurately reduce southern 
hemisphere satellite observations. In this regard Wood (273) 
discussed the astrometric program at the Sydney Australia 
Observatory (including a discussion of minor planet observa- 
tions). Kustaanheimo (140) considered another central 
problem in the determination of satellite coordinates from 
reference stars—the time derivatives of the components of 
proper motion. A star atlas good to at least 1 min of are and 
constructed from Baker-Nunn camera films was discussed by 
Wackernagel and Brandenberger (263). 

In regard to the space vehicle rather than the star image, 
Hynek (109) deduced that as long as satellites must be ob- 
served through an appreciable atmosphere, the use of short 
flashes for geodetic purposes would tend to increase rather 
than decrease the positional error. Yachoonsky and Gourko 
(275) discussed the variations of the albedo of the first arti- 
ficial Earth satellite resulting from the action of external fac- 
tors. Some 59 measures of the positions of comets and of the 
minor planet Icarus were reported by Jeffers and Gibson (114), 
and Neirinck (193), and Bajcér and Visnovcové (6) discussed 
observations of Russian space probes as did Gourko (83). 
Wasel (264) investigated the refraction correction to optical 
satellite data and the construction of refraction tables that 
utilized culmination techniques. An instrument for measuring 
relative positions of objects on curved films (Schmidt camera 
image fields) has been developed by Miczaika and Nunn (167). 

On the more general side, Herrick (102b) summarized re- 
quirements for precision orbit determination; Shepler (239) 
considered the visual tracking of space vehicles; Haywood 
(96) studied the application of optical techniques to inter- 
planetary navigation; Lillestrand (151) proposed his ‘‘aberra- 
scope,” which makes use of Bradley’s aberration of light to fix 
star positions; and Dubner (54) considered the technique of 
star occultations for providing long-range detection of ob- 
jects in space. 


Electronic 


Although undoubtedly other papers exist, we were only able 
to find three papers concerned with electronic observations. 
The first paper by Cohen and Kemper (43) found a position 
error growth, owing to the noise in Doppler observations, in a 
direction tangential to the orbit. In another paper, Frye (70) 
investigated the use of Doppler data as an aid to navigation, 
but felt that exorbitant power requirements preclude Doppler 
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data as an aid to navigation for deep space missions. A third 
paper by Guier (89b) and the previously noted paper by Guier 
and Weiffenbach (90), concerned with the Doppler tracking of 
project Transit satellites, was reported by Siry of NASA. In 
this paper is shown that orbits accurate to about 1 km have 
been obtained by a group at the Applied Research Labora- 
tory. (See also section on orbit determination.) The interest 
in radio observations of Jupiter led Krampe (135) in another 
paper to list an ephemeris for the reduction of radio observa- 
tions of that planet for the years 1960, 1961 and 1962 (this 
listing was omitted from the ‘American E phemeris and 


Attitude Dynamics 


As emphasized by Roberson, attitude dynamics forms a 
very important division of astrodynamics, and in 1960 a large 
volume of papers dealing with this subject has appeared. 
One can roughly divide these papers into two groups: Those 
dealing with passive tidal stability (e.g., akin to that of the 
moon) and those dealing with active stabilization systems or 
instrumentation for such systems. 


Passive 


A number of investigators considered natural tidal stability 
from the control system engineer’s point of view. 
tion ‘Electromagnetic Forces.’’) 
applied root-locus analysis to a nonspinning body and indi- 
cated conditions for stability and instability. 
Reiter (251) noted that attitude changes take place in a com- 


(See sec-_ 
Wolfe and Arrow (272). 


Thomson and 


pletely passive spinning system owing to energy dissipation. — 


In this case spinning about the principal axis of maximum 


inertia is unstable and results in an increasing angle of pre- — 


cession. 
whose whip antenna dissipated energy, was cited. Kondurar 
(128a) found a particular solution for the motion of two bodies 


when the satellite is a spheroid and the planet is a sphere. — 
Niday (195) attacked a more general problem—that of a- 


spherical planet and an arbitrary mass distribution for the 
satellite. His treatment is claimed to be simpler than Rober- 
son’s more general analysis. Although it is usually a very 
minor effect, Moran (178) discussed the coupling between 
orbital and librational motions. Klemperer and Baker ex- 


tended their earlier research: Klemperer (124) generalized to — 
large amplitude librations and Baker (8) generalized both to 


librations of a dumbbell on a slightly eccentric orbit and to 
the plane librations of a prolate ellipsoidal shell. The follow- 
ing seven papers deal generally with the motion of an aspheri- 
cal satellite in Earth’s gravitational field and will be noted 
without specific comment: Beletskiy (17), Doolin (53), Frick 
and Garber (67), Dubosin (55), King (120), Triplett and | 
Doolin (258), ond est Goodstein (266). 


Active 


In five papers sensing devices for active attitude stabiliza- 
tion systems were investigated. Arck and Merlen (4) con- 
sidered horizon scanners; Roberson (224) looked at the prob- 
lem of sensors in general and also (225) considered gyroscopic 
sensing of satellite yaw; von Pahlen (262) investigated an 
accurate determination of attitude by optical means; Ritchie, 
Hanes and Hainsworth (221) examined the display-control 
aspects of manual attitude control in space. In several papers 
complete active attitude control systems were analyzed. 
Grubin (89a) derived limit cycle solutions, and concluded that 
for a constant disturbance the limit cycle period is the quo- 
tient of the gas impulse by the disturbance torque, and is in- 
dependent of all other system parameters. Haeussermann 
(93) compared various actuator methods for the attitude con- 
trol of space vehicles. Three research groups associated atti- 
tude control with particular satellite systems: The demands 
of precision position and attitude controls for synchronous 


As an example, the cylindrical Explorer I satellite, _ 


~ 
are 
« 
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satellites were discussed by Barnes (12); Roman (227) set re- 
quirements for attitude control of the satellite telescope; and 
Karrenberg and Roberson (117) considered guidance and 
control of the 24-hr communications satellite. Various gases 
and gas mixtures for a gas jet attitude control system were 
compared by Mueller (180). Ormsby (203) on the other 
hand, looked at a free reaction sphere satellite attitude control 
system. A number of alternate stabilization systems includ- 
ing inertial, solar scanner and geomagnetic sensor were dem- 
onstrated by Ftenakis (71)—in his oral presentation a number 
of excellent film clips of the actual stabilization process taken 
from onboard space vehicles were presented. Raushenbax 
and Tokar (217) discussed the attitude stabilization of a 
vehicle in interplanetary space. Finally, Gillespie (80) pub- 
lished some general notes on the attitude control of Earth 


satellite vehicles. 


pe. 

As indicated by Vargo, investigations in ie area of orbit 
optimization under the assumption of impulsive expenditures 
of onboard energy have proceeded modestly in the direction 
of increased generality. The usual criterion of characteristic 
velocity with no additional constraint conditions still forms 
the basis of much work, but departures from this type of 
formulation are now recognized as having immediate and long- 
term value in a fundamental understanding of the problems 
and their application to astronautical engineering. Thus a 
mass ratio criterion and miss distance size conditions were 
adopted by Breakwell (28) to determine optimal sequences of 
corrective impulses. In another paper Vargo (258) extended 
the well-known-minimum impulse theorem to a more general 
force field. The scope of orbit optimization has been en- 
larged by Vargo (259) to include sets of orbits as the primary 
variable wherein orbital interactions are paramount in at- 
taining a defined objective. Leitmann’s results (146,147) 
emphasize the necessity of careful analysis of optimum orbit 
establishment if vehicular and launch realities are inserted in 
the mathematical description of the physical situation. This 
was originally pointed out by Vargo (257). 

The route from Hohmann’s theorem to the present status 
of characteristic velocity and no constraint optimization is a 
tortuous one, noted Vargo, and the published results of the 
past year make necessary some attempt to distinguish the 
avenues from the alleys. To this end Vargo categorized cer- 
tain papers (see Table 2) on the basis of what assumptions 
were made in their problem statements. = 

1 Initial orbit assumption: 

a, circular; b, elliptic; c, other. 

2 Initial-terminal orbit plane ialuticuahip: 

a, co-planar; b, non co-planar. 

3 Terminal orbit assumption: 

a, circular; b, elliptic; c, other. 

4 Number of impulses assumed: 

a, two; b, three; c, any number. 


Orbit Selection and Transfer 


Table 2 Summary of orbit selection and transfer papers 
Assumption— 
Reference 2 3 4 
- 181 a a b 
59 4 c) a (a,b,c) 
183 b a a b 
219 a a b b 
48 a b a c 
271 b b a b 
241 b a b b 
50a b a a b 
182,184 b a a b 
107 a b a b 
155b b b b a 


With respect to the more theoretical side of variational 
analysis of optimization, Miele pointed out that the main 
difficulty is associated with the numerical solution of the set of 
equations of motion and Euler equations. Since the varia- 
tional problems of flight mechanics are boundary value prob- 
lems, the use of lengthy trial and error procedures is an un- 
avoidable necessity even when digital computing equipment is 
used. It is for this reason that the opening of new avenues 
toward numerical solutions is important. In this connection, 
the recent novel approach by Kelley (118) on the gradient 
theory of optimum paths should be mentioned, even thoug] 
it is too early to judge the relative merit of this theory witl 
respect to the indirect methods of the calculus of variations 

The sophisticated variational techniques and the simple: 
methods of the ordinary theory of maxima and minima hav: 
been developed along parallel nonintersecting lines. Yet 
according to Miele, the following questions arise: What, it 
any, is the relationship between calculus of variations and 
ordinary theory of maxima and minima for flight mechanics 
applications? Are there cases where the optimization on a 
local basis and the optimization on an integral basis lead to 
identical results Miele did supply the answer to these ques- 
tions (169). It has been shown that for a wide class of prob- 
lems in the mechanics of flight, calculus of variations and or- 
dinary theory of maxima and minima lead to identical results, 
as long as the inertia terms are negligible in the equations of 
motion. 

The application of variational techniques to problems of hy- 
pervelocity flight (in particular, re-entry problems) has prog- 
ressed considerably. Worthy of mention are the papers by 
Levinsky (149) on the re-entry of satellite vehicles and those 
of Miele and his associates (168,170,171) on hvpervelocity 
gliders and ballistic missiles. 

A recent paper by Leitmann (148) relative to rocket trajec- 
tories in a potential force field is mentioned by Miele. It has 
been shown that, for the problem of optimizing the thrust 
direction and the thrust modulus, no variable thrust sub-arc 
may appear in the comparison of the extremal arc. The 
latter is only composed of coasting sub-ares and sub-arcs 
flown with maximum engine output. 

Because of the great interest in optimization techniques, 
several survey papers have been published in the last year. 
The reader desirous of learning more on the subject is referred 
to the soon-to-be-published work of Lawden and Leitmann 
and to Miele (172). 

Leitmann (147) has also investigated the problem of mini- 
mum transfer time for a power-limited rocket and found that 
a nonlinear acceleration for rectilinear transfer is optimum, 
but that a linear acceleration is very nearly optimum and de- 
creases the transit time only a few per cent. Schindler (232) 
considered this same problem. Fiul and Braham (66) dis- 
cussed injection onto a 24-hr orbit at a specific longitude and 
optimized vehicles and trajectories for this purposes. Graham 
(84) considered the optimum manner in which to utilize initial 
and terminal impulses with coast in between. on : 

/ 
Applications 

Although many of the subjects covered in the foregoing 
divisions of astrodynamics in a sense represent applied science, 
we have collected in this section those research papers that 
exhibit direct application to space vehicle systems. Such ap- 
plied papers fall into three broad categories: Applications to 
Earth satellite, lunar and interplanetary missions. 


Earth Satellite 


By far the greatest emphasis in satellite applications in 1960 
was in the field of orbital rendezvous or intercept. No less 
than 21 investigators authored or co-authored papers in this 
area, and it is expected that as more exotic systems of geocen- 
tric satellites are designed the emphasis on satellite rendezvous 
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will continue. Duke, Goldberg and Pfeffer (58) presented a 
rather exhaustive treatment of errors in two-dimensional in- 
tercept as well as an extension to three-dimensional out-of- 
plane-errors. They demonstrated that initial injection errors 
will build up secularly with time. Clohessy and Wiltshire 
(41) discussed a terminal guidance scheme for the rendezvous 
or ‘‘assembly”’ of a number of unmanned satellites in space 
about one central control satellite. Niemi (196), who in- 
vestigated a terminal guidance scheme for a satellite rendez- 
vous, stated that in order to minimize fuel consumption the 
thrust level should be at least four times the total accelera- 
tion. Nason (192) assumed that thrust dominates gravity 
during the intercept maneuver and developed a guidance 
equation. Petersen, Swanson and Hoover (208) proposed an 
astrovehicle rendezvous guidance concept; Roberson (223) 
considered path control for satellite rendezvous; Morgen- 
thaler (179) investigated midcourse guidance as related to 
satellite interception. In another paper (248), Swanson and 
Petersen discussed the friendly rendezvous of astrovehicles; 
consideration was given here to fuel consumption, timing and 
the effects of Earth’s oblateness. It was concluded that com- 
mensurability is required between the vehicle’s nodal period 
and the sidereal day. Brunk and Flaherty (33) examined 
methods and velocity requirements for the rendezvous of 
satellites in circumplanetary orbits, and Felleman and Sears 
(63) looked at a novel guidance technique for achieving an 
orbital rendezvous. Under the assumption of tangential 
transfer and co-planar circular orbits (see the section ‘““Obser- 
vation Theory”), Skalafuris and Schiller (240) looked both at 
minimum landing power and minimum midcourse corrective 
power requirements for rendezvous. 

Passing next to the area of satellite control and maneuver- 
ing, we find that Fitz Gibbon (65) proposed an infrared scan- 
ner system and circularized an elliptical orbit by forcing the 
velocity vector to be horizontal and the speed to be the ap- 
propriate circular speed. Errors of a few miles are to be ex- 
pected, and it is felt that the proposed method is neither 
optimum from the point of view of fuel expenditure nor com- 
petitive with standard differential correction techniques. 
The foregoing problem of circularization has also been con- 
sidered by Gundel and De Bra (91). 

In a paper related to special perturbations, Pohle (211) 
utilized a computation technique to carry out an error analysis 
of satellite orbits in the presence of drag. Drag was con- 
sidered as a function of the true anomaly, and the procedure 
resembled that of Baker and Makemson [(9), chap. 9 and 
(10)]. Vargo (259) investigated patterns for the observation 
of (or communication with) Earth’s surface by satellites af- 
fected secularly by perturbations of the aspherical Earth. He 
defined the most efficient manner in which to accomplish ob- 
servation or communications missions on the basis of the 
minimum number of satellites required, and developed the 
stability of the system with respect to perturbations. Berger 
and Ricupito (20) corrected an error in their earlier work that 
was called to their attention by Gersten. The correct rela- 
tionships were derived earlier by Herrick simply and directly 
from the C and § series found in the ‘American Ephemeris 
and Nautical Almanac” and can be found in (8). Loh (154) 
looked at flight along minor circles for boost glide or orbital 
vehicles; Detra, Riddell and Rose (51) considered low alti- 
tude satellite recovery by means of a controllable drag canopy 
after an impulsive de-orbit thrust. 


Lunar 


Orbit determination schemes not only apply to the problem 
of defining an orbit on the basis of observational error, but 
also to the problem of designing orbits, given certain specific 
conditions. In this latter regard Riddell (218) utilized a 
technique for specifying azimuths, elevation angles and times 
for lunar impact trajectories, i.e., for determining lunar 
orbits. Although interesting, it is not clear that Riddell’s 
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method is a great improvement over that of Gauss, and, fur- 
thermore, it seems to be defective in that it allows azimuth 
and elevation angle to be the nucleus of the problem, whereas 
they are actually simply by-products. Sedov (235) showed — 
the development of the trajectories utilized by the Russians 
for their lunar vehicles. Kooy and Berghuis (128b) utilized a 
Runge-Kutta numerical integration technique and the Cowell 
method to compute Earth-moon trajectories—launch sites on 
the Equator and at Cape Canaveral were considered. An in- 
teresting suggestion was made by Kaempen (115) who pro- 
posed the assembly of a large lunar vehicle from a number of 
small ones while in transit from Earth to the moon. 

Six papers in 1960 dealt with various aspects of lunar land- 
ing. Lanzano (142) applied the Jacobi intgeral to the ter- 
minal guidance of space probes, and Liu (152) utilized simple 
analytical means for isolating those areas on the moon that are 
accessible to normal impacts. Theoretical calculations of 
lunar trajectories in three dimensions were published by 
Mickelwait and Booton (166); both individuals have been 
responsible for the space mechanics aspects of the United 
States’ Able lunar program. An estimate of the value of man 
in soft landing guidance of a lunar vehicle was carried out by 
Peske and Swanlund (207). They concluded that for the same | 
fuel minimization the man required the same amount of feed- 
back information as did the automatic control system. Grube 
(88) investigated terminal guidance for both lunar and plane- 
tary probes; Stehling (247) considered start-stop,throttling 
and continuous thrust for soft and hard lunar impacts. 

A paper was presented by Hoelker, Brand, Jean and 
Schwaniger (105) concerning the parameter sensitivity and — 
guidance viewpoints for circumlunar flight and return. 
Two other papers also considered the moon-Earth return 
problem; Gunkel, Klemperer and Hunter (92) found that 
only certain preferred regions of Earth are accessible to moon- 
Earth return trajectories, and Miele (174) showed that for 
every Earth-moon trajectory there exist three image trajec- 
tories, two of which are moon-Earth return. The presenta- 
of Miele was followed by a lively exchange of views in which 
Herrick commented that the Miele theorem could also be 
proved on the basis of a Taylor series, and Baker supplied a 
further proof of the theorem starting from the dynamical re- 
flectivity and reversibility principals utilized in statistical 
mechanics. In reply, Miele commented that the number of 
images existing for the restricted two-body problem is eight, 
for the three-body problem four; furthermore, with respect 
to the work of Baker, Miele noted that the number of images 
in the rotating coordinate system is one half the correspond- 
ing number of images in the inertial coordinate system. 


Interplanetary 


The approach first utilized by Herrick in 1946, of approxi- 
mating a four-body interplanetary orbit by three two-body 
orbits, formed the basis of the analysis of Breakwell, Gillespie 
and Ross (27). Battin (15) considered the problems involved 
in “turn-around” interplanetary trajectories, i.e., trajec- 
tories in which the gravitational attraction of the target planet 
is utilized to place the vehicle on the return leg of its journey. 
A further investigation of the turn-around paths of Battin 
(15) led Smith (242) to conclude that they generated more 
problems than they solved in terms of midcourse correction, 
timing and reliability. Solar sailing, low thrust and other 
propulsion systems were described by Tanguay and Hock 
(250). Roberson (222) in a general exploration of interplane- 
tary (and other) space vehicle systems indicated that the re- 
quired studies should be expected to be oriented according 
to certain mission characteristics. In order to circumvent 
the costly, difficult and time-consuming trial and error process 
of computing definitive variant orbits, Rodriguez (226) pro- 
posed the use of approximate two-body orbit analysis. Ac- 
tually, although this proposal was not specifically established 
for steering programs, Kierstead (119) did just this, and re- 
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lated error coefficients of miss distance to midcourse guidance 
uncertainty analytically by means of simple two-body form- 
ulas. Anthony and Fosdick (2) simply indicated how yu [the 
sum of the masses, see chap. 1 of (9)] can incorporate Earth’s 
second gravitational harmonic term and how this quantity af- 
fects the injection onto interplanetary orbits. 

In this subject of interplanetary navigation, recognition 
should be given to an excellent paper by Breakwell (28) (see 
section on optimization theory). Here we find a definitive 
formulation of the problem of exactly how many corrective 
impulses should be supplied and where. Breakwell discussed 
the tradeoff between too early and too late corrections and 
formed three conclusions: 


1 If launch error is large compared to other errors, then 
the first correction should immediately follow launch. 

2 Three-dimensional considerations indicate that the 
first midcourse correction should occur three quarters of the 
way to the target planet. 

3 After this correction, the next correction should occur 
two thirds of the remaining distance to the planet. 


(These conclusions are based upon a number of approxima- 
tions contained in the article.) Lorell (156) developed guid- 
ance equations based upon the point of nearest approach to 
the planet. He defined the orbit by position and velocity at 
a “guidance point’? much in the same manner as Baker and 
Makemson did [see (9), chap. 12]. Clemance (40) in a short 
article indicated that man will not have enough time for in- 
terplanetary navigation decision making, that the present 
heliocentric and planetocentric constants are adequate, and 
that the “line-of-position’” method of fixing orbital position 
will be utilized. All of these comments are subject to con- 
troversy, and his last conclusion in particular seems to disre- 
gard the differential correction procedures that Clemence, 
himself, has so aptly applied over many decades to planetary 
orbits. Just such a differential correction was proposed in 
chapter 12 of (9), by Safren (228) and by Porter (212). 
Porter’s paper was actually a critique of an earlier paper by 
Cashmore. Specifically, Porter noted the great circle inde- 
terminacy found in Laplacian orbit determination that pre- 
cludes interplanetary orbit determination from observations 
of a single planet; he discarded the descriptive orbital ele- 
ments and suggested position and velocity; finally, he notes 
that Cashmore should have proposed two different numerical 
integration intervals—a short one during thrust and a long 
one during free fall. Cashmore and Gordon’s rebuttal to the 
criticisms of Porter (and also of Lawden and Noton) is a 
gracious retreat, e.g., with respect to Lawden’s and Noton’s 
criticisms they agree “‘. . . that a good deal is to be said for not 
forcing the vehicle back on the precalculated reference trajec- 
tory...” 

Specific consideration of the sensitivities of non-Hohmann 
interplanetary orbits to initial condition errors are discussed 
by De Bra (49). In a well-illustrated and informative discus- 
sion of vehicle rendezvous with a planet, Mickelwait, Tomp- 
kins and Park (165) showed that difficult compromises among 
the conflicting demands of payload, radio communication 
guidance and available power must be made. Pfeiffer (209) 
looked at rudimentary launch guidance methods for space 
missions; Safren (228) discussed the differential cgrrection 
method of interplanetary navigation (see the foregoing re- 
marks by Clemence in this regard), and Zilezer (276) con- 
sidered space navigation by correlation. On this same ques- 
tion of guidance, Spence (244) found that ICBM guidance 
systems are adequate for a Mars launch, and that if the 
probable error of Rabe is correct (and that if the error is not 
the difference between the Rabe and Spencer Jones values for 
the solar parallax), then the error in our knowledge of the a.u. 
will still allow for placing a space vehicle in the near vicinity 
of Mars. 

Turning next to the landing phase, we find an excellent sur- 
vey paper by Hermann (99) concerning the atmospheric entry 


problem at the planets. Friedlaender and Harry (69) con- 
sidered the requirements of trajectory corrective impulses 
during this landing phase of an interplanetary mission. A 
self-contained homing instrument for midcourse was proposed 
by Stearns (246) together with final adjustment guidance in 
order to insure passage through the atmospheric landing 
corridor. Nonweiler (197) presented a survey of a planetary 
atmospheric entry problems, since they are related to inter- 


planetary navigational errors. 
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rocket as in the aircraft field. 


N THIS paper the term diieaieeaini bodies’ shall refer to 
bodies moving and propelled under water, i.e., primarily 
to submarines and torpedoes. 

The military applications of such bodies as fighting vessels 
and weapons are sufficiently understood to make further ex- 
planations unnecessary. However, in recent years the use 
of submerged bodies as cargo carriers has added commercial 
aspects to the military applications. 

Although the military submarine and torpedo have been 
well recognized for more than half a century, today’s interest 
in submerged bodies has been greatly increased by the 
advent of nuclear power in this field. The independence of 
nuclear power from atmospheric air has opened new areas 
for the application and performance of submarines, and has 
also stimulated the development of other submerged bodies 
to which nuclear propulsion may not be applicable. 

In the general examination of a technical field comprising 
a wide range of dimensions and operating conditions of fairly 
similar bodies, considerations of similarity have been found to 
constitute the simplest and most effective tool for the initial 
approach. This paper is therefore based primarily on simi- 
larity conditions (1).2_ These considerations permit the rep- 


Presented at the ARS Semi-Annual Meeting, Los Angeles, 
Calif., May 9-12, 1960. 

1 Professor and Director, Garfield Thomas Water Tunnel, 
Ordnance Research Laboratory. Member A 
2 Numbers in parenthesés indicate References at end of paper. 
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“Fundamental Astronautical Constants and 


Submerged Bodies 


The field of submerged bodies, such as combat or cargo-carrying submarines and torpedoes, is 
surveyed from the viewpoint of similarity considerations. The lift to drag ratio of the body is repre- 
sented as a unique function of the Froude number, leading to corresponding relations for the ratio 
The importance of machinery weight suggests for underwater 
va am a propulsion similar developments from the propeller to a ducted fan engine, and to the ramjet or 
For the hydraulic equivalent of the ducted fan engine, the type of 


of machinery weight to tota] weight. 


propulsion pump is related to the operating conditions. 
- ye 2 ~~ ments of low machinery weight and those of good efficiency and good cavitation resistance, which 
+ ean be partly resolved by intake from the boundary layer of the propelled body. For very high 
nh ee “Froude numbers, travel with fully developed cavitation must be given serious consideration. 
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There is a conflict between the require- 


resentation of essential characteristics of submerged bodies 


as a function of a single variable, which is a suitably defined 
Froude number of the body, as originally suggested by Dr. 
Kenneth 8. M. Davidson. The primary characteristic so 
represented is the drag to weight ratio of the body which 
determines the essential characteristics of the power and pro- 
pulsion plant, the minimum structural requirements of the 
body, and the choice between buoyant bodies and bodies 
sustained by dynamic action like airplanes. 

This paper also touches upon the problem of drag, with 
particular attention to the surface to volume ratio and the 
effects of boundary layer intake on the efficiency of propul- 
sion. Finally, cavitation is discussed as a limit of speed, 
leading to the problems of travel with fully developed cavi- 
tation. 


Drag to Weight Ratio of Submerged Bodies 


The ratio of the drag or resistance of a moving body to its 
weight is a natural measure of its efficiency as a means of 
transport. In aeronautics this ratio usually appears in tlie 


form of the familiar lift to drag ratio of airplanes. With 
respect to bodies submerged in water, it is natural to con- 
sider neutrally buoyant bodies first. 
weight W is obviously 


For such a body the 


volume [1] 


ARS JourNAL 


one 


whi 
of t 


so t 


A 
Fri 
ing 
to t 
On 
desi; 
skin 
vary 
mar 
the 
with 

Fi 
| 
0.01: 
simil 
of 1.1 
(in 
body 
bodie 
ment 
curve 


ing re 
some 
at his 


Dec! 


: 
W 
el 
tl 
a 
Th 
+ 


where y is the weight per unit volume of the water. 
The drag of the body is generally expressed in the form 


D = Co(pV2/2)A [2] 
where Cp is the drag coefficient and A may be any suitably 
chosen area connected with the hydrodynamic drag. op is 


the mass per unit volume, and therefore p is equal to y/g, 
where g = 32.2 fps?. 


Hence the drag to lift ratio is of = 
=a 
W  ~? 2 volume [3] 


For simplicity of the following operations we define 
A =rdL volume = (md?/4)L [4a] 


where d and ZL are a representative diameter and a repre- 
sentative length of the body, respectively, which may be 
determined from the definitions [4a]. Evidently these 
definitions apply primarily to fairly slender bodies of revolu- 
tion. 
With these definitions the drag to weight ratio becomes 
D V? 
[5a] 
The expression V?/2gd is obviously the square of a Froude 
number F, referred to the diameter of the body. 
If, instead of definition [4a], A is defined as the “frontal 
area”’ 


Ay; = md?/4 [4b] 
one arrives, for the drag to weight ratio, at the expression 
D 

W = Coy QgL = Coy F,? [5b 


where Cp, is the drag coefficient referred to the frontal area 
of the body, which is obviously related to Cp by a 


wd L = Co, 
— 
Co, = 4 Co(L/d) [6] 


Although the Froude number referred to the body length 
F’, is closest to that used by naval architects, the correspond 
ing drag coefficient is not so convenient as it is proportional 
to the length to diameter ratio of the body which is variable. 
On the other hand, the drag coefficient Cp is for a well- 
designed body of revolution not much greater than the pure 
skin friction (or flat plate) coefficient, and as such does not 
vary over wide limits at the high Reynolds numbers of sub- 
marines or torpedoes. In this paper we shall therefore use 
the Froude number referred to the diameter (F2) together 
with the flat plate coefficient Cp and Equations [4a and 5a]. 

Fig. 1 shows a double logarithmic plot of D/W against 
F? according to Equation [5a], using for 4C’p a value of 
0.015. The points shown along this line represent the drag of 
similar bodies of revolution with nominal body diameters d 
of 1.6, 16 and 64 ft, respectively, traveling at the velocities V 
(in fps) noted on the points. It is evident that the smallest 
body represents a torpedo-like vehicle, whereas the two large 
bodies are vessels of about 1000 and 64,000 tons displace- 
ment, respectively. Also shown in Fig. 1 is an approximate 
curve for surface vessels. Since the submerged body has 
no wave drag, it is better than the surface craft except when 
the latter is surface planing at very high Froude numbers. 

Of particular interest is a band of D/W values correspond- 
ing roughly to the lift to drag ratio of subsonic airplanes or 
some fully submerged hydrofoil configurations. This means 
at higher drag to weight ratios the buoyant body is less effi- 
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ou 1 Drag to weight ratio of water-borne bodies 


cient than hydrofoil or airfoil supported craft. The F'.? value 
where the drag line intersects this band is about 3, i.e., above 
this value a buoyant body has more drag than a corresponding 
craft supported by dynamic lift. This is true for a body the 
size of a torpedo at practically all useful speeds, and for a sub- 
merged 1000-ton vessel at speeds over 33 knots, whereas a 
1000-ton vessel at 30 fps (18 knots) or a 64,000-ton vessel at 
practically any foreseeable speed will be more efficient than 
the airfoil or hydrofoil craft. It is also of interest to observe 
that for bodies in the size range of torpedoes the drag begins 
to exceed the weight of a buoyant body at speeds over 50 
knots. This variation in the efficiency of travel results from 
the increase in resistance (as well as of other hydrodynamic 
forces) with the square of the speed and from the dependence 
of these forces on the surface rather than the volume of the 
body. 

The natural thought of replacing the buoyant body at — 
Fa> V3 by a completely submerged hydrofoil craft could be > 
effective in water only if the density of the body could be a | 
high multiple of the density of water, because only then — 
would the resulting reduction in wetted friction surface be — 
sufficient to overcome the trend depicted in Fig. 1. It is 
therefore easy to see why in air the airfoil craft has practically 
eliminated the buoyant airship, whereas the hydrofoil craft — 
has only its wings under water and its payload carried above © 
the surface in a much lighter medium, 


Machinery Weight Problem of Submerged Bodies — 


The changes in drag to weight ratio of submerged bodies | 
described in the preceding section lead immediately to the 
most critical general problem of all submerged bodies (as 
well as of aircraft and other vehicles), namely that of ma- 
chinery size and weight at increased speeds. The foregoing 
presentation shows that it is actually not the speed, but the 
Froude number that determines the propulsion requirements. | 

The most common understanding is that the weight of a 
powerplant is proportional to its power. This is certainly 
true for the power-producing or rejecting elements, such as 
the boiler and condenser, and with some reservations it is 
also true for the nuclear reactor. 

If the entire weight of the propulsion plant were propor- 
tional to its power, the power to weight ratio of the entire 
vessel would be the proper variable to consider. It could be ~ 
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derived from Fig. 1 and Equation [5a] by multiplying the 
drag D and the weight W by the velocity of travel V. Thereby 


WwW 


which might be expressed in nautical units as horsepower per 
ton knot. 
Multiplying through by V leads to 


power 
i = constant- J 
weight of the vessel 


If the machinery weight (W,,) is assumed to be proportional 
to the power, it follows with consideration of Equation [5a] 
that 


7 73 tar: 

thus demonstrating the criticality of machinery weight with 
respect to speed. 

The weight of the force-transmitting or force-producing 
parts of the machinery, however, cannot be expected to be 
proportional to the power transmitted, as long as the ve- 
locity of these parts is not limited by other considerations, 
such as centrifugal stresses. In such parts it is the force and 
not the power that determines the limiting force-carrying cross 
sections. It is, of course, known that the weight of such 
machinery cannot be proportional to the forces alone, for 
example, to the thrust of the propulsor, since the thrust is 
known to increase with the square, and the weight with the 
cube, of the linear dimensions. 

The so to speak “missing”’ linear dimension in the relation- 
ship between force and machinery weight, specifically thrust 
and weight, is obtained from the fact that the weight of a 
force-carrying structure is for a constant force proportional 
to the linear dimensions of the structure. Dimensionally 
this means that the weight is proportional to a moment. 
The dimensional consistency of this relationship can be 
shown as 

W.» = constant F = constant o»-d? [9] 


where 


W,, = machinery weight considered 


Ym = weight per unit volume of the machinery structure 

Om = stressin the machinery material 

d = any suitable linear dimension of the machine con- 
sidered 

F = any representative force transmitted by the ma- 


chine, for example the thrust of a propulsor 


Equations [9] may immediately be rewritten in the form 


W. = constant 2" F-d = 


Tm 


Ym 
-Ons 
constant 1 [10] 


where M simply denotes the product of a force times a length, 
but may be a moment transmitted by the propulsion ma- 
chinery, for example, a shaft torque. This product F-d = M 
does indicate that for the same forces transmitted the ma- 
chinery weight increases linearly with the dimension of the 
system for the same force action; i.e., the weight goes up as 
the force action is distributed over a propelling mechanism of 
increased dimensions. It will be seen that this fact is to some 
extent in conflict with the requirements of high efficiency of 
the propelling machinery. 

Putting the force F equal to the drag D, the machinery 
weight relation expressed by Equation [10] leads to the fol- 
lowing ratio of machinery weight to the weight of the vessel 


We Y= 
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square of the speed and is independent of the absolute size o/ 
the vessel. 

Propulsion machinery may be grouped in various ways. 
One form of distinction is that between machinery where the 


propelling force is applied in the direction of the motion of 


the body, and machinery where the propelling force is 

applied at right angles to the motion of the propelled body 

Listing i in each group the propulsion devices in the order o 
decreasing size and weight for the same thrust, we obtain: 

1 Machinery with force action in the direction of th: 
body motion: The paddle wheel, the piston, the pulsejet 
the ramjet and the rocket. 

2 Machinery with the propelling force action at righ 
angles to the motion of the body: The fish, the propelle: 
the ducted fan engine and possibly the turbojet engine. 

It is seen that with the exception of the paddle wheel an: 
the fish, the first group represents propulsion devices withou 
primary rotating machinery, whereas the second group com- 
prises primarily rotating propulsors. 

With respect to machinery weight it will be recognize: 
immediately that intermittent propulsion devices, such a: 
the piston, the pulsejet and the fish, should be at a distinct 
disadvantage, since the maximum forces are available fo: 
propulsion only for a fraction of the total time, whereas th: 
machinery has to be designed to withstand these maximum 
forces irrespective of the duration of their application. 

In the first group the largest and heaviest propulsor, the 
paddle wheel, can probably be eliminated from considera- 
tion, since it does not have any known advantages for the 
propulsion of completely submerged bodies. 

In the second group the largest propulsion device, the 
“fish,’’? may deserve some consideration since it really does 
not add to the submerged body any external propulsive means 
because the entire body appears to act as propulsor. 

All propulsion devices of the second group may be said to 
comprise primarily a number of fluid deflecting surfaces or 
hydrofoils moving relative to the body at right angles to the 
direction of the body motion of travel. Figs. 2 show why 
this description is believed to apply also to the ‘fish,’ in 
comparison with the propeller represented here by the de- 
velopment of a coaxial cylindrical section through this type 
of propulsor. The important variable of all propulsors of the 
second group is the ratio of the forward velocity of travel V 
to the velocity U of the propulsor at right angles to the direc- 
tion of travel, which is called the “advance ratio” J 


n = number of revolutions per sec ae 
d = diameter of a rotating propulsor - 
The greater the advance ratio, the greater will be the ratio of 
the force Fy exerted in the direction of the motion of the 
propulsor to the useful thrust 7. Since the force Fy deter- 
mines the required strength and thereby the weight of the 
driving machinery, it is seen that an increase in J will increase 
the machinery weight. For a rotating propulsor this means 
that the lower the speed of rotation the greater the weight 
of the driving machinery, a fact well-known by marine and 
aeronautical engineers. It is also seen from Figs. 2a and 2b 
that the higher J is for a given velocity of travel, the lower 
will be the velocity of flow w relative to the propulsor. A low 
relative velocity is favorable for good resistance against cavi- 
tation, but one pays for this advantage in increased machinery 
weight. 

Finally it is seen that propulsion by motion of the pro- 
pelled body (like a fish) is likely to be limited to fairly smaii 
angles 8 and thus to large advance ratios. Thus fish-like 


demonstrating that the machinery weight ratio of the force- 
producing and transmitting machinery increases with the 
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Fig. 2a Propulsion by motion of the propelled body 
Fig. 2b Coaxial cylindrical section through a rotating, axial flow 
propulsor 


propulsion involves relatively large forces to be exerted by the 
propelling body and thereby relatively heavy actuating 
machinery. 


Efficiency and Weight of Underwater Propulsion 


All means of propulsion that are likely to be used for free- 
running submerged bodies derive their propulsive force from 
the momentum of a stream or jet of fluid. The greater the 
mass rate of flow in this stream, the lower will be the jet 
velocity required to produce a certain thrust. Since the 
momentum of a given mass increases with the first power of 
its velocity, whereas its kinetic energy changes with the 
second power of the velocity, it is clear that the greater the 
mass flow of the propelling jet the lower is the energy (or 
power) required to produce a certain thrust. This fact is 
expressed by the well-known formula for the “propulsive 
efficiency’’’ of “water or air-breathing” propulsors 


1 
where 
V = velocity of travel 
AV = excess of the jet velocity relative to the craft over 


the velocity of travel 


Since the thrust 7 = pQAV, where Q is the volume rate of 
flow in the jet, it follows that 7, is lowered as Q gets smaller, 
as stated previously. A large Q, on the other hand, requires 
a large and heavy propulsor, which confirms the foregoing 
statement that the requirements of good propulsive effi- 
ciency are in conflict with the requirements of low machinery 
weight. In subsonic air travel the requirements of low ma- 
chinery weight have prevailed, since the usual jet engine has 
propulsive efficiencies in the neighborhood of 50 per cent or 
less at subsonic velocities of travel. The criticality of ma- 
chinery weight in connection with bodies submerged in water 
suggests that underwater propulsion may develop in the 
same direction as aircraft propulsion. 

Considering the previously listed spectrum of propulsors, 
the fish-like propulsion would have the best propulsive effi- 
ciency since it works on the largest rate of flow, i.e., the flow 
surrounding the entire body. However, its machinery weight 
may well be prohibitive with present-day means, and prob- 
lems of its mechanical construction have never been attacked 
seriously. 

3 The “propulsive efficiency’’ does not include the hydrodyna- 


mic efficiency (n,) of the propulsor as a machine which is about 
0.9. Therefore 


drag times velocity of travel a 
= = 0.9+n, 


shaft power in ft lb per sec 
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The propeller of a submerged body can easily be made 
sufficiently large to keep its propulsive efficiency well over 
90 per cent. There is no reason, therefore, to seek for higher 
propulsive efficiencies than those of large, slow running 
propellers. Their weight, particularly that of their driving 
machinery, however, will eventually become prohibitive at 
increased speed of travel. 

It should be noted that fairly large, yet fast running pro- 
pellers, i.e., propellers with very low advance ratios and high 
circumferential velocities, are feasible for certain applications 
and will lead to reduced machinery weight. The available 
information on such propellers, however, is not sufficient to 
permit detailed consideration of this possibility in this paper. 

Since an increase in speed of rotation will, because of the 
resulting reduction in shaft torque, reduce the machinery 
weight, the question must be answered whether one could 
not simply reduce the propeller diameter for a given thrust, 
and thus trade propulsive efficiency for reduced machinery 
weight. This is a serious possibility, but limited by the 
maximum thrust per unit area of the propeller ‘disk.’”’ The 
upper limits of ‘disk loading”’ have not yet been established 
for open propellers. The use of two counterrotating pro- 
pellers behind each other is obviously a step in this direction 
and can be helpful in reducing machinery weight. 

The next step toward a smaller, high speed propulsor is 
obviously a device analogous to the turbojet engine of air- 
craft. Because of the nature of water, only a hydraulic 
counterpart to the “ducted fan engine’’ (2) has been con- 
sidered, which is logically called a “pumpjet.”’ The essential 
step is that here the rotating element of the propulsor is 
placed into a close fitting casing or shroud. In fact, a 
“shrouded propeller’ (using a rotating shroud attached to 
the propeller blades) should be the first step from the open 
propeller toward the pump in a fixed casing. The pumpjet 
certainly should remove all limitations of disk loading, since, 
if desired, a multistage arrangement could be used. Thus 
the change in velocity AV through the jet unit could be in- 
creased at will, and thereby the rate of flow and the diameter 
of the propulsor decreased. This will permit an increase in 
the speed of rotation and corresponding reduction in ma- 
chinery weight, limited only by the sacrifices in propulsive 
efficiency one is willing to make. 

The next step, the ramjet, depends for its operation on an 
increase in the volume of the flowing medium after its pressure 
has been increased by means of the “ram’’ pressure of the 
incoming stream. In water this is possible in principle by 
the introduction of gas or vapor or by vaporizing the entire 
stream. The latter gives the most powerful propulsor, but 
its propulsive efficiency would be quite low for foreseeable 
speeds of underwater travel, since the discharging jet would 
be a vapor jet having a velocity at least 10 to 20 times the 
velocity of travel. According to Equation [13] this will lead 
to propulsive efficiencies in the neighborhood of 10 per cent. 
Nevertheless, the compactness and simplicity of such a device 
might make it highly desirable for applications where effi- 
ciency is not of decisive importance. Of course, partial 
vaporization or addition of finely divided gas may lead to 
better efficiencies, but the flow problems involved have not 
yet been explored sufficiently to insure the feasibility of such 
a process. 

The last step in the spectrum of underwater propulsion is 
obviously that to a simple rocket. As in the aeronautical 
field, this is the simplest and lightest propulsion plant for a 
given thrust, but also the least efficient. At moderate 
depths of submergence the jet velocity of a rocket might 
conceivably be two to three times as great as the jet velocity 
of an underwater ramjet with gaseous exhaust. Thus the 
propulsive efficiency of the underwater rocket may be in the 
neighborhood of 5 per cent or less, which obviously limits 
its use to very short-range applications. 

There exists a possibility of increasing the efficiency of 
underwater jet propulsion by means of some form of “aug- 
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mentation.” Fig. 3 shows the ideal ratio uf thrust aug- 
mentation by means of devices working on the principle of a 
jetpump. The ideal ratios of augmentation given are calcu- 
lated on the basis of the conservation of momentum during 
a process of complete mixing between the primary jet and 
the induced ‘‘secondary”’ flow. These curves neglect the 
effects of skin friction within the augmenting nozzle, as well 
as well as those of incomplete mixing. It is known that 
actual augmentation ratios fall considerably short of the ideal 
values. Nevertheless, if only half or two thirds of the ideal 
augmentation were actually obtainable, the efficiency of a 
device, such as a ramjet with vaporized exhaust, could be 
increased to values between 20 and 30 per cent. Since the 
efficiency of the primary jet nozzle is in this case much higher 
than that obtained with the internal machinery driving a 
rotating propulsor, it is conceivable that propulsion without 
primary rotating machinery may become competitive because 
of its unparalleled simplicity, compactness and reliability. 


The cavitation limits of a submerged body are usually 
represented in terms of a cavitation number defined by 


h + 32 aw 


V*/29 


p = absolute static pressure in the water away from the 
disturbing effects of the body (but at the same 
submergence) 

P» = vapor pressure 

V = velocity of travel, fps 

h = submergence, ft 

32 = average value of the atmospheric pressure minus 
the vapor pressure expressed in feet of sea water 


Cavitation Limits and Design 
of Rotating Propulsors 


pV2/2 


[14] 


Assuming cavitation sets in when the vapor pressure is 
reached, o becomes equal to the coefficient of minimum pres- 
sure on the body, i.e., the greatest pressure reduction below 
ambient pressure divided by the velocity pressure of the speed 
of travel. 

It is certainly possible to design submerged bodies and 
their appendages so that o does not exceed unity, and such 
values have regularly been obtained with models. However, 
for large bodies the actual cavitation number is found to be 
considerably higher than the value obtained by model tests. 
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Fig. 3 Approximate values of the thrust augmentation by 


This may be a result partly of genuine scale effects on cavi- 
tation [see (3)], but most of all of the fact that the actual 
body has surface irregularities which can increase the loca! 
cavitation number appreciably. The relationship between 
the pressure coefficient C,, of the parent body without any 
irregularities, the cavitation number of the irregularity alone 
;. and the resulting cavitation number op of the combination 
is 

or = Cy, + (1+ C;,,)o; [15 


The cavitation number oa; of the surface irregularity alone is 
primarily a function of the height of this irregularity h com- 
pared with the local thickness of the boundary layer 6. An 
approximate curve for this relationship is shown in Fig. 4 
for isolated sharp edged surface irregularities in the forn 
of a ridge normal to the direction of the flow (4). Althoug! 
the values of this cavitation number vary as a function of th: 
Reynolds number of the irregularity, the form of the irregu- 
larity and perhaps other variables, Fig. 4 nevertheless 
does permit an estimate of this important effect of surface 
roughness on incipient cavitation. To use this curve under 
practical conditions requires a knowledge of the local thick- 
ness of the boundary layer. Although this thickness is a 
function of the shape of the body and must be determined in 
individual cases experimentally, it is for the present purpose 
sufficient to remember that the ratio 6/z varies with increas- 
ing Reynolds number from about 0.035 to 0.01, where z is the 
distance of the point considered from the leading edge of 
the body or of any one of its appendages that may be subject 
to local cavitation. For example, at about 10 ft from the 
nose of a body 6/z will be about 0. ae and therefore 6 equals 
0.24in. A surface irregularity of 7g in. = 0.25 6 would have 
a cavitation number of about 0.8. If the parent body had at 
this point a pressure coefficient C,, = 0.5, the resulting cavi- 
tation number would be 0.5 + 1.5 X 0.8 = 1.7, i.e., about 
three and one half times as high than without the surface 
irregularity. If, furthermore, a propeller blade would have 
1 in. from its leading edge a sharp surface irregularity 0.002 
in. high, h/5 would be about 345 and therefore o; approxi- 
mately 0.4. If the pressure coefficient of the perfect blade 
was at this point 0.3, then the resultant cavitation number 
of the actual blade would be 0.3 + 1.3 X 0.4 = 0.82, again 
nearly three times its original value. Since genuine scale 
effects on the cavitation process itself (3) may account with 
scale changes of the order of 10 for increases in the cavitation 
number by factors of 1.5 to 2, one should expect total in- 
creases in cavitation number from the model to the full- 
scale ship from about five to 10 or more. However, only the 
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Fig. 4 Approximate cavitation numbers of sharp edged surface 
irregularities (Holl) 
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Fig.5 Relation between the cavitation number co, and the blade 
pressure coefficient C;, for rotating propulsors 


cavitation numbers of a nearly perfect model 
dicted, and these are the numbers that will be 
the following. 

The pressure coefficients of airfoils suitable as propeller 
or pump blade profiles have been explored and reported ex- 
tensively (5). Low drag airfoils of 6 per cent thickness show 
values between C;, = 0.2 and 0.25, but are very sensitive to 
changes in the angle of attack. Thicker profiles are less 


can be pre- 
discussed in 


sensitive, but have higher pressure coefficients, say between 
C, = 0.3 and 0.35.4 

For rotating propulsors the blade pressure coefficients C; 
apply to the velocity w of the flow relative to the moving blade. 
Thus the cavitation number of a rotating propulsor referred 


to the velocity of travel V is according to Equation [14] 
w? 


= 


= C, 


b 


o, = C, (1 + U?/V?) 


Therefore the lower the advance ratio, or the larger the 
peripheral velocity U for a given V, the higher will be o, for 
given values of C;. This is shown in Fig. 5, where it is as- 
sumed that V is also the velocity at which the flow approaches 
the runner blades. It is seen that an increase in advance ratio 
improves the cavitation performance while it increases the 
machinery weight, thus demonstrating the previously men- 
tioned conflict between cavitation and weight requirements. 

The design of propellers has been repeatedly described in 


* It should be mentioned that much lower pressure coefficients, 
say Cy, = 0.1, have been used successfully, but not without ap- 
preciable cavitation at the blades. 
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the literature [see (6 and 7)], so that the foregoing remarks 
should suffice. The next smaller propulsor, i.e., the pump- 
jet, has not yet been described extensively. 

If one replaces a propeller by a pumpjet with the same rate 
of flow and thereby the same propulsive efficiency, the pump- 
jet would be larger and slower running, and thus heavier 
than the propeller, because for good cavitation performance 
one selects the axial velocity of the flow through the pump 
runner slightly lower than the corresponding velocity through 
the propeller. This would not serve the purpose of easing 
the machinery weight problem. To do so one must select 
for the pumpjet a lower rate of flow and thereby a higher 
velocity increase AV and lower propulsive efficiency than for 

the propeller. The resulting higher disk loading may be 
- expressed by the pressure coefficient or head coefficient of the 
pump ¥ = 29H/U?, where H is the “head’’ or work per 
_ pound of fluid. The head is obviously related to the velocity 
increase through the propulsor 


2 
9 | 9 


_ (2V + AV) At (17) 
29 


and the flow is assumed to enter the propulsor at the velocity 


oftravel V. It follows that a 
AV\ AV 
= (2 
(2 + ) “a [18] 


This relation between y, AV and V/U is shown in Fig. 6. 
Also shown is the so-called ‘‘specifie speed”’ (n, = n VQ/H?!4) 
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which determines the type or form of pump that will satisfy 
a certain combination of operating conditions. For axial 
flow pumps n, is determined by V/U and yW[n, = constant- 
(1/p)3/4-VV/U] for a given ratio of hub to tip diameter 
which was here assumed to be 3. The range of n, suitable 
for single-stage axial flow pumys is indicated in Fig. 6. 
Much lower values require multistaging or radial flow design; 
much higher values lead to poor internal efficiencies when 
using a pump, but are quite suitable for open propellers. 

The cavitation number o, determines the flow coefficient 
V/U of a propulsion pump in a similar manner as expressed 
by Equation [16] for propellers. Therefore the o, value one 
must meet, and the AV/V value, calculated from the pro- 
pulsive efficiency one is willing to accept, will determine a 
point on Fig. 6 and thereby the type of pump or propeller re- 
quired. For example, if very good cavitation performance 
(low o» values) is required, one must select a high flow coeffi- 
cient V/U, say over 0.5. This will rule out the use of a 
single-stage axial flow pump if one desires AV/V > 1.0, but 
a two-stage axial or a radial (mixed) flow pump would be 
suitable, which unfortunately has a fairly low speed of rota- 
tion. If, on the other hand, V/U = 0.2 to 0.3 is acceptable, 
one can use a single-stage axial flow pump up to AV/V = 2. 
This unit will have an increased speed of rotation which is 
favorable for low weight. 

With respect to all cavitation numbers quoted in this sec- 
tion, it should be remembered that they apply to the design 
shape of the propulsor and will be subject to increases owing 
to surface irregularities as indicated in the beginning of this 
section. 


Boundary Layer Intake 


In the preceding sections it was shown that a reduction in 
machinery weight obtainable by reducing the size and increas- 
ing a speed of a rotating propulsor is necessarily connected 
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with a corresponding loss in propulsive efficiency. It will 
now be shown that this decrease in efficiency can at least be 
partly overcome by using boundary layer intake to the 
propulsor. 

The advantage of boundary layer intake as discussed here 
is based on the fact that by such an arrangement the inlet 
velocity to the propulsor is reduced. Since for a given rate 
of flow through the propulsor the thrust is proportional to 
the difference between the discharge and inlet velocity (V2 — 
V, = AV), whereas the work needed to produce this velocity 
increase is proportional to V.2 — V,2 = (2V,; + AV) AY, it 
is clear that a reduction in inlet velocity V; will reduce the 
work required to produce a certain thrust and will therefore 
increase the propulsive efficiency. 

With boundary layer intake the inlet velocity to the pro- 
pulsor obviously cannot be uniform. The first problem to 
be considered is whether or not a reasonably efficient intake 
is possible under these conditions. This problem was in- 
vestigated several years ago under the sponsorship of the 
Office of Naval Research, and is reported in (8). The results 
of this investigation are that boundary layer intake is possible 
with a slight retardation of the boundary layer immediately 
in front of the intake, and the ram efficiencies obtained lie 
between 80 and 90 per cent. 

To arrive at reasonably simple considerations on boundary 
layer intake, it is necessary to define a mean intake velocity 
V, [see (2)].. This mean value should represent the momen- 
tum as well as the kinetic energy of the incoming flow, lead- 
ing actually to two different definitions, namely, for the 


momentum 
A A 
where 
y = coordinate normal to the boundary 


part of the boundary layer taken into the propulsor 
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On the other hand, the kinetic energy is represented by the 
definition 


ve =f ay [20] 


Fortunately the two different definitions for V, lead, for 
normal boundary layer velocity distributions, to practically 
the same values for V;._ The ratio of this mean velocity to the 
free stream velocity outside of the boundary layer V/Vo is 
plotted in Fig. 7 as a function of Ay/6, where 6 is the total 
thickness of the boundary layer. 

The propulsive efficiency with boundary layer intake is a 


AV-V; 
(V2 — Vy2)/2 + K(V 02/2) 


1 
Vi/Vo + (AV/2Vo)(K/2)(Vo/ AV) 


[21] 


where K is a coefficient expressing the friction losses resulting 
from the pump casing, and Vi = Jp. 

Fig. 7 shows 7, as a function of the thickness Ay of the 
layer taken into the propulsor, assuming that K = 0.05. 
For comparison the corresponding propulsive efficiency 7 is 
calculated for K = 0. It should not be surprising that this 
efficiency exceeds unity, since the propulsor uses the energy 
put into the boundary layer by skin friction with respect to 
an absolute frame of reference. “Efficiencies’’ over 100 per 
cent have actually been measured on bodies of revolution with 
the propulsor working in the wake of the propelled body. 

The relation between the propulsion stream Ay and the 
boundary layer thickness 6 is obtained by the assumption that 
75 per cent of the total drag is effective in front of the intake, 
i.e. 


= 0.75-pAyV, AV 


where @ is the momentum thickness. 

Fig. 7 also shows for comparison an efficiency 7,, calculated 
from the preceding relationship by putting V; = Vo, i.e., 
for intake from outside of the boundary layer. The corre- 
sponding efficiency 7,, was calculated for K = 0. The ratio 
of »,/nz — 1 is plotted in order to show the relative increase 
in efficiency owing to boundary layer intake. It is seen that 
this increase in efficiency is a maximum when Ajy/6 is approxi- 
mately 0.2. It is in this very range of boundary layer flow 
to the propulsor that boundary layer intake becomes signif- 
icant and increases the propulsive efficiency from values in 
the neighborhood of 0.6 without boundary layer intake to 
0.8 with boundary layer intake. The former value may well 
be unacceptable for many applications, whereas the second 
value will usually be considered as adequate. 

Finally, Fig. 7 presents AV/Vo as a function of the inlet 
thickness. This curve shows, for example, that one can 
increase AV/V» from less than 0.4 without boundary layer 
intake to 1.2 with boundary layer intake without losing 
propulsive efficiency, which is 0.8 at the two corresponding 
points on the 7, and the y, curves. Such an increase in 
AV /Vo means a reduction of the rate of flow through the pro- 
pulsor by a factor of 3, and a corresponding increase in its 
speed of rotation, or reduction in torque, by V3. Boundary 
layer intake therefore appears as a rather significant means 
for reducing the weight of the propulsor and its driving ma- 
chinery. 

The diagrammatic sketch of a submerged body with 
boundary layer intake shown in Fig. 8 is intended to illus- 
trate another advantage of this arrangement. If properly 
developed, boundary layer intake can be used to prevent 
separation behind a fairly blunt body. The ellipsoidal body 
shown there has a length to diameter ratio of 2, which would 
have only half of the wetted surface as a more conventional 
body with a length to diameter ratio of 8, both having the 
same displacement. Within the approximations used in this 
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paper this would, in accordance with Fig. 1, reduce the drag 
to weight ratio of this body by a factor of 2 for the same 
speed. 

Travel With Fully Developed Cavitation 

Up to this point cavitation has been treated as a limit of 
the fluid velocities to be used. The question arises whether 
travel or operation with substantial cavitation is feasible, thus 
eliminating this limit just as the old “acoustic barrier’’ has 
been eliminated by supersonic airplanes, compressors and 
turbines. 

There is no known reason why travel or the operation of 
propellers, pumps or turbines with extensive cavitation should 
not be feasible, although the flow problems involved may be 
considerable. With respect to turbomachinery, the general 
problem is outlined in (9). The theoretical treatment of two- 
dimensional flow through vane systems with fully developed 
cavitation is described in (10), but has recently received added 
attention. From (9) it is clear that the problems involved | 
should be easier for a propulsor than for a pump, since a — 
runner with fully developed cavitation will primarily acceler-_ 
ate the flow which should be quite sufficient for propulsion — 
purposes. 

With respect to the submerged body itself, travel with — 
fully developed cavitation has not yet received the attention 
it may deserve. Since at very high Froude numbers the drag 
of a fully wetted, buoyant body will exceed its weight, one 
should expect that under these conditions travel within a 
cavitation void offers a real possibility of more efficient travel 
than the fully wetted state. Yet the theoretical difficulties | 
of describing this form of flow are doubtlessly very great, as" 
the weight of the body “planing in the bubble’’ poses a three- 
dimensional flow problem with constant-pressure boundaries. 
So far only plane, two-dimensional problems of this type 
have been solved. Not even the axially symmetric problem 
has been conquered. In addition, there are unsolved prob- 
lems of cavitation and two-phase flow as such, so that a solu- 
tion of the overall problem of travel with fully developed | 
cavitation appears to be a truly difficult one. Nevertheless, 
since propulsion with thrust values exceeding the weight 
of the vehicle is regularly accomplished in the rocket field, 
the possibility of underwater travel at speeds many times 
greater than those presently used can no longer be ignored. 


Conclusions 


1 The drag to weight ratio of buoyant submerged bodies 
increases proportionally to V2/2gd, whereas for surface ships 
it increases faster because of the wave drag. Under similar 
conditions large bodies can travel faster than small bodies. 


Propulsor 


Fig. 8 Diagrammatic sketch of a submerged body with boundary 
layer intake to the propulsor 
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For V2/2gd > 3, the buoyant body is less efficient than a 
completely submerged hydrofoil. If V*/2gd > 100, the drag 
of a buoyant body will be larger than its weight. 

2 As V2/2g9d increases, the propulsion machinery weight 
becomes more and more critical, forcing a progressive change 
from the propeller to the pumpjet, to the ramjet and to the 
rocket. The weight of force-producing or transmitting ma- 
chinery is proportional to the propelling forces and the linear 
dimensions of the machine, suggesting the use of small pro- 
pulsors with high speeds of rotation. The efficiency of the 
propulsor and its resistance against cavitation, however, go 
down as the machinery size is reduced and its speed increased. 

3 The conflicting requirements of efficiency and weight 
of rotating propulsors can be reconciled to some extent by 
using boundary layer intake in connection with relatively 
small, high speed propulsors. Boundary layer intake also 
permits the use of much lower length to diameter ratios than 
are conventional, and thus a reduction in drag. For pro- 
pulsors with gaseous jets, thrust augmentation by momentum 
exchange between the primary jet and a secondary stream 
may be feasible. 

4 Travel and propulsion with fully developed cavitation 
deserve serious consideration at high values of V?/2gd. 


ome 4 wage jay ny 


T HAS — felt that underwater jet propulsion systems 
could develop and progress in much the same manner 
as has aerial jet propulsion. As a result, the Office of Naval 
Research and the Bureau of Naval Weapons have spon- 
sored research and development programs in this field, this 
sponsorship dating from about 1944. For the purposes of 
this paper, rockets are included in this classification. Both 
the intermittent and steady jet engines have been studied: 
These include the hydropulse, the hydroduct, or ramjet, 
the hydroturbojet and the rocket. 

The problems under water are largely due to high resistance 
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of Underwater Jet 
Systems 


v arious underwater jet Sn era hydroduct, or ramjet, the hydropulse, the hydroturbojet, 
_ and the rocket, both solid and liquid—are described and discussed for underwater applications. 

_ Because of the increasing importance of deep operation, the solid rocket is becoming more sig- 

_ nificant. The other systems offer specific advantages of greater range, etc., for shallower operation. 
_ The limit of speed imposed by the body cavitation phenomenon is discussed, as well as that imposed 
_ by power density limits. From the latter standpoint it is shown that rockets and the other jet 
systems offer much greater underwater speeds, for reduced ranges, than are attainable with rotating 
These systems are also free of the problem of propulsor cavitation at all depths. 
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to forward motion, high ambient pressure, cavitation and 
availability of suitable propellants. For moderate speeds 
(less than 90 knots), propellers turned by conventional under- 
water engines are superior but, just as in air, for high speeds 
we must have systems capable of higher power—the jet 
engines. 

A brief description of some of these engines is given, fol- 
lowed by a general discussion. The first system to be de- 
scribed, the hydropulse, is included largely as a matter of 
historical interest. 


Hydropulse 


The hydropulse (see Fig. 1) consists, in its most successfu! 
embodiment, of a straight pipe barre] which enlarges slight] 
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‘Notes on the Power-Speed-Weight Relation- 
ship of Ships,’’ Experimental Towing Tank, (Davidson Laboratory), Stevens 
Institute of Technology, N. J., Technical Memo no. 97, March 1951; 
see also Davidson, K. S. M., ‘‘Ships,’’ General Lecture IX, Internationa! 
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at the forward end and joins a reed-type inlet check valve. 
Forward of the valve is a diffuser. The axis is normally 
aligned with the direction of motion. The fuel, which in one 
case is a liquid alkali metal, is injected intermittently into the 
water-filled duct downstream of the valves. The injection of 
the fuel, which is finely atomized by a spray nozzle, is ac- 
companied by a quick reaction at high pressure and tempera- 
ture, with the generation of much hydrogen gas and steam. 
The gas expands, pushing out the barrelful of water rearward 
as a piston. When the water has been expelled, the gas may 
be overexpanded for the case of static operation. This 
causes new water to be drawn in the forward end through the 
valves, and the cycle repeats. At high forward speeds over- 
expansion is not necessary as ram pressure fills the barrel. 
Thus there is no limit on the amount of gas generated during 
each cycle. To illustrate, a barrel 4 ft 6 in. long and with a 
4-in. ID, refilled and fired at 15 cps at a forward speed of 64 
knots, produced about 1000 Ib of thrust. The nature of the 
pressure variation in the barrel at the forward end can be seen 
in Fig. 2. The tests were conducted in the ring channel 
facility of Aerojet-General. 

This system ran powerfully and efficiently, but it was found 
to have certain inherent drawbacks for high speeds. The 
intermittent intake through the diffuser affected the external 
flow, and above 60 knots at shallow depth external cavitation 
flashed cyclically over the body. Using multibarrels from 
the same diffuser helped some, but not enough. The time to 
refill the duct even at high speed occupied 75 per cent of the 
cycle (see Fig. 2). Therefore, to keep the average thrust high, 
the reaction pressures in the duct had to be high and rise 
steeply on fuel injection. This made it very difficult to keep 
the valves and their cage together, and massively heavy con- 
struction had to be resorted to. The pressure was not only 
cyclic, which forced the fatigue limit stress to be used, but it 
was suddenly applied, which cut the allowable stresses about 
in half again. Also, as a matter of practicality, the ducts oc- 
cupied too much space inside the body. 

Specific fuel consumptions operating with molten lithium 
were as low as 5 lb of fuel per thrust-hp-hr. Allowing for a 
propulsive efficiency of about 60 per cent, this is equivalent to 
about 3 lb per shaft-hp-hr on a rotating powerplant. This is 
severalfold better than is being achieved with underwater 
chemical powerplants at the present time, but it is not good 
enough to override the practical difficulties mentioned. 

These factors, together with the lack of a good water- 
reactive fuel, caused this project to be halted. 


Hydroduct 


The hydroduct is an underwater ramjet which uses a heat 
source to flash sea water to steam at ram pressure. The 
device developed by Aerojet-General uses a specially de- 
veloped nongassing propellant which has since found wide use 
in rocket igniters. This propellant provides a highly con- 
centrated source of heat energy, and in addition the steam 
formed by sea water in contact with its flame contains 
virtually no uncondensables. This feature was provided so 
that the thermal cycles using it eventually could be made 
closed and therefore depth-insensitive. Sea water enters the 
nose diffuser and passes around the propellant, which burns 
endwise, and falls into the flames. Turbulators are provided 
to make the mixture of steam and combustion products uni- 
form. The steam generation is at nearly the ram pressure 
above ambient from which the expansion is rearward through 
a nozzle. Like other ramjets, the device must be boosted to 
running speed. 

Free-running missiles using this propulsion system have 
been fired successfully many times on the ranges of the Naval 
Ordnance Test Station. Speeds in excess of 100 knots have 
been attained. As can be seen, then, the bydroduct is sen- 
sitive to depth, since for a given speed the pressure ratio de- 
creases as depth increases. 
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Fig. 2 Direct hydropulse pressure-time record 


|! 
> 
Attempts have been made to make this system depth- 
insensitive by putting a jet condenser on the exhaust, and 
although successful operation was obtained in the test pit, the 
free-running missile had starting difficulties and did not 
achieve successful operation. 


Hydroturbojet 


A variation of the hydroduct—the hydroturbojet—was 
also worked on. This unit used a turbine driven pump to 
force free sea water into the combustion chamber. A hydro- 
fuel (water-reactive chemical) was used, in this case molten 
lithium. By using a pump in addition to ram pressure, the 
device would operate at slower speeds and accelerate from 
zero speed. 


ai 


The rocket has a high thrust per unit frontal area, which 
is determined by the diameter of the combustion chamber or 
casing and nozzle. This makes it admirably suitable for 
high underwater speeds, i.e., above about 90 knots. The 
liquid rocket is not suitable for deep running because of limi- 
tations on the chamber pressure and hence the pressure ratio. 
However the solid rocket can operate at very high chamber 
pressures, and hence operates well at great depths. 

Although underwater rocket speeds are necessarily lower 
than those above water, the high power and thrust available 
can be very effective. Because of the high drag, the coast 
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Fig. 3 External expansion of hydroduct exhaust from critical- 
pressure nozzle 
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Fig.4 External flow confinement of exhaust jet expanding from 
critical pressure 
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Fig.5 Plot of limiting cavitation speeds vs. depth for bodies of a 
basic form at different fineness ratios py 


distance of underwater missiles is also limited, so that the 
rocket motors are used as sustainers, and a high-acceleration 
boost phase is necessary for trajectory stabilization. 

The exterior ballistics of these high speed missiles can be 
very complicated, but solutions to the problems have been 
found and low dispersions and reproducible trajectories have 
been obtained. 

The effect of depth introduces complexities in regard to the 
nozzle exit to throat area ratio. If the area ratio is proper for 
shallow running it will be too great for deep running, an 
considerable drag increase and loss of thrust will occur at 
depths. However it has been shown by tests in the Aerojet 
ring channel that an underexpanded nozzle will work as » 
correctly expanded one if the missile body is streamlined, an«| 


if it is not cavitating seriously so that the flow can follow th» 
body contours. The flow pattern for this case can be seen in: 
Figs. 3 and 4. The external flow can be seen to expand th: 


jet cross section, and in so doing the pressure distribution 01 
the propelled body is altered to recover the thrust associate: 
with the further expansion. There is a limit, however, to th- 
pressure at the nozzle exit which can be expanded this way; 
ie., the excess pressure at the nozzle exit above ambient 

‘annot be too large in relation to the ram pressure. In the 
tests so far this pressure excess has been about one half of the 
ram pressure. 

By properly employing this effect the performance of : 
rocket nozzle can be optimized over a wide range of depths 
It should be added that the thrust of a properly expanding 
rocket is the same under water as it would be in air at the 
same back pressure. 


Underwater Speed Limits Set by Cavitation 
and by Power Capabilities 


Because of the relatively high friction of water, the maxi- 
mum speeds of underwater bodies will be limited. First, it 
should be pointed out that the cavitation phenomenon puts 
its own limits on speed, depending on depth, fineness ratio 
and fairness of the body under consideration. Eventually 
cavity-running bodies may be developed which can run be- 
yond the cavitation limiting speeds, but variable depth re- 
quirements make this problem very difficult, and for the pres- 
ent at least, controllable cavity-running underwater bodies 
must be considered unfeasible. 

A curve of cavitation limiting speeds for a family of fair 
bodies in which the fineness ratio and depth are varied is 
shown in Fig. 5. It can be seen that, except for very shallow 
operation, body cavitation sets fairly high limits and pro- 
pulsor cavitation will usually set far lower limits. 

From a power capability standpoint the maximum speeds of 
steadily propelled underwater bodies can be considered in 
terms of the size of the system which is presumed to occupy 
the same fraction of the propelled body as the size is varied. 
For illustration, a family of similar bodies is assumed to vary 
in size and over a range of representative fineness ratios. For 
rotating powerplants, the law of variation of power with 
size can be given by 


P = [1] 


where P is the power, D the diameter of the body (hence a 
characteristic dimension of the machinery) and K is a propor- 
tionality factor. Powerplants which follow this law are, for 
example, multicylinder reciprocating engines where the size of 
the cylinder is held constant and the size of the powerplant is 
increased by adding more cylinders. Other cases are the 
addition of more of a basic power unit to increase size. 
Other powerplants vary in capacity as 


P = K,D? [2] 


An example of these is a reciprocating engine where the num- 
ber of cylinders is fixed and the cylinders are simply scaled up 
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in proportion. Turbine engines where the size of the parts is T T 
scaled also follow this law. 

When these two equations are applied to the propulsion of a 
family of similar bodies, the fineness ratio is found to have no 
effect above a certain minimum value and the power limiting 
speed is found to vary with size, as is shown in Fig. 6, where 
the rising curve reflects Equation [1] and the horizontal line 
represents Equation [2]. The curves intersect at 80 knots 
and a 21-in. diameter, as representing an attainable point for 
both equations. The curves are shown shaded to indicate 
that the values are approximate. The effect of cavitation on 
rotating propulsors which would impose limits on speed ex- 
cept when very deep is neglected in this instance. 

The motors of rocket or modified rocket type which includes 


SPELO-eNGTS 


imum 


Man 


the hydroduct offer higher speeds for underwater missiles = 

than motors which utilize rotating machinery. It is true, ae OUMETER=PCET 

however, that efficiencies and ranges will be considerably less. Fig. 6 Maximum speed of an underwater body powered by 
Consider a family of underwater bodies in which the rotating machinery based on the limiting powerplant capacity 


rocket motors are as large as practicable, i.e., limited by the 
ody diameter and assuming chamber pressures attainable by 
the solid rockets. Table 1 gives the limiting steady speeds for 
underwater rocket propelled missiles. It applies, down to 
very great depths, for solid rockets and only for shallower 
depths for liquid rockets and the hydroduct. A further Table 1 
characteristic of these powerplants is that they are free from 


which can be fitted inside the hull 


Limiting steady speeds for underwater rocket 
propelled missiles 


propulsor cavitation limits at all depths. V-knots 900 640 450 320 
Table 1 neglects the influence of the cavitation limiting Fineness ratio ° 10 20 40 
- 
ory: 
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component characteristics. 
a “dead space’”’ coefficient. 


determined for: 


Arbitrary Speed and Range 


INCE the weight of a vehicle is equal to the sum of the 
weights of its components, the following equation can 
be written 


Wr = We_+Wi+ W2+ W, [la] 
Substituting from below, the following equation is obtained 


[1b] 


Collecting terms and arranging in descending order of powers 


of L 
(4B) ps _ (Fe) 


2] 


where 


a characteristic length, and 


7 


and ¢ is a shell weight coefficient. Shell weight is propor- 
tional to vehicle volume, other things being equal. Equation 
[4] follows, since the weights of both propulsion components 


are proportional to drag, and hence the vehicle surface area 


or two-thirds power of vehicle volume (W7/y)?'?, other things 
being equal. (The influence of Reynolds number variation, 
etc., is compensated in the iteration process.) 

It can be shown that the solution to Equation [2], a 
standard cubic equation? in L, is 


Presented at the ARS Semi-Annual Meeting, Los Angeles, 
Calif., May 9-12, 1960. 
MHD 


1 Consultant. Presently Director of Engineering, 
Research, Inc., Newport Beach, Calif. Member 

2 See “Handbook of Chemistry and Physics”’ (Algebra, Cubic 
~' mena Chemical Rubber Pub. Co., p. 295, 39th ed., and 
others 
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Factors Influencing the Size and 
Weight of Underwater Vehicles 


and associated shell, and minimum propulsion system volume. 
speed ratios for minimum propulsion system weights are near 1.2. 
for two shaft driven vehicles, one with an electrical and one with a thermal propulsion system. 


R. C. BRUMFIELD! Davie 


U. S. Naval Ordnance Test Station 


Pasadena, Calif. 


ms Algebraic equations have been written to give direct expressions for the total weight and volume 
of an underwater vehicle of arbitrary speed and range in terms of performance requirements and 
The two equations are related through an average specific weight anc 
Equations for weight and volume influence coefficients of 15 desigi: 
parameters are derived. Optimum running speeds and ranges for stern-chase vehicles have bee: 
Minimum propulsion system weight, minimum total weight of propulsion system 


For several cases studied, running 
Sample calculations are giver 


where 


P(e) = 1/4[(V1 + + + + 
(V1 + 6,2 — 1)2/3]3 


The general function F(@) is plotted in Fig. 1. Values 
taken from the graph are sufficiently accurate for most 
design calculations. 

A direct expression for K, + Kz can now be written for 
the shaft powered vehicle with arbitrary? speed V and 
range’ R 


5507 550n 
or 
K; = + 9] 
where 
n = drag power + - shaft power, a propulsive coefficient 


[11] 
Therefore 
Then 


3 The speed i is arbitrary in contrast to the conditions discusse: 
in the section on propulsion system weight. 

‘The analysis covers families of vehicles each operating at 
its rated speed and rated power output for which the unit weights 
of powerplant and propellant consumption are known. 
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Equations [13, 7, 6 and 5] can be used (in that order) to 
solve for the required vehicle weight, provided all the design 
parameters are known, including an average specific weight 
y. This parameter is related to component specific weights 
and dead space (‘‘unused usable volume’’) by the volume 
equation (see the following sections). 


Equation for Vehicle Volume, With 
Arbitrary Speed and Range 
Since the volume of a vehicle is equal to the sum of the 
volumes which may be assigned to its components plus the 
unused volume (dead space), the following equation can 
be written 
We Wee, Wi, Ws 
YPL ve Y2 


7 


Ys 


or 

— = 0 [16 
1—a— (16) 


where L = (Wr/y)'/3, a characteristic length (Eq. [3]). 
K, and K; are defined by Equation [4]. 
The solution to Equation [16], a standard cubic in JL, is 


which is the same function as Equation [6] (see Fig. 1) and 


For shaft powered vehicles of arbitrary speed and range, 
K, and Ky, have the values given by Equations [9 and 12], 


respectively, and hence 

K, , Ke + R 

20 


[19] 


1100 


Equations [20, 19, 18 and 17] can be utilized to calculate 
the required vehicle volume (W7/y) provided the required 
design parameters and coefficients are known, including the 
dead space coefficient a. The weight and volume equations 
can be combined to determine the relationship between the 
average specific weight of the vehicle (7), and the dead space 


coefficient a. See the following section. 


Dead Space, Average Specific 
Weight and Buoyancy 


Equations [5 and 17] can be equated to give a direct 
relationship between y and a. However, the resulting 
equation must be solved by a process of iteration, and hence 
a graphical solution is recommended. The resulting graphs, 
as well as the solution, are instructive. The procedure is 
as follows: 

1 Plot Wr/y vs. y from Equation [5]. This is possible 
since, other things being equal, 0. = 4.(y). 

2 Utilizing the same ordinate scale, plot Wr/y vs. a, 
from Equation [17]. This is possible since, other things 
being equal, 6, = 96,(a). 

Corresponding values of y and a can then be read directly 
from the curves, and y vs. a plotted as shown in Fig. 2. 
Other things being equal, the maximum average specific 
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Wr Wet 
= F(0,) = K.F (0, 17 
Y — a — ud 
where 
= 1/4[(V1 + 88 + 1) + (VIF — 
[18] 
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Fig. 1 The mathematical function F (@) 
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Fig. 2. Graphical solution of vehicle weight and velume equa- 
tions for the relationship between average specific weight and 
the dead space coefficient 


weight (ymax) for the vehicle corresponds to full volume 
utilization or a = 0. For this condition, the (minimum) 
value of 6, can be calculated from Equation [19], and in 
turn, Wr/ymax, the minimum vehicle volume from Equation 
[17] (with a = 0). In most designs it will not be possible to 
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completely eliminate dead space but only to minimize it; Concept and Use of Weight and Volume 


hence, a is greater than zero. The plots suggested indicate Influence Coefiicients 
the vehicle weight and volume penalties associated with 
less than perfect packaging of the components. Weight Influence Coefficients 


The b he vehicle is given by , ; 
In accordance with the preceding argument, the weight 


be Vw Wr of a vehicle can be described as a function of several design 
B= 1) cs (Yu — [21] parameters and coefficients. Stated mathematically 
which in level flight must be supported by up or down Wr = Wr(du dx, 92. -90) ~ , 


“Tift.” If a high angle of attack is required to compensate where n arguments are involved. ad 
for the buoyancy (or forces induced by maneuvering), the The total differential of Wr is 


drag coefficient will be increased. Since this is usually a » ow o * 
second-order effect, it can be compensated in the iteration dWr= 
process. i=1 
2 For small variations 69;, we have, approximately 
n Wr 


We can introduce weight influence coefficients C;. for 
each parameter g; such that ~ 
‘ 


qi OWr 
iv = 2 
Wr 04: 
As seen from Equation [25], the weight influence coefficient Th 


of any parameter indicates the degree to which small changes 
in that parameter influence the total vehicle weight. 

The solution for the weight equation can be written, from 
Equation [5] 


r(e) 


or 
= — 0) [28 | whe 
The weight influence coefficients are, from Equation [26] 
+ 0? ~ (Wi + + where fun 
- coe! 
The general function f(@) is plotted in Fig. 3. (See Ap- oe 
pendix A.) ol 
The weight influence coefficients can easily be calculated U 
from the foregoing equations. Sample calculations and wits 
tables of weight influence coefficients appear in Appendix B. br 
Volume Influence Coefficients = 
a 
An analogous argument holds for the volume influence 
By : coefficient of any parameter. It indicates the degree to 
—_ which small changes in that parameter influence the total! Ste: 
Se . vehicle volume. The following equations can therefore be : 
presented Ir 
= 
\ opt. — = K,F(6,) = dv [31 | 
where K, is given by Equation [17]. Thus 
whe! 
Se 
chasi 
inter 
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which is identical in form to Equation [30]. The general 
function f(6) is plotted in Fig. 3. These volume influence 


coefficients can easily be calculated. See Table 4. 


Propulsion System Weight and 


Up to this point, the analysis has been concerned with 
vehicles required to cover an arbitrary running range at an 
arbitrary running speed. In some cases the running speed 
and range can be related analytically. A simple case is 


that involving stern chase of another vehicle. im oi 


=!» 


Stern Chase 7 


In the case of simple stern chase from an initial range Ro, 
the vehicle running range required for overtake is given by 


Ro Ror 


y= [38] 


r 


where 
= targetspeed = 
r = speed ratio 


Several optimization criteria can be introduced for stern- 
on the 


nd 


chase vehicles. The criterion chosen will depe 
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Propulsion System With Minimum Weight 

For a shaft powered vehicle, the weight of the propulsion — 
system is given by 

Wes = Wi + We = (Wr/y)?*(Ki + Ke) [39] 

For stern chase, we have, from Equations [38, 39 and 13] 


It is convenient to introduce the term 


Hence = 
Cl + Op (Wr?! 
_ + Dp (Wr 42 


For the constant volume vehicle, Wr/y is fixed, and the 
minimum propulsion system weight corresponds to the 
condition 


(OW wo 0 [43] 


which leads to 


Tw = (Bw — 1 + + 1)/3aw [44] 


thus defining an optimum speed ratio for the vehicle. Equa- 
tion [44] is plotted in Fig. 4 in general form. 
It is thus apparent that this optimum running speed ratio 
is a function of a, only. The running range required at 


optimum speed is given by hs 


R/Ry = 2+ + 1 [45] 


as shown in Fig. 5. 
The weight of this optimum propulsion system at minimum 
is given by 


Pin, « x 


11007 
(V 3a. + 1 + 2)! 
3a,(V 3a, + 1+ 1)? 


[46] 
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Fig. 5 Running range ratio for optimum speed ratios 
Therefore 
| 
4 
4 
1155 
| 


indicates the penalty in propulsion system weight paid for 
nonoptimum running speed. The ratio reduces to 


(Wes)o 

as shown in Fig. 6. 


Propulsion System With Minimum Total Weight 
Including Associated Shell 


Since the shell weight is proportional to volume, we can 
assign a portion of the shell weight in proportion to the volume 
of the propulsion system. We therefore require that the sum 


be at minimum with respect to the running speed ratio, and 
that the total vehicle volume be held constant. Following 
the argument of the preceding section, it can be demonstrated 
that this corresponds to 
Bae’ — 1+ +1 


= [49] 


based on constant total volume for the torpedo, where 
1 + 


and a, is taken from Equation [41]. 


T 


b 


Fig. 6 Weights of propulsion systems for off-optimum speed 
ratios 
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The ratio of weights given by Equations [42 and 46] | 


Propulsion System With Minimum Volume 


Given a set of design criteria and parameters, there is an 
optimum running speed for a stern chase which will give a 
minimum volume propulsion system. We hold the total 
torpedo volume constant and require that the sum 


Wi Ws 
v1 v2 


[51] 


be at a minimum with respect to vehicle running speed v or 


speed ratio r. Following the forementioned argument, it 
can be demonstrated that 


ry = (8a, — 1 + V3a, + 1)/8a, [52] 
based on constant total volume for the vehicle, where 
Vo (h/¥1) tr 
a, = — ay = 3600 — (53 


The form of r, is again similar to that for r, and ry’ (Eqs. 
[44 and 49]) and Fig. 4 can be used. The argument amounts 
to putting h and T on a volume basis. 

In this section, several criteria for optimum running speeds 
are developed. For purposes of differentiating the equations 
involved, these are based on constant drag area (W7/y)*”, 
and hence apply to propulsion system variations alone. 
The values y and a therefore vary to hold drag area constant. 
The required variations in y and a can be calculated from 
the solutions to the weight and volume equations (Eqs. [5 
and 17]). 

Optimum running speeds based on minimum overall 
vehicle weight or volume can be calculated from the equations 


resulting from setting C,. and C,, to zero. 


Application of the Method 


The following examples serve to illustrate the use of the 
method, the type and form of the results, and the manner 
in which variations in design parameters influence weight 
and volume. Two families of small underwater vehicles 
are studied for weight characteristics. These are shaft 
powered and driven by thermal and electric prime movers, 
respectively. The following base values of parameters and 
coefficients are used for illustration: 


Speed of pursued vehicle, »% = 50 fps. 

Stern-chase initial range, Ry = 3000 ft. 

Sea water density, p = 2 slugs/ft*. 

Drag coefticient based on (volume)*”’, C, = 0.0237. _ 

Tail drag increment, t = 0.2. 

Propulsive coefficient for thrust producer = 0.85. 

Shell weight coefficient, ¢ = 16.1 lb/ft’. 

Average specific weight, y = 64.4 lb/ft® (neutrally 
buoyant). 

Payload weight, Wp, = 100 lb. 

Thermal engine weight/hp = 1.13 lb/hp. 

Thermal propulsion system specific weight of range 
dependent components = 18 Ib/hp hr. 

Electric motor weight/hp = 2.5 Ib/hp. 

Electric propulsion system specific weight of range 

7 dependent components = 40 lb/hp hr. 


The parameters shown in Table 1 were calculated for the 
two vehicles assuming design for minimum weight of pro- 
pulsion system. 

The corresponding weight influence coefficients were 
calculated and tabulated in Table 2. These are listed in 
descending order of importance (absolute value). Negative 
values indicate decrease in weight with increase in value of 
the parameter. The totals indicate the sensitivity of the 
weight to design improvements on a relative basis. Thus, 
favorable changes of 1 per cent in each of the parameters listed 
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would decrease the vehicle weights by 6.2 and 12.5 per cent, 
respectively. 

In addition to these calculations, certain of the design 
parameters were allowed to vary independently over wide 
ranges. The results are shown graphically in Figs. 7 to 10. 
The influence coefficients are the slopes of the weight curves 


(on a dimensionless basis). en 


Appendix A: Mathematical Functions yo 
F(6) and f(6) 


Function 


The solutions to both the weight and volume equations 
involve the function F(@) which is plotted in Fig. 1. 
For 0 < @ < 0.3 we can use the approximation 


F(0) = 1 + 1.890 + 2.856? [A-1] 
as shown in Fig. 11. 
For 0.3 < @ + 1.0 we can use the approximation 
= 0.919¢74/8 [A-2a 


which plots as a straight line on semi-log coordinates as shown 
in Fig. 11. The approximations have no special utility, 
except in visualizing the nature of the weight and volume 
relationships. 


For >> 1 (very large values of 6) 
F(0) = (27/4)68 


Accurately calculated values of F(6) are included in 
Table 3. 


[A-2b 


Function f(§) 


The combination of terms appears 


Table 1 Calculated parameters 


Thermal Electrical 

Parameter vehicle vehicle Reference 
aw 3.76 3.76 Eq. [41] 
Tw 1.22 1.22 Eq. [44] 
V, fps 61 61 ae 
R/Ro 5.50 5.50 Eq. [45] 
K, 15.6 34.5 Eq. [9] 
Ky 18.7 41.5 Eq. [12] 
Ow 0.294 0.652 Eq. [7] 
F(@) 1.82 4.20 Fig. 1 
2.22 2.38 Fig. 
6f(6) 0.65 1.56 Fig. 12 
Wr 242 558 Eq. [5] 


Values 


Parameters Thermal Electrical 
Vo (target speed) 1.59 3.83 
y (av. specific weight) —0.92 —1.72 
W pr (payload weight) 0.78 0.52 
0.65 1.56 
—0.65 —1.56 
r 0.36 0.85 
Ro 0.36 0.85 
h 0.29 0.71 
t 0.11 0.26 
Total (absolute values) 6.19 12.54 
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peatedly. This has been defined as f(@) and plotted in 
Fig. 3. Thus 

—— A-3 
and 
= 


3{(0/V1 + 63) [(V1 + + + (V1 + — + 1} 
(V1 + + + (V1 + — + 8] 


[A-4| 

Table 3 Numerical values of F'(6), and 

| 

0.00 1.00000 1.88989 0.00000 
0.05 1.09447 1.95297 0.09765 
0.10 1.21523 2.00808 0.20081 
0.15 1.34557 2.06553 0.30983 
0.20 1.49396 2.12061 0.42412 
0.25 1.66326 2.17266 0.54316 
0.30 1.85614 2.22070 0.66621 
0.35 2.07662 2.26420 0.79247 
0.40 2.32774 2.30101 0.92040 
0.45 2.61337 2.33181 1.04931 
0.50 2.93774 2.35557 1.17778 
0.55 3.30817 2.37121 1.30417 
0.60 3.72450 2.37976 1.42786 
0.65 4.19548 2.38051 1.54733 
0.70 4.72488 2.37304 1.66176 
“0.75 5.31889 2.36044 1.77033 
0.80 5.98245 2.34067 1.87254 
0.85 6.72128 2.31530 1.96800 
0.90 7. 54036 2.28518 2.05666 
0.95 8.44657 2.25097 2.13842 
1.00 9.44192 2.21366 2.21366 
1.05 10. 538319 2.17351 2.28219 
1.10 11687928 2.13344 2.34678 
1.15 13.04079 2.08840 2.40166 
1.20 14. 461078 2.04444 2.45333 
1.25 16.0012 1.99999 2.49999 
1.30 1766324 1.95561 2.54229 

| 
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Fig. 11 The function F(@) showing exponential and series 
approximations 
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@ Accurate values of @f(@) are included in Table 3. 


Table a Influence coefficients* for arbitary speed and range” 
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Ro/vo(r — 1) for stern chase. 


Note that the exponential approximation of Equation [A-2] 

implies that f(@) = %. This is only a rough average value 

for the center portion of the curve from 0.3 to 1.0 (Fig. 3). 
For 6 >> 1 (very large values of 6) 


= 3/0 [A-5] 


= 3 A-6] 

Accurately calculated values of f(6) or @f(6) 
in Table 3. Curves of f(@) and 6f(@) are shown on Figs. 3 
and 12, respectively. 


wey 
Appendix B: Formulas for Weight and 
Volume Influence Coefficients oe. 


Weight Influence Coefficients for Vehicles of r 
Arbitrary Speed and Range 
These are calculated from Equation [29], using aidilie 
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arbitrary speed and range, with Ro/v(r — 1) substituted for 

cient of the speed ratio, r. 
Equation [44]. 


[28] for K,. Ten parameters and coefficients are of interest: 
Co, t, n, h, T, R and v. The results are given in 
Table 4, first column. 


Volume Influence Coefficients for Vehicles of 
Arbitrary Speed and Range 


These are calculated from Equation [36], using Equation 
{17] for K,. Fifteen parameters and coefficients are of 
interest; the 10 listed plus: a, y:, yrz, Yi and y2. See 
Table 4, second column. 


Weight Influence Coefficients for 
Stern-Chase Vehicles 


These are, in general, the same as those for vehicles of 


In addition, we can develop a weight influence coeffi- 
If this is set to zero we obtain 
See Table 4, third column. 


1159 


— 
ip 
= 
i= 
3 
4 
| 
: 4 
« 
( 
| 
i 
i 
a 


Volume Influence Coefficients for function of @ appearing in the influence coefficient 


funet 
Stern-Chase Vehicles (Fig. 3) 
= f(0) ford = 
These are, in general, the same as those for vehicles of S(O.) = f(0) for 0 = Ow 


function of @ appearing in the solution to the cubic 


arbitrary speed and range, with R/»(r — 1) substituted for F(@) 
equation “ag 5 and 17] 


R/v. In addition we can develop a volume influence coeffi- 


cient for the speed ratio r. If this is set to zero, we obtain 
Tobie 4 h = W, + shaft output (unit weight of power- 
plant), lb/hp 
Nomenclature K,, K, = weight coefficients, lb/ft? oo 
K. = vehicle volume parameter, ft* a 
a = dead space (unused usable volume) coefficient K, = vehicle weight <a Ib 
B = vehicle buoyancy (net), lb L = characteristic length for the vehicle, ft 
Cie = vehicle volume influence coefficient for parameter q; 4: = typical independent variable 
Ci = vehicle weight influence coefficient for parameter q; r = vehicle running speed ratio = v/v 
C, = body drag coefficient based on (volume) 1 rf = speed ratio for minimum volume of propulsion system 
D = vehicle drag force, lb fe = speed ratio for minimum weight of propulsion system 
- : i = speed ratio for minimum weight of propulsion system 


ee 


R = required running range of vehicle, ft 

Ro = initial range at start of stern chase, ft 

t = tail drag increment 

v = vehicle running speed, fps 

Up = target running speed, fps 

V = unused usable volume (dead space) in vehicle, ft* 
WwW, = total weight of range independent powerplant com- 


ponents (prime mover, power transmission, bearings 
and accessories), lb 

W, = total weight of range dependent propulsion com- 
ponents (propellant and tankage, etc.), lb 


Wp, = weight of payload, lb 
Wrs = W, + Ws, total weight of propulsion system, lb 
W, = total weight of shell, including tail, lb 
Wr = weight of entire vehicle, lb 
ay = modified form of aw 
Qw = prime mover/energy source weight ratio parameter 
Oy! = modified form of ay 
¥ = average specific weight of entire vehicle, lb/ft* 
"1 = average specific weight of Wi, lb/ft? 
v2 = average specific weight of Ws, lb/ft* 
pL = average specific weight of W pz, lb/ft* 
Ys = average specific weight of W,, lb/ft* 
Yw = specific weight of ambient medium (sea water), lb/ft? 
ly = W, per shaft horsepower hour (unit weight of energy 
source), lb/hp hr 
” = propulsion coefficient 
Ow = @ for solution to vehicle weight equation ‘_ 
Oy = 6 for solution to vehicle volume equation i 
p = density of ambient medium (sea water), slug/ft? 


special symbol for vehicle volume in solution to vehicle 
volume equation [=K,F(6,)], ft? 


Fig. 12 The product 6f(6) 
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Selection of High Performing 
Propellants for Torpedoes | creer 


Texaco Experiment Incorporated 


selection of chemical propellants for heat engines of underwater missiles is discussed on the 
basis of the maximum performance available. Here it is evident that reactions in which *‘free”’ sea 
7 - water is both oxidant aad diluent provide the best performance, and aluminum and zirconium may 
- oh be the best fuels. Possible methods are mentioned for bringing about the necessary chemical re- 
J actions. 


HE NUMBER of candidate propellants for use in pro- 
pulsion of underwater missiles or torpedoes by means of a 

heat engine is possibly as great as is the number considered for 
use in rocketry. Also, in general, similar considerations apply 
in choosing the optimum system for a given application in 
either space or underwater applications. These include the 
important factor of performance capability, based on the 
intrinsic thermodynamics of the propellant, and practical 
considerations of suitability under use conditions. This 
paper is concerned primarily with providing an understanding 
only of the intrinsic performance capabilities of the propel- 
lants. 

In contrasting the processes involved in this country in the 
actual selection of propellants for underwater and rocket 
usage, it is notable that in underwater applications greater 
emphasis has been placed on practical considerations and that 
the attainment of a truly high performance capability 
generally has been relegated to a secondary role. In this 
regard, there has been a far greater tendency to apply “blue 
sky’’ techniques to achieve spectacular success in rocketry 
than has been true of underwater studies. In fact, quite 
often underwater propellant development waits on successful 
accomplishments in rocketry and, where possible, adapts these 
successes to torpedo use. 

Set forth are the factors involved in selecting chemical 
propellants for use in heat engines of underwater missiles 
based only on the maximum performance available, per- 
formance figures on several candidate fuels, and means of 
carrying out the propellant reactions. 


Discussion 


The heat engines used in torpedoes have a combustion 
chamber in which the propellants are reacted to liberate their 
energy. This chamber closely resembles a rocket combustion 
chamber and may use liquid or solid propellants. Liquids 
generally are pumped mechanically to the chamber and may 
be oxidizer, fuel, monopropellant or water, the latter taken 
from an internal tank or pumped from the environment. The 
hot gases produced may be used in a mechanical engine, 
either turbine or reciprocating, which subsequently powers a 
propeller or a pump; or they may expand through a nozzle 
to provide thrust as does a rocket. 

Considering now the carried propellant alone (i.e., excluding 
water taken from the environment), there is available in a 
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given torpedo a certain space where propellant can be stored. 
Thus a logical propellant performance parameter is a numeri- 
cal representation of the effect on performance of filling this 
space with the various propellants considered. This com- 
pares performance on an energy per unit-volume basis and 
neglects weight effects. Fortunately, propellant weight 
effects generally are not definitive in torpedo applications for 
reasons next discussed. 

In the first place, the proportion of the torpedo weight that 
is the carried fuel is relatively low, probably never exceeding 
50 per cent and usually less. (This is in decided contrast to 
rockets where ratios exceeding 90 per cent may be realized.) 
Thus the total weight of the missile is not critically dependent 
on the propellant weight. Also, weight effects in underwater 
applications are diminished by the bouyancy produced by the 
water. 

Another factor contributing to lower the importance of 
weight is that most torpedo “flights” are of the constant speed 
or sustainer type with speed changes occurring only over a 
small part of the movement. Thus, through most of the 
flight the torpedo powerplant is concerned only with counter- 
acting drag. At high speeds drag is much more a function of 
external configuration than of weight, or in other words 
“drag due to drag”’ is far greater than is “drag due to lift.’’? 

A further consequence of the sustainer type of flight is that 
work, or more succinctly, power output becomes the definitive 
criterion. This is because the powerplant is operating against 
a given force, drag, in sustaining the system, which is tanta- 
mount to the classical concept of the production of work by 
displacing an object against a force. 

Thus, in light of the previous discussions where it is shown 
that volume and not weight effects are important, and that 
work or power output is involved in sustainer applications, the 
propellant parameter of most usefulness in characterizing 
underwater propellant performance is power output per unit 
volume or, in customary units, hp-hr per ft*. This parameter 
is of definite value in grading new propellants for underwater 
applications, although it should be understood that any final 
choice between propellants requires a rigorous study at the 
exact conditions of the contemplated flight plan. 

The hp-hr per ft* parameter may be calculated from the 
energy available from the chemical reaction and the type of 
working gases produced. This calculation is similar to rocket 
specific impulse calculations and involves the determination 
from thermodynamic considerations of the enthalpy drop on 
isentropic expansion of the exhaust products from a 
chamber to ambient conditions. is 


2 Or ‘‘drag due to dive’’ in certain 
i 161 
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Fig. 1 Approximate theoretical power available on weight and 
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Fig. 2 Approximate theoretical power available on weight and 
volume bases for some water-reacted propellants vs. moles of 
free water (expansion ratio, 1000/14.7 psia; open data points 
indicate combustion temperature about 1800 F) 
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Fig. 3 Approximate theoretical specific impulse and density 
impulse of some water-reacted propellants vs. moles of free 
water (expansion ratio, 300/14.7 psia; open data points indicate 


combustion temperature about 1800 F) 
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Fig. 4 Approximate theoretical specific impulse and density 
impulse of some water-reacted propellants vs. moles of free 
water (expansion ratio, 1000/14.7 psia; open data points indicate 


combustion temperature about 1800 F) 
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Performance parameters for a number of propellants are 


gathered in Table 1 and in Figs. 1 through 4. The data were 
calculated by an approximation technique and are probably 
correct to within 10 per cent. In Table 1 they are arranged 
in descending order of power output per unit volume for 
reasons previously discussed. In Figs. 1 and 2 performances 
in terms of power are given for various amounts of “free 
water” used (water removed from the environment); in 
Figs. 3 and 4 the data are given in terms of rocket parameters 
for comparison. 

It is evident from the table and figures that high perform- 
ance with respect to both weight and volume, or specific im- 
pulse and horsepower, generally does not characterize any 
single propellant, emphasizing the need for a selection of a 
correct performance parameter for use in making the com- 
parisons that have been discussed. It is also evident from 
Table 1 that highest performance on a power per volume basis 
is attained with fuels that can react energetically with free 
water, and lesser performances result when the oxidizer is 
carried. Also, certain metals, such as aluminum, zirconium, 
magnesium and lithium provide particularly energetic per- 
formances when used with free water or even with a carried 
oxidizer. 

With certain metals established as the most energetic per- 
formers, the problem remains to achieve reaction at the neces- 
sary high efficiency within the torpedo combustion chamber. 
This often requires unconventional approaches, some of which 
are mentioned. 

One feasible method of attaining reaction under the neces- 
sary controlled conditions is to provide means in the torpedo 
for premelting the metal, which then is pumped as a liquid.* 
The necessary heating for melting the metals may be 
accomplished electrically before firing the torpedo or by more 
novel means based on heat transfer from the combustion 
chamber to the metal.4 Metals which might be considered 
here include aluminum, magnesium or lithium. An appro- 
priately injected spray of these metals may be reacted with a 
liquid oxidizer or with water drawn from the environment. 
The latter becomes feasible even with ordinarily passive 
metals because of an increase in reactivity when the metals are 
at or above the melting temperature. 

Another means for bringing about controlled reaction is to 
incorporate the metals as finely divided particles in a slurry, 
which is then pumped to the combustion chamber. The 
slurry might be toothpaste-like to prevent settling, or a stirring 
mechanism could be incorporated in the torpedo. Generally, 
it is envisioned that such slurries would be reacted with a 
potent liquid oxidizer, such as nitric acid or hydrogen per- 
oxide, and would not prove sufficiently labile for reaction 


3 Work along this line was performed at the U. 8. Naval Ord- 
nance Test Station, China Lake, Calif., about 1940. 

* Suggested by Texaco Experiment Incorporated, Richmond, 
Va., 1955. 


Table 1 Approximate theoretical performances of some 
underwater propellants. (Formulations diluted with 
water to decrease combustion temperatures to about 
1800 F; results in horsepower for 20.4/1 expansion ratio) 


Carried hp- hp- 
No. Fuel oxidizer hr/ft?  hr/Ib 
1 114. 0.70 
2 106 0.27 
3 LiClo, 78 0.57 
4 71 
5 LiClO, 60 0.55 
6 BrF; 59 0.29 
7 HO (90%) 55 0.55 
8 H.02 (90%) 50 
9 LiClO, 45 0.7 
10 ae. 41 1.27 
11 H.2O2 (90%) 40 0.68 
12 H.0, (90%) 28 «0.36 


(90%) 28 0.35 
HO. (90%) 26 0.34 


13 diesel fuel 
14 alcohol (92.5%) 


15 hydrazine mono- 23 =0.28 
propellant? 
16 typical solid pro- ... 23 «(0.31 
pellant 
17 Na LiCl,O 21 0.31 
18 12 0.14 
20 Na eae 5 0.24 
21 alcohol (92.5%)° compressed 


air 
“ Free water used as oxidizer as well as diluent. 
» Liquid-monopropellant formulation, 8% N»H,, 
N2H;NO;, 20% 
© Diluent water carried. 
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directly with free water. However, it is possible that ex- 
ceedingly fine metal particles (in an appropriate dispersion) 
might be spontaneously reactive with water.4 

As another means of promoting reactions, the metal may be 
incorporated as a powder along with an oxidizer in a solid 
propellant grain. The metal-oxidizer proportions within the 
solid propellant may be stoichiometric, in which case reaction 
of the metal would take place with the solid oxidizer; or an 
excess of metal could be used. In the latter case, the metal 
would be preheated to a labile temperature as a result of par- 
tial reaction within the solid charge and so would react 
spontaneously with water. These systems resemble hybrid 
rocket propellants in which liquid and solid propellants are 
used together. 

Needless to say, these methods of utilizing metals in the 
combustion chambers of torpedoes are not simple, common or 
conventional. However, for reasons of exceptionally high 
performance they may be worth developmental study. 


° 
t 


s 
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ERODYNAMIC ablation is the process by which the 

thermal energy generated in the layer adjacent to a solid 
body in high speed flight is absorbed by the solid, which 
melts and/or evaporates and is transported downstream and 
away from the solid surface. It is well known that Earth is 
protected from the bombardment of cosmic fragments by 
aerodynamic ablation. Solid bodies falling into Earth’s 
atmosphere from outer space attain very high speeds, be- 
come meteors because of aerodynamic heating, and usually 
dissipate before reaching Earth’s surface. The same phe- 
nomenon, aerodynamic ablation, must be engineered effi- 
ciently to protect long-range ballistic missiles re-entering 
the atmosphere, or vehicles returning to Earth from outer 
space. 

Theoretical and experimental investigations have shown 
that the injection of relatively cool gas into the boundary 
layer is an efficient way of alleviating aerodynamic heating. 
If there were no mixing or heat transfer between the injected 
gas and the outside air, the injected gas would form a heat- 
insulating, low temperature layer. However, heat is trans- 
ferred from the hot air to the injected gas just as it would be 
transferred to the solid surface, and also the hot air and cool 
gas are mixed by molecular diffusion and/or turbulent eddy 
motion. Consequently an idealized insulating layer does not 
exist, but the temperature gradient near the surface is con- 
siderably reduced, and the rate of heat transfer to the wall 
is less than it would be without injection. It is clear from 
these considerations that the effectiveness of injection cooling 
is increased with an injected gas of large heat capacity, low 
heat conductivity and low diffusivity. 

In the ordinary mass injection scheme, the rate of injection 
is an independent parameter. However, in aerodynamic 
ablation involving evaporation of surface material the rate of 
evaporation, and hence the rate of injection of gaseous mate- 
rial into the boundary layer, is determined by the rate of heat 
transfer, which in turn is reduced by the injection of gas. If 
one borrows the terminology of electrical engineering, the 
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‘ome Model 2 at 


_ tunnels is examined on the basis of available laminar boundary theories. 
such experiments may yield basic results useful for the understanding of aerodynamic ablation. 
With such materials as water ice, dry ice, naphthalene and camphor, a fairly wide range of import- 
ant parameters of ablation phenomenon may be covered. An experiment was performed with a 
- } water ice, hemisphere-cylinder model in the GALCIT M = 5.8 wind tunnel at a supply temperature 

of 300 F, and the results show: That the ice melted around the stagnation point when the suppl) 
__- pressure was above 34 psia, and the ice sublimed everywhere on the model when the pressure wa: 
less than 34 psia; and that the rates of ablation heat transfer were 0.69 + 0.11 cal per cm? sec at the 


TOSHI KUBOTA? 


Guggenheim Aeronautical 
Laboratory, 

California Institute of Technolog) 

Pasadena, Calif. 


The possibility of an experimental study of aerodynamic ablation in moderate temperature wind 


The results indicate that 


These experi- 


a ss stagnation point, and from 0.15 to 0.03 cal per em?’ sec on the cylindrical afterbody. 
mental results compare favorably with theoretical predictions. 


input function is the total enthalpy difference between the 
external flow and the gas adjacent to the wall; the output is 
the heat conduction into the solid, and the evaporation forms 
the negative feedback loop in aerodynamic ablation. Hence 
the output increases less rapidly with increasing input than 
it would without this feedback effect. Thus Lees (1)? finds 
that the effective heat of sublimation increases with increasing 
rate of ablation, or increasing enthalpy potential. 

Extensive experimental investigations of ablation cooling 
utilizing high temperature facilities such as an arc-heated jet 
are being performed in search of suitable materials. Un- 
doubtedly high temperature facilities provide much valuable 
information. On the other hand, thermodynamic and fluid 
mechanical aspects of the phenomenon depend on certain 
similarity parameters, and the question arises as to whether 
a sufficiently wide range of these parameters can be simulated 
in moderate temperature, steady flow wind tunnels under 
carefully controlled conditions by choosing suitable materials. 

This paper examines the possibility of such an experimental 
study in the GALCIT hypersonic wind tunnels, and discusses 
some results obtained with water ice models. 


Aerodynamic Ablation—Sublimation 


It is assumed that the boundary layer is laminar, and that 
the gas in the boundary layer is a binary mixture of air and 
the vapor of the material subliming under the action of aero- 
dynamic heating. The analysis given in (1) is simplified 
for the case of a binary mixture without chemical reaction. 
From the conservation of mass of the subliming material 


— (pDi2) w (OK /dy) + Pw wK w =m {1] 


where K is the mass fraction of the vapor of the subliming 
material. From the conservation of mass of the air molecules 


(pD2) (OK /Oy) w Pww(1 = 0 [2] 


since the concentration of the air molecules is 1 — K and the 
net flux of air molecules at the wall is zero. From these two 


equations, one obtains 


= = (pDj2) /Oy) w 
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Fig. la Mass addition parameter and surface temperature, 
water ice 


From the conservation of energy 


K 
Pwlwhw = dw — 
w 


oT 
E oy — oy 
[4] 


where hy, he and hy are static enthalpies of the vapor, air and 
the mixture, respectively. The first term on the left-hand 
side of Equation [4] represents the “conduction” of energy 
from the gas to the wall, and the second term represents the 
convection away from the wall. On the right-hand side, 
the first term is the conduction of heat away from the surface 
in the solid wall, and the second term represents the con- 
\ of energy toward the surface in the solid. 

Le, = 1 


k oh 


and Equation [4] becomes 


then 


hw) m(hw CoT'w) = b,w [5 


Furthermore, if Pr = Le = 1, then the concentration K and 
the total enthalpy h, satisfy identical equations, and there- 
fore 


K= Kulhs, = hs)/(Rse hw) 
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Fig. lc Mass addition parameter and surface temperature, 
camphor 


Then from Equation [3] 


and from Equation [5] 
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If Kw, q@.. and Lg are known as functions of T.. then the 
surface temperature 7’, may be determined from Equations 


[8 and 9]. For that purpose, a mass addition parameter 
B’ is defined such that 
B’ = m/puCu [10] 
Then Equations [8 and 9] become fit 
= [11] 
= (he — hew)/(Le + Lr) [12] 


where Lr = @.u/m 

If one assumes that the vapor pressure of the subliming 
material is equal to the equilibrium value corresponding to the 
surface temperature, Equation [11] becomes 


and the surface temperature is given as the slaicas at which B’ 
given by both Equations [12 and 13] are equal. The corre- 
sponding value of B’ is determined at the same time. Several 
such plots are shown in Figs. 1. 

A more refined method of determining the surface tem- 
perature on a subliming flat plate, without the assumption 
of Pr = Le = 1, is given in (2). The present approximation, 
however, applies to any geometry. 

Once B’ and T,, are determined, the ablation rate m is deter- 
mined as follows: From Equation [10] 


~ M21 — (p./p) 


pelle Ho Cu Cx [14] 


Ww B= 
the injection rate of foreign gas into the boundary layer, and 
Cxu/C xo is the ratio of the rates of heat transfer with and with- 


mh/peueC y, is the commonly used parameter for 


out mass injection. The ratio Cy/Cx, is ordinarily given as 
a function of the parameter B. Symbolically 


Cu/Cu, = f(B) [15] 


The variation of Cy/Cxz, with B’ may be obtained by sub- 
stituting Equation [14] into Equation [15] and solving for 
Cu/Cx, as a function of B’. When the relationship between 
Cx/Cx, and B is given graphically, the solution is obtained 
by drawing straight lines Cx/Cz, = B/B’ for fixed values 
of B’ and reading the values at points where the curves of 
Cx/Cx, vs. B and the lines Cx/Cz, = B/B’ cross each other. 
For the laminar stagnation point flow the relation may be 
written with good approximation as 


The dependence on the molecular weight was found by cor- 
relating the theoretical results for a flat plate boundary layer 
in (2), and is assumed valid for a stagnation point flow. 
Then m is determined from Equations [12 and 14] as 


= Go/(Le + Lr)(1 + BB’) 


where 

qo = peeC x (he 
B = 0.64(M,/M2)—'/ 


Equation [16] may be writtenas 


Le + Lr + — hew) [18] 


This relation shows that the effective heat of sublimation is a 
linear function of the enthalpy potential as experimentally 
verified in (3). 

This assumption, Equation [13], regarding the concentra- 
tion near the surface, needs a little more critical examination, 
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since in general the vapor pressure is not at the equilibrium 
value when vapor is being removed from the surface. The 
difference between the vapor pressures under the equi- 
librium condition and under the condition of finite mass trans- 
fer is as given in (4 and 5) . 


(Prog Pv) w = (tn/a)V 24ReT 


where a is the coefficient of evaporation representing th 
fraction of the impinging molecules that do condense. Unde: 
the present test condition of po = 15 psig, T> = 300 F, an 
with the water ice model of 0.5-in. nose radius, a calculatio: 
based on this assumption yields an ablation rate of 0.12 > 
10~* gm per cm? sec at the stagnation point. The correspond 
ing difference between the equilibrium vapor pressure and th: 
nonequilibrium vapor pressure is 0.08 mm Hg for a = 1 
which is about 2 per cent of the equilibrium vapor pressure 
For a = 0.1, the difference amounts to 20 per cent of th: 
equilibrium vapor pressure, and the correct wall tempera 
ture will be higher than the value calculated by the presen 
method. However, since the slope of the curve of vapo 
pressure vs. temperature is rather large, the error in th 
temperature amounts to only a few degrees, provided a> 0.1 


Equipment, Test Procedure and Data Reduction 


The experiment was performed in the GALCIT hypersoni: 
wind tunnel leg 1, which is of the closed-circuit and con- 
tinuous running type. The tunnel has a contoured solic 
nozzle which gives a nominal test section Mach number oi 
5.8. The maximum supply air condition in leg 1 is at 100 
psi and 300 F. During most of the present test, the tunne! 
was run at supply pressures of 20 and 15 psig and a supply 
temperature of 300 F. 

As shown by Fig. la, the rate of mass addition with water 
ice is small for the supply conditions of the M = 5.8 tunnel. 
Therefore, techniques can be developed for this case, wher 
the fluid mechanical aspect of the phenomenon is fully under- 
stood, and then applied to conditions with large mass transfer. 

The water ice models were hemisphere-cylinders of 1-in. 
diameter and approximately 5-in. length. In order to avoid 
machining ice, and also to facilitate the installation of thermo- 
couples, the ice was frozen onto a micarta tube in a brass 
mold. Several thermocouples were installed on the micarta 
tube for measuring the ice temperature (Fig. 2). A con- 
siderable amount of effort was spent to obtain clear, bubble- 
free ice in the freezing process. First, ordinary tap water 
was frozen in a mixture of dry ice and acetone, which yielded 
white waxlike ice. Freezing the water at a temperature 
slightly below the freezing temperature resulted in clear ice, 
but it contained a large number of bubbles in the core of the 
ice. The process finally adopted was as follows: A jar full 
of distilled water was de-aerated by boiling under vacuum; 
the de-aerated water was introduced into the mold, and the 
mold was put into a leaktight vessel; the vessel was then 
evacuated to eliminate the gas absorbed by water during the 
transfer from the jar to the mold, and also to minimize aera- 
tion during freezing; finally, the vessel was placed in a freezer 
which maintained a temperature of about 20 F. 


Test Procedure 


The tunnel was started about half an hour before the test 
and warmed up, in order to minimize the transient period of 
the tunnel supply temperature. After the warmup period, 
the air was short-circuited through the bypass line, the mode! 
was installed in the test section, and then the supply air was 
again switched through the tunnel for the test. 

It took about 1 min for the supply temperature to reach the 
steady-state value. Hence, at exactly 1 min after the supply 
air was switched from the bypass line to the tunnel, the pho- 
tographic and temperature recordings of the ablation were 
started. The 35-mm photographs of the back-lighted mod: | 
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were taken at either 1- or 3-min intervals depending on the 
rates of ablation. The temperatures at several locations in 
the model were read on a Brown self-balancing potentiometer. 


Data Reduction 


Using an optical comparator, the rates of ablation were 
determined from the 35-mm photographs. The _ posjtions 
of the stagnation point were determined with respect to a 
fixed point from a series of photographs, and the ablation 
rate was determined by plotting the nose position against the 
time. Except for an initial period of 3 or 4 min, the ablation 
rate was almost constant with time (Fig. 3). Body diameters 
on the cylindrical afterbody were measured at several loca- 
tions fixed with respect to the initial nose position. Hence 
the ablation rates on the cylindrical part were the average of 
those at the top and the bottom of the model. 

The errors involved in the measurements of linear dimen- 
sions were at most +0.0005 cm on the 35-mm negatives, which 
corresponds to +0.003 cm in the actual dimensions. The 
estimated probable errors in the ablation rates were +6 X 
10-4 em per min which amounts to relative errors of +0.9 
per cent at the stagnation point and +14 per cent along the 
cylindrical afterbody. The temperature measurements were 
accurate within +1 F. 


Experimental Results and Discussion 


Melting Around the Stagnation Region > 


As may be seen from Fig. la for a given recovery tempera- 
ture, there exists a value of pressure at which the wall tem- 
perature is just at the melting point of the ice. For pressures 
above this critical pressure, the rates of heat transfer from the 
air to the wall are sufficiently high to bring the surface above 
the melting temperature. When the pressure at the wall is 
less than the critical value, the ice does not melt but sublimes. 
This conclusion was confirmed in the test, as shown in a series 
of photographs (Fig. 4). At M = 5.8 and 7) = 300 F, this 
critical condition corresponds to a reservoir pressure of 20 
psig. When the tunnel supply pressure is much higher than 
20 psig, the ice around the stagnation point melts, and forms 
a craterlike depression. Since the water has low viscosity it 
is easily convected away from the stagnation region and re- 
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5.0" (approx.) -| 
Thermocouple Positions 


from nose 


@ 0.7 in. 

1.2 in. 

1.7 in. 
@) 2.2 in. 

6 3.2 in. 

4.2 in. 


Fig. 2 Model construction 


freezes in a low pressure region, and thus forms the ridge — 


around the stagnation region. As the pressure is lowered, | 


the extent of the “crater” diminishes, and finally at po = 20 © 
psig no depression is discernible. Although it was difficult — 
to determine precisely at what value of the pressure the ice — 
ceased to melt around the stagnation point, the experimental — 
results and the simple theoretical prediction seemed to agree 
quite well.4 
4 A similar experiment was carried out at M = 8, 7) = 900 F, 

and the result confirmed this prediction. 
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Fig. 3 Nose position vs. time as measured on photographic 


film 
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Fig. 4 Photographs of ice models in 1 o = 5.8 flow 


Temperature on Cylindrical Afterbody 


The thermocouples were installed on the outside surface of 
the micarta tube, and they were 0.125 in. below the initial 
surface of the ice. The temperature gradients normal to the 
surface in the ice were inferred from the measured tempera- 
ture drop across the micarta tube wall, and they were insig- 
nificantly small in the region not too close to the nose. Hence 
it was concluded that the temperatures on the micarta tube 
were equal to the surface temperature of the ice. 

When the ice models were frozen and installed in the test 
section, the temperatures were about —7 C. After the test 
section was evacuated to a pressure of about 1 psia prior to 


p, = 20 psig, T, = 300°F 


the start of a test, the temperatures dropped to about —17 C 
because of the increased evaporation at the lower ambient 
pressure. Immediately after the flow was started in the test 
section, the ice temperature decreased and reached steady- 
state temperatures in the range —30 to —40 C, 3 min after 
the flow started. Then the temperatures remained un- 
changed until the thermocouples became exposed to the air. 
Thus the model actually became cooler in spite of the aero- 
dynamic heat transfer, because the equilibrium surface tem- 
perature was determined mainly by the balance between this 
heat transfer and the heat carried away from the surface by 
ablation. The measured ice temperatures are plotted in 


2 


m Le ,cal./cm. sec. 


° 


X/D 
Fig. 6 Ablation heat transfer rate along cylindrical surface 
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Ablation Rates 


The measured ablation rate at the stagnation point of the 
ice model was (1.10 + 0.17) X 10-8 cm/sec at po = 2.0 atm 
and 7) = 149C. The corresponding value of the vaporization 
heat is 0.69 + 0.11 cal/em? sec. Under these conditions, 
the theoretical value of the heat transfer without mass in- 
jection is 0.81 cal/cm? sec at the stagnation point of a 1-in. 
diameter sphere. The theoretical value of B’ is 0.05, and 
the heat conduction in the ice at the stagnation point is 0.05 
cal/em? sec. Therefore, the theoretical value of the heat 
available for the evaporation becomes 0.73 cal/cm? sec, which 
compares favorably with the measured value. 

Although the individual measurements can be made with 
high accuracy, the scatter of the data from several tests is 
about +15 per cent. It was noticed during the experiment 
that some models were ablated nonsymmetrically around the 
nose, and the results from obviously defective models were 
discarded. This behavior may have been caused by the 
difficulty of making ice models with the same consistency and 
homogeneity. 

The rate of heat absorption by sublimation along the 
cylindrical afterbody is shown in Fig. 6, together with the 
theoretical values. The large part of the scatter in the 
measured values was caused by the surface roughness pro- 
duced by surface melting during the installation and starting 
period. 


Conclusions 


1 Calculations based on the available laminar boundary 
theory indicate that experiments performed in moderate 
temperature wind tunnels yield basic results useful for the 
understanding of aerodynamic ablation. With judicious 
choice of materials a fairly wide range of parameters may be 
covered. 

2 The experiment performed with water ice showed fair 
agreement with the theoretical prediction. Experiments 
with other materials such as dry ice and camphor are desir- 
able to study the effects of various physical properties of 
material and larger rates of mass addition on ablation. 
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Nomenclature 
Cx = thermal enthalpy transfer coefficient 
Cp = specific heat at constant pressure 
co = specific heat of the solid body 
Di. = coefficient of binary diffusion 
= static enthalpy of mixture 
hk, = static enthalpy of vapor 
= static enthalpy of air 
h, = total enthalpy = h + (u?/2) 
k = coefficient of thermal conductivity : 
K = mass fraction of vapor of the subliming material _ : 
Le = pDyc,/k, Lewis number 
Le = specific heat of sublimation 
M Mach number 
M; = molecular weight of 7th species 
m = rate of ablation “45 6* 
p = static pressure 
Po = tunnel supply pressure 
= vapor pressure 
Pr = uc,/k, Prandtl number ar 
q = rate of transfer of enthalpy 
qo = rate of heat transfer to interior of solid 
T = temperature 
= tunnel supply temperature 
u,v = components of velocity parallel and normal to surface 
x = distance along axis of cylindrical body wy 
y = distance normal to surface - 
= absolute viscosity 
p = gas density : 
()e = values at the edge of boundary layer 
( )w = values at the solid surface 
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Technical Notes 


Launch Parameters for Interplanetary 
Flights 


W. C. RIDDELL! 


Convair (Astronautics) Div., General Dynamics Corp., 
San Diego, Calif. 


T IS assumed that the heliocentric transfer orbit from 

Earth to the target planet, both considered as massless 
points, is known. This may be determined by using two- 
body equations to describe transfer ellipses and an iterative 
process to match flight time with planet positions. The 
initial heliocentric velocity vector, then, is known. This is 
transformed to a geocentric vector by subtracting the ve- 
locity vector of Earth. The direction of this vector is called 
the target direction. The magnitude is the hyperbolic excess 
of the departure orbit, and the square of the magnitude is 
taken to be the energy of this orbit. Intervals of time and 
distance which occur between launch and the practical 
attainment of the vector are small compared to the inter- 
planetary distances and times involved. 

The development is carried out with the idea that there will 
be no powered change of plane. Hence, the missile plane 
relative to a nonrotating Earth is defined as that plane which 
contains the Earth-centered target direction vector and the 
launch site position vector at the instant of launch.? For 
any time of day such a plane is possible. With each plane a 
launch azimuth is associated. The relationship between 
launch azimuth and launch hour is obtained. 

If departure is made from a circular coasting orbit, the 
point of entering the hyperbolic orbit is easily found, since, 
for tangential impulse, it will be the perigee of the hyperbola. 
The coasting orbit should be circular for greater efficiency. 
If the coasting orbit is elliptic, iteration is needed, since the 
elements of the hyperbola depend upon the point of departure 
which, in turn, depends upon the required energy and asymp- 
tote direction of the hyperbola. For the direct shot, itera- 
tion may again be needed, since the launch time, and there- 
fore launch azimuth, will depend upon the burnout conditions 
of downrange, altitude and flight path, and these in turn de- 
pend to some extent upon the launch azimuth. 


Method 


The initial velocity vector of the heliocentric transfer orbit: 
is referred to the center of Earth by subtracting Earth’s 
velocity vector. Its direction now is the target direction 
from Earth. The expression of the vector may take various 
forms. A typical and convenient sequence of steps is the 
following 


Received March 11, 1960. 
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2 This neglects the generally unavoidable deviations from this 
plane that occur during the powered launch phase and which are 
Sass of Earth’s rotation, guidance commands, aerodynamic 
effects, etc. 


= (VsinT — V, sin I.) cos w — 
(V cos — V, cos I.) sin u 

v, = (VsinT — V, sin I.) sin w + 
(V cos cos I — V, cos cos u 
v, = VeosI'sin] 


v = (V?+ V2 — 2VV, (sin sin I, + 


v 
1 
vy cos € — v, sin € Jequatorial 
. 
Vy SIN € + COS € 


cos A cos 6 
l, = {| sin A cos 6 


sin 6 


Equating Equations [2 and 3] gives the right ascension and 
declination of the target direction. 

A launch into a plane containing the target vector may be 
made at any time of day, the launch azimuth varying with 
launch hour. The relation between azimuth and launch hour 
may be obtained by equating the expression of the normal to a 
plane passing through the launch site with a given azimuth 
to that for a plane passing through the launch site and the 
target vector. Itis 

sin (A — H) cos 6 
cos a = = 
sin 
sin (A — H) cos 6 
{1 — [eos (A — H) cos 1 cos 6 + sin 1 sin 6]?} 1/2 


= 


[4] 


Letting a = 0 gives 
cos (A — H) = tan 6/tan 1 [5] 


Thus when |/| > |6| a due east launch may be made. 
In case |l| < |6|, the least |a| is found, from spherical 
trigonometry, to be obtained when 


cos (A — H) = tan l/tan 6 i (6) 
and then 
cos a = cos 6/cos | [7] 


In both cases the plane is that having the smallest inclina- 
tion to the Equator. In the first case the inclination is equal 
to the latitude of the launch site, and in the second case to 
the declination of the target vector. The launches would be 
into coasting orbits, since it is highly unlikely that other 
conditions would be right for a direct entry into the hyper- 
bola. 

The point of initiating the hyperbola may be found. The 
true anomaly of the asymptote is given by 

cos = —l/e me = —(e? — 1)!/? 


For the circular coasting orbit the injection point is the 


Eprror’s Note: The Technical Notes and Technical Comments sections of ARS JourRNAL are open to short manuscripts describing 
new developments or offering comments on papers previously published. Such manuscripts are usually published without editorial! 
review within a few months of the date of receipt. Requirements as to style are the same as for regular contributions (see masthead page). 
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perigee and thus is 7. back from the target vector. 
case 


= (3 2)" [8] 


In case an elliptic coasting orbit occurs, a tentative initial 
point is chosen. The flight path angle and radius of this 
point are determined from the orbital elements as found 
from the burnout conditions of the surface launch phase. 
Now 


The true anomaly in the hyperbola of the initial point is found 
from 


= v y2 1/2 
tan = cos +2) (9) 


(v?/v.2 + 2) sin 27; 
a 2[(v2/v.2 + 2) cos? y; — 1] 


If the selected point is at an angle wu in the forward direction 
from the target vector, then it required that 


+ (Na = 360 deg 


The process is repeated until a point is found for which this 
occurs. 

In the case of the direct shot, kinematics and geometry are 
interrelated. The angle (ns — 7) is found from Equations 
{9 and 10]. Adding to this the downrange angle of the power 
phase gives @. H and a are then found from 


cos (A — H) = 
cos a = sin (A — #) cos 6/sin d 


[10] 


tan 7: = 


(cos @ — sin! sin 6)/cos/ cos 6 [11] 


Since launch is from a rotating Earth, the power phase is 
affected by a, so that a cyclic iteration may be needed. Also, 
Equation [11] may give an imaginery (A—#), so that the 
power phase may have to be revised to obtain a compatible 
¢. A low ¥; will usually require a large |a|. 

The hyperbolic injection velocity for all cases is given by 


v2 1/2 
(= 2) [12] 


The launch azimuth relative to a rotating Earth may be 
approximated by 


Uo/Ve 
cot a’ = [13] 


The distance of the asymptote from the target vector is 
given by 


p = + 2/(v/v.)*}” [14] 


As a simple example let w = 0, T = I, = 0,7 = 1.5 deg, 
V/V, = 1.187, € = 23.5 deg, 1 = 28.5 deg, and let launch be 
into a 100-nautical mile altitude circular coasting orbit. 


From Equation [2] 


= (0.0819/0. 1136) = 35.8 deg 


0 
1, = “(0.1136 
 \0.0819 


Using Equation [3] 
A = 90 deg 


+41.1 deg AD 
$90 deg + 41.1 deg = 131. 1 
(90 deg — 41.1 deg = 48.9 deg 


hr 
3.26 hr sidereal ae 


= +22.8 deg 


From Equations [6 and 7] we find 
(A—H) 


Equation [13] gives a’ = +24 deg. For minimum |{a’| 
then, launch should be made at approximately 3.26 hr 
sidereal and 24 deg north of east, or at 8.73 hr sidereal and 
24 deg south of east. 

Using Equation [lb], v/v. is equal to 0.535, and from 
Equation [8], 7. is equal to 141.1 deg. The perigee, there- 
fore, lies 141.1 deg back from J, in the missile plane. At 
this point, from Equation [12] 


= 1.514 


From Equation [11] it is found that ¢ = 
angle from launch to hyperbolic initiation is 


35.3 deg. The 


360 — 141.1 + 35.3 = 254.2 deg a 
for the northeast launch, or 

360 — 141.1 — 35.3 = 183.6 deg V« 
for the southeast launch. =... 


From Equation [14]. p is about 10,000 nautical = or 
some 3 Earth-radii. 


Nomenclature 

l, = target direction vector, Earth centered 

v = velocity relative to Earth at initiation of planetary transfer, 
i.e., the hyperbolic excess velocity 

V = velocity relative to the sun at initiation of planetary trans- 
fer 

V. = orbital velocity of Earth at time of initiation of planetary 
transfer 


lr = initial flight path angle of interplanetary orbit 

I. = path angle of Earth’s orbit at initiation of interplanetary 
orbit 

I = inclination to the ecliptic plane of the interplanetary orbit 

w = longitude of Earth at the initiation of the interplanetary 
orbit 

e = inclination of the Equator 

A = right ascension of the target vector 

6 = declination of the target vector 

a = launch azimuth relative to a nonrotating Earth, measured 
from due east, positive to north 

‘ = launch azimuth, but relative to a rotating Earth 

H = sidereal hour angle of the launch site at instant of launch 

1 = latitude of launch site 

nq = true anomaly of the hyperbolic asymptote 


e = eccentricity 
ve = local circular velocity 
vo = Earth surface velocity 

y = flight path angle 


: = angle in the missile plane from launch site to target vector 
= distance from target vector to hyperbolic asymptote 
Subscript i irefers to the initiation of the hyperbola 
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Effect of Earth’s Obie on the 
Calculation of the Impact Point of 
Ballistic Missiles 


CHING-SHENG WU! 


Jet Propulsion Laboratory, 
California Institute of Technology, Pasadena, Calif. 


HE EFFECT of Earth’s oblateness on the calculation of 

the trajectory of long-range ballistic missiles may be con- 
sidered from two different aspects: The dynamic effect and 
the geometric effect. The former comes from the non- 
Newtonian gravitational force which has a potential (1)? 


_ GM [Ro {1 


= constant of gravitation 
M = Earth’s mass 

Ro = Earth’s equatorial radius 
r = radius measured from the center of Earth 
@ = latitude 

J = constant, 1.638 1073 


Generally speaking, this dynamic effect is quite small and is, 
therefore, neglected in most of the existing analyses concern- 
ing missile trajectories. The second effect is a result of the 
nonsphericity of Earth. This effect may be of interest in the 
calculation of the impact point, the point where the missile 
trajectory intersects Earth’s surface. 

In the following discussions, the dyamic effect will be ig- 
nored. The geometric effect, however, will be taken into 
account in order to generalize the analy sis given by Singer 
and Wentworth (2). 

To calculate the impact point, we shall consider the entire 
trajectory as having two successive stages, namely the ballistic 
phase and the re-entry phase. The ballistic phase starts 
from the burnout point and terminates at a hypothetical 
“re-entry point.’ In this phase, the missile will move along 
an elliptical trajectory whose geometrical configuration de- 
pends entirely upon the magnitude and direction of the burn- 
out velocity and the burnout altitude. The aerodynamic 
drag is neglected during the calculation of the ballistic flight. 
However in the re-entry phase, the performance characteristics 
of the missile will be entirely different. After the missile 
re-enters Earth’s atmosphere, it will experience tremendous 
aerodynamic friction because of its extremely high flight 
speed. As a consequence, aerodynamic heating will be a 
serious problem. [This point will not be discussed here; 
further information may be obtained from (8).] In the 
analysis of the re-entry phase, the terminal conditions of the 
ballistic phase may be used as initial conditions. Again, 
since the time of flight during the re-entry phase is short com- 
pared to the time of the ballistic phase, the rotation of Earth 
may be neglected. The calculation of the re-entry phase 
trajectory can be readily performed by applying the analysis 
of Allen and Eggers (3). 


Analysis 


If we consider Earth to be nonspherical and the atmospheri- 
cal shell similarly oblate, the effect of nonsphericity on the 
determination of missile impact point may be seen to be of 
interest only in the ballistic phases. This follows because the 
oblateness of the atmospherical shell will shift the re-entry 


point of the rocket when the altitude of the atmospheric shell 


is fixed. 
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_ missile will be determined if the burnout velocity vector 


According to Singer and Wentworth (2) 


+ [r+é + - sin? cos*( gz — |= 


The elliptical trajectory during the ballistic phase of the 


and burnout altitude of the rocket are given. The trajectory 


__ may be described by a polar coordinate system designated 
__ within the plane of the missile trajectory with origin at the 


center of the Earth as shown in Fig. 1. Again, if p is the 


} polar distance, the ballistic trajectory is given by the well- 


known equation 
= a(1 — e*)/(1 — e cos ¢) oa 
where ~ a 
¢ = polar angle measured from the perigee (see Fig. 1) 
@ = semimajor axis of the orbit ellipse 
e = eccentricity of the orbit ellipse = (1 — 6?/a?)!/2 
b = semiminor axis of the orbit ellipse 


Wag (2GM/Rx — | 


(GM)? 


where dow ses = 
= burnout velocity 
R, = radial distance from the center of Earth and the 
burnout point 
G = gravitational constant 
M = mass of Earth 


The international ellipsoid of Earth may be represented 
also in the terrestrial coordinates 


(1 — €*)(r? sin? ¢ — Ro?) + r? cos? g = 0 [2] 


where 
€ = (equatorial radius — polar radius) /equatorial radius 
go = 7/2 — latitude 
Ro = equatorial radius 


Hence, for a given latitude ¢, the local radius of Earth may 
be determined from 


= (1 — — — ¢]!/? [3] 
with « ~ 1/297. Since ¢€ is small, it is possible to write 
r = [1 — (¢/2)cos* [4] 


where = €(2 — €). 

If the orientation of the trajectory plane and the orienta- 
tion of the equatorial plane are different by an angle ®o, then 
the relation between the angle ¢ and the angle of advance 
of the missile 8 may be determined such that (4) 

cos? = sin? cos?B [5] 


where @ is measured from a plane containing the North Pole 
and Earth’s center, but perpendicular to the plane of tra- 
jectory (see Fig. 2). Since it is possible to write 


B= 


where Apo is known if the perigee of the “orbit ellipse”’ is fixed, 
we therefore have 


= [1 — (€/2)sin? Bp cos?(g — Ao) [6] 


The angle gg at the re-entry point may be determined by 
equating Equations [1 and 6], i.e., by solving 


a (1-e) 


Ro 1 + cos 7] 


where / is the thickness of the atmospherical shell. 
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Fig. 1 Elliptical trajectory of a missile 


If the difference of altitude between the burnout point and 
re-entry point for a spherical atmospheric shell with radius 
(Ro + h) is denoted by k then 


k = Ri — (Ro + h) 


Thus after taking into account the oblateness of Earth, the 
difference of altitude between the burnout point and re- 
entry point k becomes 


k = k + Ro[(é/2)sin? Bp cos*(gz — Ao) ] 


Since in general R is small compared to Rz, we can obtain 
Oz (5) the angle between the trajectory and radial line (as 
shown in Fig. 2) at the re-entry point from the relation 


sin*®p cos*(yz — Xo) |x 


T k é 
4 


cos? 6[1 + tan 6 — cos? 


Again the velocity at re-entry point vg and the total time 
of flight during the ballistic phase tg can be determined as 
follows 


5 sin? dy cos*( gz — ra) |x 


GM 
[1 + sin? sin 20 (1 cos? 5) | 


k Roz. 
te =t+ E 5 sin? cos*( |x 


T.J/1— ele + cos ¢(Rz)] 1 


Qn + ecos 1+e” 
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Fig. 2 Geometrical relation between the trajectory plane and 


sin 


2 tan (1/2)(Rx)\]\ , 
ea | 


l+e f 


T, = wd 


After @z, ve and Rez are known, the calculation of the re- 
entry phase may be readily started. Discussion of the pro- 
cedure is not considered here. 


Remarks 


It is seen that the effect of oblateness is generally small, es- 
pecially for small ®. However, according to a rough esti- 
mation, the spherical assumption may introduce an error 
about 50 km in the calculation of impact point when a = 
1.1Ro, e = 0.4 and ®) = Ao = 7/2 are assumed. 

Again, since the coordinates are assumed to be fixed in 
space, in order to determine the angle ®, the local linear 
velocity at burnout point owing to Earth’s rotation should 
be considered. In order to calculate the correct latitude of 
the re-entry point, the rotation of Earth must be taken into 
account. This may be easily calculated, since the total time 
of flight is known. Since no new feature will appear in this 
case, discussion along this line is omitted. 
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Spin Axis Attitude 


Radio Corp. of America, Princeton, N. J. 


HE GRAVITATIONAL torque acting on an Earth satel- 
lite has been analyzed in numerous publications. The 
long-term effects of differential gravity on a satellite have 
been discussed by Doolin (1)? and Manger (2). However, 
Doolin’s paper contains results for only a few orbits, and 
Manger’s uses an approximation which averages the gravity 
torque over the orbital period. This note reports the results 
of a computer study the specific purpose of which was to 
determine accurately the secular variation in spin axis atti- 
tude resulting from gravitational torque. The interesting 
dynamics revealed by this brief study would seem to warrant Fig. 1 Orbital geometry 
an extensive analysis of the problem for any “spin stabilized”’ : 
satellite application. 
The gravitational torque (2,3) is taken as 


M = X 5) [1] 
If the total angular momentum vector of the satellite be 


considered as coinciding with the spin vector @, the vector 
equation to be solved is 
Tw(d3/dt) = 3a? Al (7-5) (F X 5) [2] 
The orbital geometry is shown in Fig. 1. 
Equation [2] is then expanded to give 


ds;/dt = k(rzs8, + 1y8y + 1282) (Ty8: — [3] 
= + ry8y + 148.) (tebe — 1282) [4] 
ds,/dt = k(r,s, + + r8,) [5] 


cos 2 cos 8 + sin Q sin i COS a [6] 


ll 


sin 2 cos 8 — cos Q sin B cos a [7] 
rT, = sin B sin a [8] 


Equations [3, 4 and 5] were solved computing 36 points per 
orbit while holding 2 constant for 10 orbits, and then ad- 
vancing 2 by the appropriate AQ for computing the next 10 
orbits. The computation was carried out for the following 
conditions: 

Orbital altitude = 1000 miles 

Spin angular velocity = 10 and 100 rpm 

Total weight of satellite = 150 Ib 

Orbital inclination = 66.5 deg 


k= 7.365 X 10-7 


A direct orbit was used, with initial conditions selected to pro- 
vide for the initial attitude of the spin axis to be normal to the 
ecliptic plane. 

Of particular interest for this study was the sun angle A 
defined in Fig. 2 as 


A = cosR-$ 
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Fig. 4 Projection of spin axis unit vector on ecliptic plane— 
100 rpm 


The components of the spin axis unit vector 5 referred to 
the inertial (sun-centered) coordinate system, may be written 
as 


Sx= 8, 
8, cos 23.5° — s, sin 23.5° 


8, sin 23.5° + s, cos 23.5° 


The sun angle A then becomes 


A = cos7! (8, cos wet + s, sin wet cos 23.5° — 
8, SIN wet sin 23.5°) 


The results are plotted in Figs. 3 through 6. It is clear 
that the spin axis of a spin stabilized satellite is not so stable 
as at first appears, unless the spin rate be sufficiently high. 
In any case, the secular variation in spin axis attitude must 
certainly be accounted for. 
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Fig. 5 Sun angle for 10-rpm spin 
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Fig. 6 Sun angle for 100-rpm spin 


7 = unit vector along the radius from Earth’s center 
Is = spin axis moment of inertia ; 
w = spin angular velocity 
2 = position of ascending node 
B = orbital position of satellite = 
x,y,z = Earth-fixed coordinate system 
X,Y,Z = inertial (sun-centered) coordinate system 
R = unit vector along radius from the sun to Earth 
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Stabilizing a Spherical Satellite by 


Applied Physics Laboratory, 
The Johns Hopkins University, Silver Spring, Md. 


ROBERT R. NEWTON? 


have been many proposals for using radiation 
pressure for attitude stabilization. In two recent pro- 
posals, Sohn (1)* discusses a weathervane-type of stabilizer, 
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and Frye and Stearns (2) discuss the use of a large cone ap- 
pended to the main satellite. Both of these involve awkward 
geometry. This paper shows that stabilization satisfactory 
for many purposes can be achieved without departing appre- 
ciably from spherical geometry. 

The success of a stabilization method depends upon the 
torques which must be stabilized against. If the principal 
upsetting torque is gravitational, the magnitude of this 
torque (8) is 


(89/R) I xz [1] 
where 
R= distance of the satellite from Earth’s center 
g = gravitational acceleration at that distance : 
Ixz = product of inertia 


We assume that / xz will not exceed 10~% of a principal mo- 
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RADIATION 


Hf Fig. 1 Absorption on any surface 


(b) 


Fig. 2 Torque on a sphere, of which one half is reflecting and 

one half is absorbing. (a) Viewed normal to the radiation; (b) 

viewed parallel to the radiation. The torque arises from the ab- 
sorption on the cross-hatched area 


Fig. 3 Torque on a sphere with a girdle which shadows the rear 


hemisphere. The net torque arises from absorption on the 
cross-hatched area 


(a) 


Fig. 4 Torque on a sphere with a girdle which does not entirely 
shadow the rear hemisphere, assumed to be absorbing. The 
net torque arises from absorption on the cross-hatched areas 


7 


Fig. 5 Dimensionless torque 37'/2Pa* as a function of @ for: 
(a) a sphere with an absorbing hemisphere; (b) a reflecting 
sphere and a girdle with a,/a = 1.1; (c) a reflecting sphere and 


a girdle with a,/a = 1.2; (d) a pessimistic estimate of maximum i 


gravitational torque 


ment of inertia. For a hollow spherical shell, of 100-kg 
mass and 0.5-m radius, Jxz < 107% (2/3)(100)(0.5)? = 
(5/3) X 10-*kgm?. If R = 7.38 X 10°m (an altitude of 1000 
km), we get from Equation [1] 


= 5 X at:.m. [2] 


Before calculating the torque from radiation pressure, we 
shall first make three general points. 

1 With spherical geometry, no torque arises from external 
radiation pressure on a perfectly reflecting surface. With 
perfect reflection, the resultant force is normal to the surface; 
hence it passes through the center of the sphere and exerts 
no torque. 

2 With spherical geometry, no torque arises from the 
pressure of emitted radiation. The force on any surface 
element from emitted radiation is normal to the element and 
hence exerts no torque. With nonspherical geometry, it is 
just as necessary to consider the emitted as the absorbed 
radiation; over a long period of time, just as much energy 
is emitted as absorbed. 

3 The torque arising from total absorption on any surface 
is the same as that arising from total absorption in the shadow 
of the surface on a plane normal to the radiation. From Fig. 
1, we see that the same energy and momentum would be ab- 
sorbed by either A; or Ay. Hence the forces on A; and A: 
are equal, have the same line of action, and exert the same 
torque about any point. 

We now look at three important configurations: 

(a) In the first configuration, we imagine a sphere of 
radius a, with one hemisphere totally absorbing and one 
totally reflecting (see Fig. 2). By virtue of the general state- 
ments, the only torque present arises from that part of the 
absorbing surface which lies in the hemisphere toward the 
sun. The direction of this torque is such as to rotate the 
satellite clockwise as drawn in the figure, so, in the equilibrium 
attitude, the reflecting hemisphere faces the sun. Further, 
as indicated in the figure, the torque is the same as that pro- 
duced by absorption in the lunate area lying between the 
semicircle of radius a and the semi-ellipse having semi-axes 
aand acos 6. 

(b) In the second configuration, we add to the sphere of 
configuration (a) a girdle, whose outer radius is a, of absorb- 
ing material, lying in the plane which divides the reflecting 
and absorbing halves of the sphere. The absorbing material 
in the girdle is placed so that it all faces the same direction 
as the reflecting hemisphere (see Fig. 3). Further, we assume 
that the sphere is not tilted far enough for the shadow of the 
sphere to intersect the outer edge of the girdle. 

From Fig. 3, we see that the torques from the parts of the 
girdle which are not cross-hatched balance each other. Fur- 
ther, the absorbing hemisphere is completely shadowed and 
contributes no torque. Hence, the only unbalanced torque 
arises from the part of the girdle which is cross-hatched, and 
the torque is the same as in configuration (a). 

(c) This is the same configuration as (b), except that the 
shadow of the sphere intersects the outer rim of the girdle 
(see Fig. 4). The unbalanced torque now arises from the 
cross-hatched part of the absorbing hemisphere, and the total 
torque is the same as before. 

We see that the torque is the same for all these configura- 
tions, and is the same as that from the lunate region shown 
in Figs. 2-4. This torque is therefore 


T=PS% defi y dy [3] 


where P is the radiation pressure. y- and y., the values of y 
on the circle and ellipse, are 


Yo = Va? — x? ye = cos 0 Va? — 2? [4] 
Substituting Equations [4] into Equation [3] and integrating 
T (2/3)P a’ sin? 6 [5] 
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The radiation pressure P estimated by Roberson (4) is 
P = 0.43 X 1075 nt/m? [6] 
For the satellite used in obtaining Equation [2] pet - 


T = 3.6 sin* 6 X 1077 [7] 


This torque will balance the maximum gravitational torque 
given by Equation [2] when the angle 6 equals 


6 = sin“(5 X 10-*/3.6 X 10~7)!/2 = 22 deg [8] 


This angle is satisfactory if the purpose of stabilization is to 
secure efficient operation of solar cells. The effective area 
of the solar cells at this angle is cos 22 deg = 0.93 of the total 
area. 

We conclude that radiation pressure is adequate to over- 
come the destabilizing gravitational torque, if we use rea- 
sonable care in the mass balance of the satellite. 

The configurations just discussed may present a serious 
heat balance problem: They present a reflecting face to the 
sun, and a “black” face away from the sun. They may be a 
material which can overcome this by having high absorption 
for the solar spectrum and low emission for thermal] radiation, 
to use for the back hemisphere. If not, we can probably 
get by with using the girdle only, and having high reflec- 
tivity over the entire sphere. This, however, decreases the 
torque available for large values of 6. 

As long as the shadow of the sphere does not intersect the 
outer edge of the girdle, as in Fig. 3, the rear hemisphere is 
completely shadowed, and its absorption is immaterial. 
However, when @ > cos~!(a/a,), only the cross-hatched part 
of the girdle gives torque; the sphere no longer gives torque 
if it is reflecting all over. Calculating the torque by the 
method already used 


(5) Pa’ sin? @ 


ay 
(;) Pa’ sin sin? 8 | cos 


The factor in braces vanishes when 6 = 90 deg, so there is a 
maximum torque available. Torque as a function of @ is 
shown in Fig. 5 for a;/a = 1.1 and 1.2, and also for the con- 
figuration having an absorbing rear hemisphere. The gravi- 
tational torque given by the equality in Equation [2] is 
shown as the horizontal line. It may be safe to operate with 
a,/a = 1.1, but more detailed estimates would be required 
to decide this. 

A new problem is introduced when the entire sphere is re- 
flecting. If, by some accident, the satellite should ever be 
facing the wrong way, it could not right itself. Hence, a 
compromise is needed. There must be enough absorbing 
material around the pole of the rear hemisphere to allow 
righting the satellite if it should face the wrong way, but not 
enough to cool it unduly. The amount of material needed 
will require careful computation. It will also be necessary 
that the back face of the girdle have low absorption. 
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Some Possible Operations on the Moon 


F. ZWICKY' 


California Institute of Technology, Pasadena, Calif. 


LASSICAL astronomy concerned itself with two basic 

activities: The observation of celestial bodies and of 
events by means of Earth-bound instruments, and the inter- 
pretation of these observations in terms of the scientific 
knowledge which has been accumulated on Earth. 

The exploration of extraterrestrial space and of its material 
contents was masterminded by men like Aristarchus (about 
265 B.C.) and Hipparchus (190-125 B.C.), Giordano Bruno 
(1548-1600) and Knut Lundmark (1889-1958), who con- 
ceived of the ever widening universes of the Solar System, 
the Milky Way stellar system (our own galaxy) and the realm 
of the extragalactic stellar systems. 

Because of the development of powerful, active radar and 
radio transmitters, and particularly because of the develop- 
ment of rockets and of other spacecraft, modern astronomy 
intends to embrace two additional activities: The direct 
experimentation with extraterrestrial bodies and phenomena, 
and the possible reconstruction of celestial bodies, starting 
with those of the Solar System. wer 
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Some suggestions as to the goals and operations involved 
have been given in other places (1).2 A more complete 
survey of this future field of endeavor is badly needed. In 
the present note we restrict ourselves to some considerations 
concerning the moon. 


Conditions on the moon 


Because of its lack of atmosphere, its temperature char- 
acteristics and its exposure to various types of dangerous 
corpuscular radiations, the moon appears as a most inhos- 
pitable place. 

On the other hand, the availability of an extremely high 
vacuum, as well as of sharp differences in temperature in 
front and in back of suitable baffles, which shade the radi- 
ations from the sun, will make possible the operation, with 
very high efficiency, of all sorts of thermomechanical and 
thermoelectric devices for the generation of mechanical and 
of electric power, as well as for the production of various 
types of most useful materials, both through the use of dis- 
sociation or of synthesis of the moon’s various natural con- 
stituents. From a generalization of our knowledge gained 
about the composition of the Earth, of the sun, of the at- 
mospheres of some of the planets, of the meteors, meteorites 
and comets, of the stars and of interstellar dispersed matter, 
we may assume with confidence that all of the essential 
chemical elements will be found on the moon. Some, of 
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course, like hydrogen in the water of crystallization of the 
various minerals, may have “‘evaporated’’ out of the surface 
layers and may have to be looked for at some moderate 
depths. 

We do not yet know the exact physico-chemical structure 
and composition of the moon’s surface layers, but it is hoped 
that some direct tests of the type described further on will be 
made in the near future. 

The moon offers another advantage, that is its low gravity, 
which is only about one sixth of that on Earth. This means 
that an astronaut will be able to erect large structures on the 
moon with relative ease, structures with far lighter construc- 
tion than on Earth. While thinking about constructions on 
the moon, attention should also be called to the intriguing fact 
that, because of the availability of an extreme vacuum, we 
shall be able to produce chemical elements and compounds of 
a purity not attainable in our laboratories. This will not 
only be of extreme interest to chemists, physicists and metal- 
lurgists, but, as Professor Pol Duwez of the California Insti- 
tute of Technology has observed, it will make possible entirely 
unheard of new means of construction. Indeed, in the ab- 
sence of an almost instantaneous formation of protective layers 
of oxides, malleable metals like aluminum, magnesium, etc., 
can be welded, without heating or other auxiliary means, 
by simply compacting them together. 
Goals to be achieved 

Astronauts who intend to “live off the land’’ on the moon, 
must produce out of its constituents oxygen and nitrogen for 
breathing, water and certain elementary foodstuffs for living, 
shelters for protection and mechanical and electric power 
for all of the desired activities of construction and of locomo- 
tion. 


Required and available power sources 


The simplest and most readily available power source is the 
radiation from the sun. We shall here be mainly concerned 
with its use. At some future date, when the problem of the 
ignition of nuclear fusion reactions will have found universal 
and easily applicable solutions, the use of solar energy on the 
moon may become less important. 

If we use an optical image furnace to concentrate the solar 
radiation, we shall have available a maximum flux of energy 
A X o, where A is the aperture of the lens or mirror used and 
o = 1.4 X 10% ergs per cm? sec = 1/30 cal per cm? sec. 

The maximum requirements for energy stem from our need 
to liberate O2, Ne and other elements from compounds such 
as H,O and CO:, whose heats of formation are respectively 
68.3 and 94.5 keal per gm mole. For 1 gm mole of O2 we 
thus require energies of the order of 100 keal. Since a mirror 
of 1-m? aperture collects a radiation flux of 1/3 keal per sec, 
32 gm of oxygen or 22.4 liters at standard pressure and tem- 
perature can be obtained in 300 sec. This amounts to about 
270 liters per hr, which covers the need for about 15 persons, 
since the amount of oxygen used up in breathing, per indi- 
vidual and hour, is about 15 liters at SPT or 21.4 gm. 

Nitrogen, which, for reasons as yet unknown to the present 
author, seems to be essential as an admixture to oxygen for 
sustained breathing, can be generated in sufficient quantities 
once and for all at the beginning of the operations on the 
moon. If the astronaut is lucky and finds the proper min- 
erals, he will be able to generate 1 gm mole of either oxygen 
or nitrogen at an expense of energy smaller than 100 kcal. 


Degrees of dissociation at solar furnace temperatures 


From the considerations made, it follows that small solar 
furnaces of about 1-m? aperture intercept and concentrate 
an ample flux of radiation from the sun to produce oxygen, 
nitrogen and water in sufficient quantities to sustain the life 
of several astronauts. It is also certain that the available 
energy can be trapped at an absolute temperature 7 which is 
actually elevated enough to decompose into the desired ele- 
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ments almost any of the rocks or sands which make up the 
surface layers of the moon. Theoretically, T could be as high 
as the temperature of the sun, that is about 5700 K. For 
practical purposes, however, it is quite sufficient to achieve 
temperatures 7” of the order of 4000 to 5000 K. In the sun’s 
atmosphere most molecules are dissociated, and some of the 
remaining molecules, like TiO, exist only in very small rela- 
tive concentrations. In fact, not only is the concentration 
of molecules small, but most of the atoms with ionization 
potentials less than about 9 ev are ionized. Only hydrogen, 
the noble gases and other atoms with higher ionization po- 


tentials exist in the sun’s photosphere as neutral atoms. 7 q | 


Rates of dissociation and of evaporation 


In order to complete our inquiry into the production of 
oxygen, nitrogen, water and other necessities for sustaining 
astronauts on the moon, we must inquire into the rates, that is 
into the kinetics of the processes of sublimation, evaporation 
and dissociation of the rocks on the moon. Actually, this 
is a field which requires further research, although it is cer- 
tain that, at temperatures between 2000 and 3000 K the 
water of crystallization will escape from the rocks rapidly, 
and carbon dioxide and nitrogen will also be produced at fast 
rates. It is also known from laboratory tests that the 
evaporation of solids like graphite, at a temperature above 
3000 K, is very “‘lively.”’ Nevertheless, concentrated re- 
search on rates of evaporation, sublimation, dissociation and 
ionization at temperatures above 2000 K is in order. 


Modes of operation of a solar furnace on the moon 


Low Temperature Operation: The water of crystalliza- 
tion in various minerals is generally not so tightly bound as 
are the atoms in the most stable molecules. This water can 
therefore be sublimated from the rocks at very much lower 
temperatures than can actually be achieved in a good solar 
furnace. We remark in this connection that the igneous 
rocks on Earth contain on the average around 1 per cent of 
water of crystallization, but that there may be as much as 
8 per cent contained in some minerals. Furthermore 1 mole 
of CaCO;, for instance, may be dissociated into CaO and 
CO, with an expenditure of only 43 keal. Or, as another 
typical case, 2 moles of Ca(NOz)2 can be dissociated into 
2CaO + 2N2 + 30, with an input of only 145 keal, thus fur- 
nishing (2N2 + 30.) gm moles or 44.8 liters of nitrogen and 
67.2 liters of oxygen. We may thus use a solar furnace of 
mediocre optical quality to produce temperatures of between 
2000 and 3000 K, and with it extract water, nitrogen and car- 
bon dioxide from the rocks on the moon. 


High temperature solar furnace 


At a next step we may aim to extract pure oxygen for 
breathing as well as perhaps pure hydrogen for purposes of 
obtaining fuels and propellants for internal combustion en- 
gines, rockets and for propulsive powerplants in general. 
Indeed, even at T’ as low as 3500 K, and atmospheric pres- 
sure, 85 per cent of the carbon dioxide will be dissociated into 
carbon monoxide and oxygen. At lower pressures this frac- 
tion is still greater. The fraction of H,O dissociated at a pres- 
sure of 1 atm and at 7’ = 3000 K is 0.164. There is thus 
no difficulty in producing water, oxygen and nitrogen. 


Generation of electrical and mechanical power 


Using optically high-quality solar furnaces, their beams 
can be focused within reaction chambers of, for instance, 
tungsten or tungsten carbide and titanium carbide, en- 
closing the minerals which are to be reacted. If these con- 
tainers have the shape of rocket chambers, gases at high speeds 
v will be ejected from their nozzles. These jets will contain 
positive ions of most of the metals, such as Al, Mg, Ca and Fe. 
These, when traveling through a magnetic field H, will be de- 
flected and separated from the electrons and negative hy- 
droxyl ions. From electrodes placed at right angles to both 
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vy and A, an electric current can therefore be drawn. The 
current thus generated by the ionic induction jet can be used 
as a secondary power source. It can also be used to increase 
the strength of the magnetic field originally applied to the 
jet induction generator, and, finally, it may be conducted 
through the melt in the reaction chamber to augment the 
production of ions through direct electrolysis. 

A thermomechanical power generator might be built in the 
shape of a closed cycle gas turbine, whose working gas is 
heated through heat exchange from a reaction chamber located 
in the focus of the solar furnace. The necessary condensers 
for the system could be placed underground, where the tem- 
perature will be sufficiently low. 


Auxiliary operations 


Among these are the production of “bricks” from the molten 
slags, and the welding together of stones and bricks through 
glazing of the joining surfaces by means of the focal beam of a 
primitive solar furnace. A temperature of 3000 K is ample 
for melting any of the metal oxides which are of interest in 
this connection. 

A further operation to be considered is the compression to 
high pressures of some of the gases produced out of the rocks. 
Such pressures could be achieved either directly through con- 
version of the momentum of the jets described in the fore- 
going, or by the more conventional means of pumps activated 
by the auxiliary electric and mechanical power available. 
This power is a fraction of the 1.4 kw of solar radiation flux 
intercepted per square-meter aperture of the solar furnaces 
used. 


Foodstuffs from the moon 


The production of foodstuffs from the materials of the moon 
should also be possible as soon as high pressures can be 
achieved. Operations would begin with such syntheses as 
CO + 2 H. = CH; OH, and then up the line, so that to the 
dismay of some, the astronauts would get their first internal 
calories from drinks of ethylalcohol diluted with water! 


Solar furnaces 


Summarizing, it may be said that the astronauts’ most 
immediate first needs could be filled by small solar furnaces 
of optically intermediate quality and capable of producing 
temperatures up to 3000 K. The optical system of a furnace 
might consist of a plane and a concave mirror. The latter 
could be spherical, if a simple (plastic) Schmidt correction 
plate were added. It would then be sufficient to hand crank 
the system every half hour or so to follow the course of the 
sun. The focusing of the rays could take place on the sur- 
face materials of the moon, and these be covered by a bell 
jar jammed into the ground, with piped outlets leading to 
buried collecting vessels. Water would condense in the first 
trap, whereas carbon dioxide and nitrogen might be separated 
by the ordinary means (temporary absorption or adsorption 
of CO), or CO. might conceivably freeze out at the tempera- 
tures encountered sufficiently far underground. 

The high temperature furnaces which are required for the 
dissociation of all minerals, as well as for the production of 
the maximum amount of ionization for the generation of elec- 
tric power, would have to be of better design than the medium 
temperature furnaces. 


Creation of an atmosphere for the moon and nuclear 
fusion 


While speculating on how to endow the moon with an 
atmosphere, the following thoughts occurred. Since the 
gravity on the moon is too feeble to retain an atmosphere of 
oxygen and of nitrogen in the face of the sweeping action of 
the sun’s radiation, the gravity would have to be increased 
if one insisted on an atmosphere which contains enough 
oxygen to sustain life. Increasing the mass of the moon for 
this purpose by a transfer of matter in a wholesale recon- 
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struction of the planetary system is thinkable. More feasi- 
ble, however, will be the shrinkage of the moon to about half 
its diameter through the ignition of nuclear fusion reactions 
within its interior, once such reactions can be handled on a 
large scale. 

With nuclear fusion reactions available, one may think. 
however, of endowing the moon, as it is, with an atmosphere 
which it can hold. This could, for instance, be an atmosphere 
of carbon dioxide. Such an atmosphere would act as a shield 
against dangerous extralunar radiations. Through suitable 
catalytic devices, CO. might be rendered breathable, or, if 
the lunar population lived in shelters, the atmosphere could 
be circulated through subsurface gardens in which plants 
would assimilate the carbon and release oxygen. 

All of these schemes depend on the solution of the problem 
of the nuclear fusion ignition of all important elements in the 
periodic system below iron. As is well known, far reaching 
attempts to solve a part of this problem through the investi- 
gation of the dynamics of electromagnetic plasmas are in 
progress. Attention should be called to the view that the pre- 
occupation with plasma physics may not be the only or even 
the most suitable way of searching for solutions of the prob- 
lem of nuclear fusion ignition. The first thought that 
emerges from a morphological study of the problem is that, 
instead of accelerating a stream of ions to very high speeds, 
one might attempt to do so with streams of neutral particles, 
or with droplets and with pellets of colloidal size, or perhaps 
even of greater size. On closer scrutiny, several thermo- 
mechanical and electromechanical modes of propulsion sug- 
gest themselves for the purpose of accelerating small pellets 
up to velocities of the order of 1000 km per sec. Pellets of 
the corresponding kinetic energy of 5 X 10" ergs per gm 
should, on impact with suitable solids or liquids, be capable of 
easily igniting macroscopic nuclear fusion reactions. Among 
a number of possible methods for the acceleration of pellets 
to 1000 km per sec, the construction of vastly improved 
shaped charges appears to be the most promising. 

Success may be predicted if developments are pursued along 


the following lines: 


1 Improved design of the shaped charges. 

2 Use of coruscative, self-reacting inserts (2). 

3 The use of more powerful explosives, especially those 
of the type of frozen radical or of metastable states. 

4 Increased velocity of detonation of the explosives. 
Recent work indicates that all of these developments are 
feasible, and that, if compounded, they will eventually 
achieve the desired results. Success with nuclear fusion 
would greatly facilitate explorations and constructions along 
the lines described. 


Tests of the surface of the moon 

The spectroscopic observation of flashes of sufficiently large 
and fast natural or artificial meteors impacting on the moon, 
as well as of the resulting dust clouds, would provide one of 
the most direct means for the analyses of the physico-chemical 
structure of the moon’s surface layers. Natural meteors are 
unfortunately not predictable, and consequently the neces- 
sary large telescope cannot be harnessed to the task. I have 
therefore attempted, first unsuccessfully from a V-2 rocket 
on Dee. 17, 1946, to launch artificial meteors into interplane- 
tary space, onto the moon and to the planets. Belated 
success came on Oct. 16, 1957, when J. Cuneo and I succeeded 
in launching a red hot pellet, which consisted mostly of the 
oxides of titanium and aluminum, at a speed of 15 km per 
sec, from an Aerobee rocket at an altitude of about 80 km 
(3). This pellet, which was about 1 cm in size, is the first 
manmade object to have been propelled permanently away 
from Earth. Proposals to launch larger particles to the 
moon with the aid of large shaped charges have not yet 
found any sponsor. 

In this connection it may also be asked why the impact 
of the Russian rocket on the moon was not observed. The 
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answer to this question is, in part, that the event could not be 
viewed or photographed from the Palomar Mountain Ob- 
servatory, since the moon at that time was below the eastern 


horizon. On the other hand, I know of no answer as to why 
the Russians, the Europeans and the American astronomers in 
the East did not succeed in recording the impact of Lunik IT 
on the moon. In any event, observations of future impacts 
on the moon are badly needed if we are to work out in detail 


some of the possibilities described in this note. 
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namic of Re-Entry 


Vehicle’ 
KENNETH WANG? and LU TING? 
Aerodynamics Laboratory, Polytechnic Institute of 
Brooklyn, Freeport, N. Y. 


HE PROBLEM of aerodynamic heating of space vehicles 

during re-entry has been investigated by Lees (1),4 Chap- 
man (2) and others. The heating rate at the stagnation point, 
in the limiting case of thermodynamic equilibrium, for lami- 
nar flow at hypersonic speed was obtained in (1). 

In a previous note (3), the fact that the principal part of 
the change in velocity took place in the part of the trajectory 
with higher atmospheric density was utilized to derive the 
analytic expressions for the velocity and the angle of inclina- 
tion. In the present note, it is found that the maximum 
heating rate takes place at the point of the trajectory very 
near to the minimum elevation. Again, the fact that the 
major portion of the heat input also takes place in the same 
part of trajectory as the velocity is used to obtain the analytic 
expression for the total heat input. 

The numerical results for the total heat input check closely 
with the machine calculations. 

The notations used in this note are the same as those in 
(4) except where noted. 

The stagnation point heating rate may be written in the 


form (1) 
Qo p'2V3/Ro (1] 
where 
K = constant 
Ro = radius of curvature of the surface at the stagnation 


point 
With the relation between the velocity and the density in 
(3), Equation [1] gives the heating rate Qo at the stagnation 
point as a function of density or elevation. Fig. 1 shows the 
heating rate vs. the logarithm of atmospheric density ratio 
for a particular trajectory. 


The condition for maximum heating rate dQ/dt = 0 yields 

which can be transformed into a similar relation as in (3) 
= (3CpA/mB) [2] 


Equation [2] together with Equation [9] in (3) determine 
the value of 3 and p at the maximum heating rate. Since 
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only a factor of 3 differentiates pai (2) from the similar 
expression in Equation [11] in (3), it can also be concluded 
that the maximum he: ating rate could be approximated by 
the heating rate at the minimum elevation. 

The total heat input at the stagnation point from the entry 
to the minimum elevation is, by integration of Equation [1] 


¢ 0 K dt 
Qo aa fi = oad. 3 [3] 


From the expression for the exponential atmosphere p = ps 
exp (—h), the following relation is obtained by differenti- 
ation with respect to time, for shallow trajectories 


V'dp 


K 


Fig. 1 indicates that the major portion of the heat input 
takes place in the part of the trajectory with higher atmos- 
pheric density, or near the point of minimum elevation. 
Therefore, a simple approximate relation of the form 


(V — Vo)? = C(o — p) [5] 


is used in the integration of Equation [4]. Note that Equa- 
tion [5] satisfies the values of the velocity and the slope dp/dV 
at the minimum elevation. In addition, the constant C is 


Qo = 
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Fig. 1 Heating rate vs. logarithm of atmospheric density 
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Equation [3] becomes 
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defined such that the curvature or 02p/0V? of the exact ve- 
locity and density relation at the minimum elevation is not 
altered. Thus, the constant Cisexpressedas == 


_ (SoA)? _ 
(24) (a —2b) 


Together with the analytic expression for the density and 
the angle of inclination obtained in (3), the integration of 
Equation [4] yields the total heat absorbed at the stagnation 


The corresponding expression for the total heat input when 
the entry velocity is less than the circular satellite velocity 
may be obtained in a similar manner, and will be reported in 
a forthcoming PIBAL report. 

Equation [6] may be modified for the evaluation of the heat 
input from the minimum elevation to the exit by replacing V, 
with V., which, as shown in (4), can be approximated by 
Vo?/V.. Also, the sign before the second term in the braces 
should be changed to minus. 


point from the entry to the point of minimum elevation 


K QV (can po? 
= — | 
BV Ry — 35 \ ing 


where 


= sin“) 


lic + 


(ton~ — 1 -— tan 2( 


Equation [6], and they deviate less than 3 per cent from the 
values calculated by numerical integration of the complete 
differential equations of motion. 


As a confirmation, numerical results were computed using 


AC pA 
mB Va — % 


Va — 3b 
mB 
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k = V—b/(a — 3b) 
F,E = elliptic integrals of the first and second kind, 3 


respectively 


In obtaining Equation [6] the entry velocity is taken to be 
greater than the circular satellite velocity (= 26,000 fps). 
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Wang, K. and Ting, L., ‘‘An Approximate Analytic Solution of Re- 


Entry Trajectory With Aerodynamic Forces,’’ ARS JourNat, vol. 30, no. 
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Ali Bulent Cambel, Northwestern University, Associate Editor 


Mechanical Properties of Intermetallic 
Compounds, edited by J. H. Westbrook, 
John Wiley and Sons, Inc., New York, 
1960, 435 pp. $9.50. 

Reviewed by D. H. WurtrmMorE 
Northwestern University 


This book is based on a symposium 
which was held in Philadelphia in May 
1959, under the sponsorship of the Electro- 
chemical Society. It represents a fairly 
extensive coverage of current research in 
the exciting field of mechanical behavior of 
that broad family of materials, the inter- 
metallic compounds. The subject matter 
of the book may be categorized as follows: 
(1) a comprehensive literature survey; (2) 
two chapters concerned exclusively with 
the effects of temperature and crystal 
structure; (3) three papers which treat the 
experimental methods that have proved 
useful in studying mechanical properties 
of these materials; (4) four chapters dis- 
cussing the role played by line and point 
imperfections, and (5) seven contributions 
dealing with phenomenological aspects of 
the mechanical behavior of particular 
compounds. 

Since space restrictions will not permit 
reviewing each of the 17 contributions 
in this volume, each of which is an 
entity in itself, the most satisfactory de- 
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scription of the book is provided by a 
résumé of its contents. The fundamental 
problems, as well as the complexities of the 
mechanical behavior of intermetallic com- 
pounds, have been extensively reviewed 
by J. H. Westbrook. A discussion of the 
relation of mechanical strength of inter- 
metallics to the electronic factor is pre- 
sented by G. M. Schwab and the effect of 
temperature on their mechanical proper- 
ties is elucidated by E. M. Savitskii who 
emphasizes some of the less well-known 
Russian researches in this area. Some of 
the useful techniques for elevated tempera- 
ture tensile and torsional tests are dis- 
cussed by L. Green, M. L. Stehsel and C. 
E. Waller. Results of fractrographic stud- 
ies of some of the nickel aluminides are 
summarized by R. W. Guard and A. M. 
Turkals, and recent success in the extru- 
sion of intermetallics is reported by D. L. 
Wood. N. Brown reviews some of the ef- 
fects of dislocations in alloys exhibiting 
long-range order and observations on the 
deformational modes in particular face- 
centered cubic superlattice alloys are dis- 
cussed by T. Tavka, K. Yasukochi and R. 
Honda. Use of the electron microscope in 
observing domain structure of the super- 
lattice in Cu-Au is reported by D. W. 
Pashley and A. E. B. Presland; radiation 
effects on the properties of titanium alu- 


minide are described by E. M. Grala and 
J. B. McAndrew. The plastic propertie 
of the semiconducting compound In-Sb 
are discussed by J. J. Duga, and work on 
elastic constants, anelastic effects and 
piezoresistivity is reported by R. F. Potter, 
J. H. Wasilik and R. B. Flippen. Me- 
chanical properties of transition metal 
beryllides are discussed by A. J. Stone- 
house, R. M. Paine and W. W. Beaver; 
the plastic behavior of the y2 phase of the 
Cu-Al system is described by D. J. Mack, 
A. J. Birkle and W. L. Krubsack, and the 
high temperature properties of some inter- 
metallic compounds of the transition met- 
als are reviewed by R. D. Grinthal. Re- 
cent Russian work on the mechanical 
properties of some intermetallic com- 
pounds of iron and nickel is discussed by 
I. I. Kornilov, and E. M. Grala describes 
his investigations of the ductility, hard- 
ness and strength of some nickel alu- 
minides. 

The two stimuli for the increased re- 
search effort in the field of intermetallic 
compounds have resulted from the never- 
ceasing need for outstanding, high tempera- 
ture materials to meet complex military and 
industrial requirements and the discovery 
of semiconducting behavior, particularly in 
the III-V intermetallic compounds. These 
materials, however, suffer from a lack of 
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sufficient ductility at ambient tempera- 
tures, and this will severely limit their com- 
mercial or military applications. Basic re- 
search efforts of the kind described in 
this book should eventually overcome this 
handicap. Those readers interested in 
keeping abreast of such current develop- 
ments in intermetallics should find this 
book most interesting. 


Evaporation and Droplet Growth in 
Gaseous Media, by N. A. Fuchs, Perga- 
mon Press, Ltd., London, 1959, 72 pp. 
$5.50. 

Reviewed by S. S. PENNER 
California Institute of Technology 


This little, expensive book is a transla- 
tion by J. M. Pratt of a monograph writ- 
ten by Fuchs several years ago, and con- 
tains no references to articles written after 
1954. Because of the unfortunate timing 
of the original publication, Fuchs’ mono- 
graph does not contain an account of the 
voluminous and informative work on drop- 
let and spray burning performed par- 
ticularly in the United States and in the 
United Kingdom during the last 6 years. 
Thus the reader interested in droplet and 
spray combustion phenomena, as well as in 
evaporation and condensation processes, 
will have to regard Fuchs’s monograph as 
an introductory treatise that must be sup- 
plemented by other reference sources, 


e.g., the Proceedings of the International 
Combustion Symposia, Proceedings of the 
AGARD Combustion and Propulsion 
Colloquia, recent texts on combustion 
processes. 

Fuchs treats quasi-stationary evapora- 
tion and growth of droplets in droplets that 
are stationary with respect to the medium 
(31 pages) and in moving droplets (22 
pages), also nonstationary evaporation and 
growth of droplets (8 pages). The trans- 
lation of the text is good and the printing 
legible although not up to the usual high 
standards of Pergamon Press. The an- 
alytical approach used by Fuchs is in- 
formative and demonstrates that the au- 
thor possesses considerable insight into 
the physical phenomena he discusses. 
The analyst who likes to derive all of his 
results on the basis of an appropriate set of 
conservation equations with suitable 
boundary conditions will find Fuchs’ treat- 
ment stimulating and conducive to a great 
deal of additional work. 


Nuclear Engineering Handbook, edited 
by Harold Etherington, McGraw-Hill 
Book Co., New York, 1958, xv + 1872 
pp. $25. 

Reviewed by J. A. SCANLAN 
The University of Texas 


This monumental work covers, in a 
single authoritative volume, practically all 


phases of the present theory and practice 
of nuclear engineering. Seventy special- 
ists contributed in their particular fields. 
The discussions, of course, are more con- 
densed than in texts, but the total coverage 
gives more for one’s money, even at $25, 
than any corresponding collection of books 
known to the reviewer. A comprehensive 
index, 52 pages, and well-organized con- 
tents makes it easy to find any topic. The 
format is compact, 6 X 9 X 2.25 in., yet 
durable and easy to use and read. Many 
figures and tables are included. The 14 
major sections are: Mathematical Data 
and General Tables, Nuclear Data, Math- 
ematics, Nuclear Physics, Experimental 
Techniques, Reactor Physics, Radiation 
and Radiological Protection, Control of 
Reactors, Fluid and Heat Flow, Reactor 
Materials, Chemistry and Chemical En- 
gineering, Nuclear Power-Plant Selection, 
Mechanical Design and Operation of Re- 
actors, and Isotopes. 

Extensive references and bibliographies 
are included at the end of each subsection. 
A random but rather extensive check by 
the reviewer of tabular data turned up no 
errors. 

Every engineer whose work even touches 
upon nuclear problems would do well to 
own this handbook. No library should be 
without it. The price is rather high for a 
casual university student, but the serious 
student would find it well worth the cost. 


Jet driven helicopter rotor (2,931,441). 
J. R. Root, Independence, Mo., assignor 
to Westinghouse Electric Corp. (ARS cor- 
porate member). 

Rotor arm with detachable tip for per- 
mitting access to an engine within a 
housing. 

Solar engine balloon altitude control 
(2,931,597). C. B. Moore Jr., Bedford, 
Mass., assignor to General Mills, Inc. 

Means for dropping ballast to com- 

pensate for loss of lift due to cooling of 
gas which occurs when the sun sets, and 
for terminating ballast dropping when 
descent becomes less than a predetermined 
rate. 
Ejection seat for aircraft (2,931,598). 
G. E. Sanctuary, Wichita, Kan., assignor 
to Boeing Airplane Co. (ARS corporate 
member). 

Rocket mounted on a seat structure to 
produce a jet thrust in a different direction 
than the seat, and causing the seat to tilt 
at a predetermined angle relative to the 
direction of ejection. 

Automatic star tracker (2,931,910). R. 
H. Ostergren and L. L. Nalley, Compton, 


Eprtor’s Norte: Patents listed above 
were selected from the Official Gazette of 
the U. S. Patent Office. Printed copies 
of patents may be obtained from the 
Commissioner of Patents, Washington 25, 
D. C., at a cost of 25 cents each; design 
patents, 10 cents. 
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Calif., assignor to Northrop Corp. (ARS 
corporate member). 

Optical system for focusing light from a 

star and its adjacent field. Scanning 
means for modulating the light into 
periodic light signals in accordance with 
a reference frequency. 
Method of initiating combustion of liquid 
hydrocarbon fuels (2,932,158). W. P. 
Ter Horst, McDonough, Md., assignor to 
Olin Mathieson Chemical Corp. (ARS 
corporate member). 

Introducing an azine into the combus- 
tion chamber and contacting it with 
fuming nitric acid, instantly initiating 
combustion. 

Ignition means for rocket (2,932,162). 
E. A. Malick (ARS member), Waco, 
Texas, assignor to Phillips Petroleum Co. 

Igniter for use with a propellant charge 
having exposed burning surfaces. Igniter 
members adapted to pierce the external 
surface of the charge. 

Quick disconnect nose (2,932,252). B. 
B. Korn, Los Angeles, Calif., assignor to 

— Aircraft Co. (ARS corporate mem- 
er). 

Locking rings rotatable about the axis 
of a guided missile composed of separable 
longitudinal sections. Detents hold rings 
in locked and unlocked positions. 

Ballistic missiles (2,932,467). G. C. 
Scorgie, Reading, England, assignor to The 
English Electric Co. 

_ Missile carrying a record of data relat- 
ing to a predicted course and a gravity- 
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influenced acceleration-responsive naviga- 
tion means for deriving data relating to 
its actual course. 

VTOL aircraft (2,932,468). P.G. Kappus 
(ARS member), Cincinnati, Ohio, as- 
signor to General Electric Co. (ARS cor- 
porate member). 

Air duct extending longitudinally 
through the body, having a forwardly 
directed inlet and a downwardly and 
rearwardly directed outlet. 

Hypergolic fuels (2,932,941). A. L. Ayers 
and C. R. Scott, Bartlesville, Okla., 
assignors to Phillips Petroleum Co. 

Method of developing thrust by the 
combustion of bipropellant components 
using sulfenamides and thiosulfenamides 
as fuels. 
et fuel compositions (2,932,942). G. G. 

ke and A. J. Kolka, Clairton, Pa., as- 
signors to Ethyl Corp. 

Process for cooling the lubricating oil 
in a jet engine, using a thermally stabi- 
lized fuel consisting of a distilled hydro- 
— fuel having an endpoint of 480 
eg. 
propulsion device (2,932,943). Dr. 

. Zwicky (ARS member), assignor to 
Aerojet-General Corp. (ARS corporate 
member). 

Hydroresonator comprising a duct with 
inlet and outlet openings, and a reaction 
chamber arranged so water can flow 
around the chamber containing a hydro- 
fuel spontaneously reactive with water. 
Rocket grain (2,933,041). J. M. Ambrose 
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signor to Phillips Petroleum Co. 

Cylindrical grain stable to thermal 
cycling provided with a single longitudi- 
nally extending spiral perforation of sub- 
stantial diameter. 

Annular wing flying machines (2,933,266). 
H. P. von Zborowski, St. Antoine par 
Brunoy, France. 

Ramjet powerplant adapted to exert 
thrust along a given line, and annular 
lifting means forming the aerodynamic 
sustaining means as a_ tunnel-shaped 
wing. 


Cooling system for rocket engines 
(2,933,888). A. Africano and C. W. 
Chillson (ARS members), Packanack 


Lake, N. J., assignors to the U. S. Air 
Force. 

Members extending axially along a 
nozzle and spaced from each _ other. 
Wire is coiled around the nozzle and 
members to form a coolant passage be- 
tween. 


Rocket motor thrust cutoff (2,933,889). 
R Tolkmitt, China Lake, Calif., 
assignor to the U.S. Navy. 

Pie-shaped cutouts in the forward face 
of a plate-like member rotatable to 
another position to align with a nozzle 
and terminate thrust while permitting a 
solid propellant to continue burning. 


Automatic repeating rocket launcher 
(2,933,981). P. E. Anderson (ARS mem- 
ber), T. Q. Cromp, P. I. Evans (ARS 
member) and M. Ransom, Torrance, 
Calif., assignors to the U.S. Army. 

Single barrel launcher tube with in- 
terconnected rocket-carrying tubes suc- 
cessively feeding rockets from a storing 
magazine. 

Radio controlled missile (2,934,286). E. 
F. Kiernan, San Diego, Calif. 

Antenna array comprising a pair of 
half-wave doublet antennas, each mounted 
a half wave length from each other, and 
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(ARS member), McGregor Texas, as- 


connected to a rudder control. 
Sonde (2,934,287). J. W. Ault, Silver 
Spring, Md., assignor to the U.S. Navy. 

Means responsive to varying conditions 

of light entering windows in a spinning 
missile for generating modulation im- 
pulses. The sonde radiates wave signals 
indicating light conditions encountered 
in flight. 
Marine cargo vessel design (188,013). 
C. A. Gongwer (ARS member) and F. J. 
Walcott Jr., Glendora, Calif., assignors to 
Aerojet-General Corp. (ARS corporate 
member). 

Pontoon-type passenger and crew cabin 
elevated by means of a central vertical 
fin projecting high above a cargo sub- 
marine of high aspect ratio. 
Telescoping ramjet (2,935,946). E. A. 
Gallo, D. B. Clark and E. W. Schwartz 
(ARS member), El Cajon, Calif., assignor 
to the U. S. Army. 

Projectile with center body portion 
telescoped within the combustion cham- 
ber when the chamber is not in use. 
When fired, the chamber shell is pushed 
rearward by aerodynamic forces, expand- 
ing the projectile to operating position. 
Three-axis gyroscopic aerodynamic damp- 
ing system (2,935,947). L. T. Jagiello 
(ARS member), China Lake, Calif., 
assignor to the U.S. Navy. 

Gyroscopes mounted in control surfaces 
hingedly connected to fins to provide cor- 
rection of pitch, roll and yaw of the mis- 
sile by providing damping at the hinged 
connection. 

Rocket motor throttling injector (2,936,- 
577). J. S. Amnéus, Woodside, Calif., 
assignor to the Univ. of Michigan. 

Movement of a throttle provides an 
oxidizer lead on starting the motor, 
and an oxidizer override on shutdown. 

Jet deflector (2,936,582). I. M. David- 
son, Farnborough, England, assignor to 
Power Jets (Research and Development). 


Canopy of flexible metal cloth and flex- 

ible cords of metal strands, packed and 
stowed in a housing and capable of open- 
ing up to form a chute deployed in the 
jet stream. 
High Mach deceleration device (2,936,- 
710). Dr. W. Boliay (ARS member), 
Pacific Palisades, Calif., assignor to 
Curtiss-Wright Corp. 


Thin shell surrounding a missile, tele- 

scoped to form a pocket blocking air flow 
between the nose cone and the shell, 
thereby spoiling the lift. 
Jet sustained and propelled aircraft 
(2,936,969). A.A. Griffith and W. Shaw, 
Alvaston, England, assignors to Rolls- 
Royce Ltd. 

Body of aerodynamic sustaining form 
carring horizontal propulsion engines, 
and jet-lift-type engines with their axes 
extending vertically through the body. 
Propeller sustained aircraft (2,936,972) 
J. Zinavage, Baltic, Conn. 

Conical lifting surface with its peri- 
meter spaced above the top of a vertical 
duct having an open top and bottom. 
Air is discharged downward in the duct. 
VTOL aircraft (2,936,973). P.G. Kappus 
(ARS member), Cincinnati, Ohio, assignor 
to General Electric Co. (ARS corporate 
member). 

Tubular housing closed at one end to 
form a load-carrying surface and having 
a gas generating engine exhausting through 
a nozzle at the other end. Controls vary 
the orientation of vanes carried within 
the nozzle. 43 
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Nigam and 8. N. Singh, Quart. J. Mech. & 
Appl. Math., vol. 13, Part 1, Feb. 1960, 
pp. 85-97. 
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Rasor and J. D. McClelland, Rev. Sci. 
Instr., vol. 31, June 1960, pp. 595-604. 
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Phys.: JETP, vol. 37(10), no. 6, June 
1960, pp. 1167-1178. 
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a Magnetic Field, by Yu. A. Tarasov, 
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tached Shock Waves, by Richard A. 
Scheuing, ARS Preprint 1112-60, May 
1960, 29 pp. 

Similarity Solutions for Inviscid Hyper- 
sonic Flow Over Slender Power Law and 
Related Bodies, by Harold Mirels, ARS 


M.H. Smith, Associate Editor 


The James Forrestal Research Center, Princeton University 


Preprint 1111-60, May 1960, 46 pp. 
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in Low Density Supersonic Flow, by S. A. 
Schaaf, E. S. Moulic, M. T. Chahine and 
G. J. Maslach, ARS Preprint 1130-60, 
May 1960, 23 pp. 

High Temperature Rarefied Ultra-high 
Mach Number Flow over a Flat Plate, by 
Henry T. Nagamatsu and R. E. Sheer Jr., 
ARS Preprint 1132-60, May 1960, 13 pp. 

Effect of Surface Mass-transfer on the 
Hypersonic Strong Interaction Between 
Shock Wave and Boundary Layer, by 
Ting Yi Li, ARS Preprint 1166-60, May 
1960, 20 pp. 

Dissociation Effects upon Compressible 
Turbulent Boundary Layer Skin Friction 
and Heat Transfer, by William H. 
Dorrance, ARS Preprint 1143-60, May 
1960, 42 pp. 
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Stanford Unw., Dept. Aeron Engng., 
SUDAER 92, March 1960, 11 pp. 
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Calculation of the Internal Temperature 
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Rise Due to Aerodynamic Heating, by 
Eugene E. Covert, MIT, Naval Super- 
sonic Lab., Tech. Rep. 202, Feb. 1957, 19 
pp. 

Nonsteady Discharge of Subcritical 
Flow, by George Rudinger, Princeton 
Univ., Project Squid, Tech. Rep. CAL-80- 
F, May 1960, 44 pp. (Available only on 
microcard.) 

General Research in Flight Sciences; 
Jan. 1959-Jan 1960, Vol. 1, Part 1, 
Fluid Mechanics, Lockheed Aircraft Corp., 
Missile & Space Div., Rep. LMSD-288,- 
139, vol. 1, Part 1, Jan. 1960. 

The Effect of Body Forces in Melting 
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D. M. Tellep, Paper 1, 13 pp. 

Radiant Energy Transfer in Gaseous 

Flows, by D. M. Tellep and D. K. 

Edwards, Paper 2, 31 pp. 

A Numerical Method of Solution for 

Heat Conduction in Composite Slabs 

With a Receding Surface, by J. J. 

Brogan, Paper 3, 37 pp. 

On Mass Transfer and Shock-generated 

Vorticity, by H. Hoshizaki, Paper 4, 

34 pp. 

The Effect of Shock-generated Vorticity, 

Surface Slip, and Temperature Jump 

on Stagnation-point Heat-transfer 

Rates, by H. Hoshizaki, Paper 5, 7 pp. 

Hypersonic Flow Around Bodies of 

Revolution Which Are Generated by 

Conic Sections, by M. Vinokur, Paper 

6, 21 pp. 

Kinematic Formulation in Rotational 

Gas Flow, by M. Vinokur, Paper 7, 19 

pp. 

Hypersonic Flow Around Blunt Bodies 

of Revolution Whose Shock Waves Are 

Generated by Conic Sections, by M. 

Vinokur and R. W. Sanders, Paper 8, 

22 pp. 

Inviscid Hypersonic Flow Around 

Spheres With Mass Injection, by M. 

Vinokur and R. W. Sanders, Paper 9, 14 

pp. 

Third-order Cylindrical Blast-wave 

Theory and Its Analogy to the Flow 

Around Hemisphere-cylinder Configura- 

tions in High-speed Flight, by R. J. 

Swigart, Paper 10, 41 pp. 

Heating Rate Characteristics of an 
Oxyacetylene Flame Apparatus for 
Surface Ablation Studies, by Edwin 
Joseph Russ, Calif. Univ., Inst. Engng. 
Res., Tech. Rep. HE-150-176, Dec. 21, 
1959, 68 pp. 

On the Magnetogasdynamics of Com- 
pressible Vortices, by James E. McCune 
and Coleman DuP. Donaldson, ARS 
Preprint 1139-60, May 1960, 30 pp. 

General Research in Flight Sciences, 
Jan. 1959-Jan. 1960, Vol. 1, Part 2 
Fluid Mechanics, Lockheed Aircraft Corp., 
Missiles & Space Div., Rep. LMSD-288,- 
139, vol. 1, Part 2, Jan. 1960. 

Shock-tube Experiments Simulating 

Planetary Atmospheres, by 

R. W. Rutowski and K. K. Chan, 

Paper 1, 22 pp. 

A Note Regarding Stagnation-point 

Radiation Heat Transfer, by R. W. 

Rutowski, Paper 2, 7 pp. 

Formation of Nitric Oxide Behind 

Strong Shock Waves in Air, by J. J. 

Allport, Paper 3, 19 pp. 

Steady Two-dimensional Flow with a 

Transverse Magnetic Field, by M. 

Mitchner, Paper 4, 19 pp. 

Some Reflected Shock Parameters in 

N, and O:2, by M. Scheibe, Paper 6, 15 

pp. 

Stagnation-point Heat Transfer 

Measurements in Hypersonic, Low- 
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density Flow, “4 S. E. Neice, R. W. 
Rutowski, and K. K. Chan, Paper 7, 
17 pp. 

Comparative Studies of Vibrational 

Relaxation in Nitric Oxide, by P. R. 

Monson, J. J. Allport, and F. A. 

Robben, Paper 8, 19 pp. 

Heat Transfer and Skin Friction in 
Magnetogasdynamic Channel Flow, by 
F. D. Hains, ARS Preprint 1138-60, 
May 1960, 40 pp. 

Measured Pressure Distribution and 
Local Heat Transfer Rates for Flow Over 
Concave Hemispheres, by J. C. Y. Koh 
and J. P. Hartnett, ARS Preprint 1166- 
60, May 1960, 12 pp. 

The Aerothermodynamics of Reentry 
Trails, by Madeleine Goulard and Robert 
Goulard, A RS Preprint 1145-60, May 1960 
17 pp. 

Real Gas Flow Tables for Nondissoci- 
ated Air, by Louis G. Kaufman II, Wright 
Air Dev. Center, Tech. Rep. 59-4, Jan. 1959, 
258 pp. 

Plasma Oscillations, I, by I. B. Bern- 
stein and 8. K. Trehan, Princeton Univ., 
Project Matterhorn, Rep. Matt-42, May 
1960, 163 pp. 

Heating of an Ionized Gas Sheath by 
Microwaves, by Mahendra Singh Sodha, 
Appl. Sci. Res., vol. 8B, no. 3, 1960, pp. 
208-212. 


Uniqueness of Radioactive Solutions to 
the Schroedinger Wave Equation, by C. 
Zemach and F. Odeh, Arch. for Rational 
Mech. & Anal., vol. 5, no. 3, May 17, 1960, 
pp. 226-237. ; 

Nonlinear Parabolic Equations with 
Application to Flow of Gases Through 
Porous Media, by A. Friedman, Arch. 
for Rational Mech. & Anal., vol. 5, no. 3, 
May 17, 1960, pp. 238-248. 

Taylor Instability in Shock Acceleration 
of Compressible Fluids, by Robert D. 
Richtmyer, Communications on Pure & 
Appl. Math., vol. 13, no. 2, May 1960, 
pp. 297-320. 

Some Problems Associated with the 
Definition of the Heat of Transfer for 
Binary Liquid Systems, by H. J. V 
Tyrrell, Trans., Faraday Soc., vol. 56, no. 
6, June 1960, pp. 770-787. 

Vaporization Processes in the Hyper- 
sonic Laminar Boundary Layer, by 
Sinclaire M. Scala and Guido L. Vidale, 
Internat. J. of Heat & Mass Transfer, 
vol. 1, no. 1, June 1960, pp. 4-22. 

Experimental Methods Applied to the 
Determination of Some Temperature 
Radiation Parameters, by D. T. Kokorev, 
Internat. J. of Heat & Mass Transfer, vol. 
1, no. 1, June 1960, 99. 23-27. 

Thermal Radiation Between Parallel 
Plates Separated by an Absorbing-emit- 
ting Nonisothermal Gas, by C. M. 
Usiskin and E. M. Sparrow, Internat. 
J. of Heat & Mass Transfer, vol. 1, no. 1, 
June 1960, pp. 28-36. 


Hypersonic Research Summary, by G. 
J. Fabian, Cornell Aeron. Lab., Rep. 
Ad-1118-A-11, June 1960, 24 pp. 


Shock Wave and Flow Field Develo; 
ment in Hypersonic Re-entry, by Ronald 
F. Probstein, 4 RS Preprint 1110-60, May 
1960, 36 pp. 

Ballistic Research Laboratories’ New 
Hypersonic Tunnel, by J. Sternberg, 
Aberdeen Proving Ground, Ball. Res. Labs., 
BRL Rep. 1076, Jan. 1960, 24 pp. 

Reduction of the Boundary Layer 
Thickness in a Hypersonic Nozzle, by 
Serge A. Waiter, Univ. Southern Calif, 
Engng. Center, USCEC Rep. 56-213, May 
20, 1960, 22 pp. 

Effects of Fixing Boundary-layer Transi- 


tion for a Swept- and a Triangular-wing 
and Body Combination at Mach Numbers, 
by Louis S. Stivers Jr., NASA TN D- 
312, June 1960, 28 pp. 

On the Stability of a Heterogeneous 
Shear Layer Subject to a Body Force, by 
J. Menkes, CIT, Jet Prop. Lab., Tech. 
Rep. 32-12, March 25, 1960, 8 pp. 

High Temperature Heat Transfer from 
Gases to Cylinders and Nozzles, by I. E. 
Kanter, F. Martinek, M. L. Ghai, D. 
Namkoong and J. Loch, General Electric 
Co., Flight Prop. Lab. Dept., Feb. 1, 1958- 
Feb. 28, 1959, 128 pp. 

Spectroscopic Technique for Tempera- 
ture-density Measurements in Oxygen- 
bearing Flows, by W. H. Wurster and C. 
E. Treanor, Cornell Aeron. Lab., Rep. 
AD-1118-A-10, Dec. 1959, 40 pp. 

Some Unexpected Results of Shock 
Heating Xenon, by Per Gloersen, General 
Electric Co., Missile & Space Vehicle 
Dept., Rep. TIS R60SD364, April 26. 
1960, 53 pp. 

An Integral Method for Natural- 
convection Flows at High and Low Prandtl 
Numbers, by Willis H. Braun and John E. 
Heighway, N ASA TN D-292, June 1960. 
40 pp. 

Effect of Buoyancy Forces on Mass 
Transfer Cooling, by R. Eichhorn, Univ. 
Minnesota, Inst. Tech., Heat Transfer 
Lab., Tech. Rep. 25, Feb. 1960, 8 pp. 

Non-equilibrium Centered Rarefaction 
for a Reacting Mixture, by Luigi G. 
Napolitano, Arnold Engng. Dev. Center, 
TN, 60-129, June 1960, 63 pp. 

Reaction Rate Distributions on Catalytic 
Surfaces—II. Comparison of Exact and 
Approximate Flat Plate Solutions for 
Asymptotically Small Prandtl Number, by 
Daniel E. Rosner, Aero-Chem Res. Labs., 
Tech. Pub. 16, April 1960, 13 pp. 

Experimental Investigations of the 
Interactions Between a Shock and a 
Magnetic Field, by John Paul Barach, 
Maryland Univ., Inst. Fluid Dynam. & 
Appl. Math., TN BN-203, March 1960, 
29 pp. 

Behavior of Solutions of the Navier- 
Stokes Equations for a Complete Class of 
Three-dimensional Viscous Vortices, by 
Coleman duP. Donaldson and Roger D. 
Sullivan, Aeron. Res. Associates of Prince- 
ton (ARAP), May 1960, 22 pp. 

A Discussion of Two Incompressible 
Magnetohydrodynamic Flows, by Keith 
P. Kerney, Wright Air Dev. Div., Tech. 
Rep. 59-486, Part 3, Oct. 1959, 52 pp. 

On Fully Developed Channel Flows: 
Some Solutions and Limitations, and 
Effects of Compressibility, Variable Prop- 
erties, and Body Forces, by Stephen H. 
Maslen, NASA Tech. Rep. R-34, 1959, 
22 pp. 

Equilibrium Emissivity Calculations 
(from Spectroscopic Data) for CO, 
at 300°K and at 600°K, by M. Lapp and 
S. S. Penner, Calif. Inst. Tech., Daniel & 
Florence Guggenheim Jet Prop. Center, 
Tech. Rep. 32, June 1960, 18 pp. 

The Equation of State of an Ionized Gas, 
by George K. Bienkowski, Wright Air 
Dev. Div., Tech. Rep. 59-486, Part 2, 
Sept. 1959, 62 pp. 

The Wave Motions of Small Amplitude 
in a Fully Ionized Plasma, Part II—With 
Longitudinal Applied Magnetic Field, by 
S. I. Pai Unw. Maryland, Inst. Fluid 
Dynam. & Appl. Math., TN BN-207, 
April 1960, 28 pp. 

Equilibrium Emissivity Calculations for 
a Hydrogen Plasma at Temperatures up to 
10,000°K, by D. Olfe, Calif. Inst. Tech., 
Guggenheim Jet Prop. Center, Tech. Rep. 
33, ne 1960, 67 pp. 
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PENNSALT is one of America’s prime producers of 
AMMONIUM PERCHLORATE, the leading oxidizer used 
in solid propellants. Recent expansion of Pennsalt’s 
production facilities assures a reliable source of a 
consistently uniform material. 


SEND FoR complete specifications of PENNSALT AM- 
MONIUM PERCHLORATE. Produced as either a condi- 
tioned or non-conditioned material. Write PENNSALT 
CHEMICALS CORPORATION, Tacoma 1, Washington. 


PLANTS THROUGHOUT THE U. S., WITH INTERNATIONAL REPRESENTATION 


Pennsalt is making its contribution to the space age 
and to modern rocketry, just as this pioneer firm has 


Industrial Chemicals Division -West 


TACOMA 1, WASHINGTON 
Plants and Offices: Los Angeles and Fresno, Calif. 
Philadelphia, Pa. * Portland, Oregon * Vancouver, B.C. 


Other Pennsalt Chemicals Corporation plants and offices at Apple- 
ton, Atlanta, Aurora, Bryan, Calvert City, Chicago, Dallas, Detroit, 
Houston, Montgomery, New York, Pittsburgh, St. Louis, Wyandotte. 
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met and solved hundreds of specialized problems since 
its founding in 1850. 

Pennsalt plants and development centers are located 
throughout the United States, and service represen- 
tatives are located in several industrial areas abroad. 


PERCHLORYL FLUORIDE 
a new rocket fuel oxidizer—Pennsalt has cross- 
bred fluorine with two other very active elements— 
oxygen and chlorine—to produce the new rocket 
fuel oxidizer—perchlory] fluoride—a promising new 
development in rocketry. 


ats PROVIDING INDUSTRY WITH QUALITY CHEMICALS 


Foreign plants and 
offices: Pennsalt of Can- 
ada Ltd., Oaksville, On- 
tario. Pennsalt Chemi- 
cals GmbH, Ludenscheid, 
Germany. Industrial 
Quimica Pennsalt, Mex- 
ico City, Mexico. 
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Robert M. Wood and Ronald J. Tagliani, 
Aero/Space Engng., vol. 19, July 1960, 
pp. 32-38. 

The Variational Formulation of the 
Magneto-hydrostatic Equations, by P. C. 
Kendall, Astrophys. J., vol. 131, no. 3, 
May 1960, pp. 681-683. 

Principle of Corresponding States for 
Transport Properties, by Eugene Helfand 
and Stuart A. Rice, J. Chem. Phys., vol. 
32, June 1960, pp. 1642-1643. 

Stability of a Broken-line Jet in a 
Parallel Magnetic Field, by Philip G. 
Drazin, J. Math. & Phys., vol. 39, April 
1960, pp. 49-53. 

Compressibility Effects in Magneto- 
aerodynamic Flows Past Thin Bodies, by 
J. E. McCune and E. L. Resler Jr., J. 
Aero/Space Sci., vol. 27, July 1960, pp. 
493-503. 

An Experimental and Theoretical Study 
of Quartz Ablation at the Stagnation 
Point, by Mac C. Adams, William E. 
Powers and Steven Georgiev, J. Aero/- 
— Sci., vol. 27, July 1960, pp. 535- 
43. 

Radiation Cooling of Finite Heat- 
conducting Solids, by P. J. Schneider, J. 
Aero/Space Sci., vol. 27, July 1960, pp. 
546-548. 

A New Concept for Calculating Thermal 
Boundary Layers, by Ernest Wilhelm 
Adams, J. Aero/Space Sci., vol. 27, July 
1960, pp. 549-550. 

Viscous Hypersonic Similitude for a 
Dissociating Gas, by Robert J. Whalen, 
J. Aero/Space Sct., vol. 27, July 1960, pp. 
550-551. 

Hydromagnetic Effects on Heating and 
Shear at a Three-dimensional Stagnation 
Point in Hypersonic Flow, by Nelson H. 
Kemp, J. Aero/Space Sci., vol. 27, July 
1960, pp. 553-554. 

Inviscid Hypersonic Flow over Un- 
yawed Circular Cones, by L. Pottsepp, J. 
Aero/Space Sci., vol. 27, July 1960, pp. 


558-559. 


Orbit Determination from Range and 
Range-rate Data, by Robert M. L. 
Baker Jr., ARS Preprint, 1220-60, May 
1960, 48 pp. 

General Research in Flight Sciences; 
Jan. 1959-Jan. 1960, Vol. 3, Flight Dy- 
namics and Space Mechanics, Lockheed 
Aircraft Corp., Missiles & Space Div., 
Rep. LMSD-288, 139, vol. 3, Jan. 1960. 

Topological Structure of the Space of 

Trajectories for the Restricted Problem 

of Three Bodies, by Sze-Tsen Hu, Paper 

1, 21 pp. 

Minimum Transfer Time for a Power- 

limited Rocket, by G. Leitmann, Paper 

2, 32 pp. 

On the Anomalies in the Movements of 

the Perigee Due to Oblateness and Air 

Drag for Critical Values of Inclination, 

by J. V. Breakwell, Paper 4, 7 pp. 

The Local Accuracy of Doppler Orbit 

Determination, by J. V. Breakwell, 

Paper 5, 6 pp. 

On the Motion of a Satellite About a 

Planet With Rotational Symmetry, by 

C. M. Petty, Paper 6, 39 pp. 

Researches in Interplanetary Transfer, 

by J. V. Breakwell, R. W. Gillespie and 

S. Ross, Paper 7, 38 pp. 

Mechanics Problems of Spaceflight, 

by W. E. Jahsman, G Cline Jr., 

H. V. Hahne and W. Nachbar, Paper 

9, 89 pp. 
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Protection by- Ablation, by 


Spin-induced Forced Resonant Be- 
havior of a Ballistic Body Re-entering 
the Atmosphere, by J. S. Kanno, 
Paper 10, 28 pp. 

Spin-induced Forced Resonant Be- 

havior of a Ballistic Body Having Stable 

Nonlinear Aerodynamic Properties, by 

J.S. Kanno, Paper 11, 46 pp. 

Optimization of Trajectories, by J. V. 

Breakwell, Paper 12, 32 pp. 

Evaluation of the Motion and Aerody- 

namic Characteristics of High-fineness- 

ratio Reentry Configurations, by R. L. 

Nelson, Paper 13, 17 pp. 

Meteroids Versus Space Vehicles, by 
R. L. Bjork, ARS Preprint 1200-60, May 
1960, 5 pp. 

Motion of an Orbiting Vehicle About its 
Center of Mass Due to the Gravity 
Gradient, by J. P. King Jr., ARS Pre- 
print 1222-60, May 1960, 9 pp. 

Gradient Theory of Optimal Flight 
Paths, by Henry J. Kelley, ARS Pre- 
print 1230-60, May 1960, 13 pp. 

Planar Motions About an Oblate 
Planet, by Maurice L. Anthony and 
George Fosdick, ARS Preprint 1231- 
60, May 1960, 22 pp. 

The Choice of Unperturbed Orbit in 
the Use of Encke’s Method for the 
Effects of Oblatness and Drag, by Nor- 
man S. Hall, H. G. Galowicz and E. J. 
Wallman, ARS Preprint 1232-60, May 
1960, 13 pp. 

General Research in Flight Sciences; 
Jan. 1959-Jan. 1960, Vol. 1, Part 2, Fluid 
Mechanics, Lockheed Aircraft Corp., Mis- 
siles & Space Div., Rep. LMSD-288,139, 
vol. 1, part 2, Jan. 1960. 

Charge and Magnetic Field Inter- 
action with Satellites, by David B. 
Beard and F. S. Johnson, Paper 
17 pp. 

Attitude Drift of Space Vehicles, by 
William T. Thomson and Gordon S. 
Reiter, J. Astronaut. Sci., vol. 7, no. 2, 
Summer 1960, pp. 29-34. 

Normal Lunar Impact Analysis, by 
Anthony Liu, J. Astronaut. Sci., vol. 7, 
no. 2, Summer 1960, pp. 35-39. 

Error Analysis of Satellite Orbits in 
the Presence of Drag, by Frederick V. 
Pohle, J. Astronaut. Sci., vol. 7, no. 2, 
Summer 1960, pp. 40-44. 

The Intersection of Coplanar, Con- 
focal Conic Sections, by Geza S. Gedeon, 
J. Astronaut. Sci., vol. 7, no. 2, Summer 
1960, pp. 45-46. 

Two Maneuver Ascents to Circular 
Orbits, by J. F. Wolfe and D. DeBra, 
J. Astronaut. Sci., vol. 7, no. 2, Summer 
1960, pp. 47-48. 

Generalized Tables for the Calcula- 
tion of Trajectory Curves for Bodies 
Moving in Air, by M. M. Currie, Canada, 
Nat. Res. Council, LR-277, April 1960, 
10 pp. 

Boundary Layer Separation Effects on 
the Static Stability of a Flared-tailed 
Missile Configuration at / = 2 to 5, 
by J. Don Gray, Arnoid Engng. Dev. 
Center, AEDC-TN-103, June 1960, 27 pp. 

Effect of Eccentricity of the Lunar 
Orbit, Oblateness of the Earth, and Solar 
Gravitational Field on Lunar Trajec- 
tories, by William H. Michael Jr. and 
Robert H. Tolson, NASA TN D-227, 
June 1960, 37 pp. 

Analysis of Trajectory Parameters for 
Probe and Round-trip Missions to Mars, 
by James F. Dugan Jr., NASA TN D- 
281, June 1960, 68 pp. 

Motion of Low-thrust Vehicle in a 
Strong Gravitational Field, by W. A. 
Moser and J. B..Denniston, North Amer- 
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flames, and scientific interpretetion of 
results, mathematical analysis of rocket 
exterior ballistic problems, such as per- 
formance stability, etc. Good thecret- 
ical background in fluid dynamics and 
thermodynamics. in mechanicel, 
chemical. or aeronautical engineering 
with a thorough knowledge of all 
phases of heat transmission. 


Rocket Project Engineer 

To work on prototype development 
of solid propellant rocket propulsion 
systems. Involves coordination of pro- 
grams with mechanical design, ballistic, 
instrumentation and test engineers. 
Requires personeble men with interests 
and demonstrated abilities to carry sys- 
tems from inception through test and 
qualification. 

B.S. or M.S. in mechanical or chem- 
ical engineering with the ability to 
work effectively with others. 


LABORATORY DIRECTOR 


Hyperthermal Gasdynamics Laboratory 


DESCRIPTION OF POSITION RESPONSIBILITIES 

Initiate, plan, solicit, conduct and supervise advanced re- 
search activities, both theoretical and experimental. These 
activities conducted for government agencies and _ private 
corporations include: | 


To supervise assembly and loading 
of test rockets and customer units. 
Manage static firing fecility. including 
scrutiny of static test setups to insure 
proper alignment and support, of in- 
strumentation and calibration systems 
to insure data accuracy, and regular 
equipment checkouts. Cooperate with 
development project enaineers in test 
scheduling, review test results from 
computation group, and contribute to 


| Rocket Test Supervisor 


¢ Simulation of space vehicle re-entry phenomena in 

hyperthermal, hypervelocity plasma tunnels 

Investigation of the state of high temperature plasma 

in arc chambers, mixing chambers, and nozzles 

Design and development of instrumentation techniques | 

and systems applicable to the improvement of such in E, or Ch.E., aptitude in 

devices as pitot probes, calorimeters, total enthalpy mechenics and electronics, and con- 
siderable experience in production 

and mass flow probes and electronic measuring instrumenta- 

e Study of arc generator characteristics tion 

Communicate with customers, management andthe Research Aerodynamicist 
scientific community. Assume responsibility for conduct of To: work 
the Hyperthermal Gasdynamics Laboratory with specific re- Group on aerodynamic heating prob- 

gard for timely and satisfactory contract execution, personnel lems, gaseous combustion of turbulent 
alten facilities’ modernization, expansion, financing, and 
business administration. 

BACKGROUND EXPERIENCE REQUIRED 

Ph. D. in Aeronautical Engineering, Engineering Sciences, 
or equivalent with minor emphasis in Physical Chemistry, 
Physics and Mathematics. 

Minimum of five to ten years of theoretical and experi- 
mental research related to such phenomenological processes 
as abalation, gas dissociation and recombination, gas flow 
interaction with electric arc discharges, aerodynamic heating. 


Atlantic Research, the largest independent R&D or- 
ganization in the Washington, D.C. area, is active in 
all phases of solid propellant rocketry and propulsion 
systems, from basic research on combustion to the 
manufacture of shelf items. All our staff require- 
ments are for permanent positions offering the oppor- 
tunity to contribute to our expansion and your future. 


Specific background with reference to aerodynamic sim- US. citizenship is required for all pesitions. ' 
ulation concepts, similarity laws, basic physics related to ie 
ite . ice acic ac ic > 
plasma _ technology, orbital mechanics, basic gasdynamic 
Clarence H. Weissenstein, Director, 


principles, shock tube and hypersonic flow research. 
Please send resume to Pete De Paolo, Personnel Manager. 
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MOSPHERE, John V. Becker (TN).... 
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Gas TeMperatures, Hugh N. 
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TRANSIENT HEAT CONDUCTION IN A 
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So.uTion, R. W. Sanders (TN)..... 

ArropyNaMic HeEatTinG oF RE-ENTRY 
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ASCENDING SATELLITE VEHICLE FOR 
Crracucar Ornsirts, A. D. Cohen and H. 
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Masses, Howard D. Greyber (TN)... 
Optimum Time TO RENDEzvovus, John 
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FREE SURFACE CONDITION FOR SLOSHING 
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LaTeRAL HypRODYNAMIC FORCES ON AN 
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WATER, Theodore R. Goodman (TN). 
LAUNCH PARAMETERS FOR INTERPLANE- 
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Rocket, J. J. Rodden and R. J. Pollack 964 
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tems, O. M. Kryzhanovskii and V. M. 
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With ArsiTRARY PeRiop, INVARIANT 
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